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Methods and Apparatus for Det:ermining the Size 
and Shape of Particles 

Related Applications 

This applications claims priority of U.S. Provisional 
Patent Application No. 60/550^591, filed March 06, 2004 



TECHNICAL FIELD OF THE INVENTION 

In general, the present invention relates to systems and 
methods that analyze particles in a sample using laser light 
diffraction. More particularly, the present invention relates 
to systems and methods that analyze laser light diffraction 
patterns to determine the size of particles in a sample. 



Disclosure Of The Invention 



This disclosure describes an instrument for measuring the size 
distribution of a particle sample by counting and classifying 
particles into selected size ranges. The particle 
concentration is reduced to the level where the probability of 
measuring scattering from multiple particles at one time is 
reduced to an acceptable level. A light beam is focused or 
collimated through a sample cell, through which the particles 
flow. As each particle passes through the beam, it scatters, 
absorbs, and transmits different amounts of the light, 
depending upon the, particle ,size. So both the decrease in the 
beam intensity, due to light removal by the particle, and 
increase of light, scattered by the particle, may be used to 
determine the particle size, to classify the particle and 
count it in a certain size range. If all of the particles pass 
through a single beam, then many small particles must be 
counted for each large one because typical distributions are 
uniform on a particle volume basis, and the number 
distribution is related to the volume distribution by the 
particle diameter cubed. This large range of counts and the 
Poisson statistics of the counting process limit the size 
dynamic range for a single measurement. For example, a uniform 
particle volume vs. size distribution between 1 and 10 microns 
requires that one thousand 1 micron particles be measured for 
each 10 micron particle. The Poisson counting statistics 
require 10000 particles to be counted to obtain 1% 
reproducibility in the count. Hence one needs to measure more 
than 10 million particles. At the typical rate of 10,000 
particles per second, this would require more than 1000 
seconds for the measurement. In order to reduce the 
statistical count uncertainties, large counts of small 
particles must be measured for each large particle. This 
problem may be eliminated by flowing portions of the sample 
flow through light beams of various diameters, so that larger 
beams can count large count levels of large particles while 
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small diameter beams count smaller particles without the small 
particle coincidence counts of the large beam. The best 
results are obtained by using multiple beams of ever 
decreasing spot size to improve the dynamic range of the 
count. The count vs. size distributions from each beam are 
scaled to each other using overlapping size ranges between 
different pairs of beams in the group, and the count 
distributions from all of the beams are then combined. 



BEST MODES OF CARRYING OUT THE INVENTION 

Light scattered from the large diameter beam should be 
measured at low scattering angles to sense large particles. 
The optical pathlength of this beam in the particle sample 
must be large enough to pass the largest particle of interest 
for that beam. For small particles, the interaction volume in 
the beam must be reduced along all three spatial directions. 
The beam crossection is reduced by an aperture or by focusing 
the beam into the interaction volume. The interaction volume 
is the intersection of the incident light beam and the field 
of view for the detector which measures scattered light from 
the particle. However, for very small particles, reduction of 
the optical path along the beam propagation direction is 
limited by the gap thickness through which the sample must 
flow. This could be accomplisheid by using a cell with various 
pathlengths or a cell with a wedge shaped window spacing (see 
Figure 9b) to provide a range of optical pathlengths. Smaller 
source beams would pass through the thinner portions of the . 
cell, reducing the intersection of the incident beam and 
particle dispersant volume to avoid coincidence counts. The 
other alternative is to restrict the field of view of the 
scattering collection optics so as to only detect scatterers 
in a very small sample volume, which reduces the probability 
of multiple particles in the measurement volume. So 
particularly in the case of very small particles, a focused 
laser beam intersected with the limited field of view of 
collection optics must be used to insure single particle 
counting. However, this system would require correction of 
larger beams for coincidence counts based upon counts m 
smaller beams. To avoid these count errors, this disclosure 
proposes the use of a small interrogation volume for small 
particles, using multiple scattering angles, and a 2 
dimensional detector array for counting large numbers of 
particles above approximately 1 micron at high speeds. 

Three problems associated with measuring very small particles 
are scattering signal dynamic range, particle composition 
dependence, and Mie resonances. The low angle scattered 
iS?ensity per particle changes by almost 6 orders of magnitude 
between 0.1 and 1 micron particle diameter. Below 
aSp^SxJmately 0.4 micron, photon multiplier tubes PMT) are 
llldld to meLure the minute scattered light signals. Also the 
scattered intensity can change by a factor of 10 between 
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particles of refractive index 1.5 to 1.7. However, the shape 

of the scattering function (as opposed to the amplitude) vs. 
scattering angle is a clear indicator of particle size, with 
very little refractive index sensitivity. This invention 
proposes measurement of multiple scattering angles to 
determine the size of each individual particle, with low 
sensitivity to particle composition and scattering intensity. 
Since multiple angle detection is difficult to accomplish with 
bulky PMT's, this invention also proposes the use of silicon 
photodiodes and heterodyne detection to measure low scattered 
signals .from particles below 1 micron. However, the use of any 
type of detector and coherent or non-coherent detection are 
claimed . 

Spherical particles with low absorption will produce a 

transmitted light component which interferes with light 
diffracted .by the particle. This interference causes 
oscillations in the scattering intensity as a function of 
particle size. The best method of reducing these oscillations 
is to measure scattering from a white light or broad band 
source, such as an LED. The interference resonances at 
multiple wavelengths are out of phase with each other, washing 
out the resonance effects . But for small particles, one needs 
a high intensity source, eliminating broad band sources from 
consideration. The resonances primarily occur above 1.5 micron 
particle diameter, where the scattering crossection is 
sufficient for the lower intensity of broadband sources. So 
the overall concept may use laser sources and multiple 
scattering angles for particles below approximately 1 micron, 
and broad band sources with low angle scattering or total 
scattering for particle size trom approximately 1 micron up to 
thousands of microns. We will start with the small particle 
detection system.. 

Figure 1 shows a configuration for measuring and counting 
smaller particles. A light source is projected into a sample 
cell, which consists of two optical windows for confining the 
flowing particle dispersion. The light source in Figure 1 
could also be replaced by an apertured light source as shown 
in Figure lA, This aperture, which is in an image plane of the 
light source, blocks unwanted stray light which surrounds the 
source spot. In the case of laser sources, this aperture may 
be used to select a section of uniform intensity from the 
center of the laser crossectional intensity profile. In all 
figures in this disclosure, either source configuration is 
assumed. The choice is determined by source properties and 
intensity uniformity requirements in the sample cell.^ So^ 
either the light source, or the apertured image of the light 
source, is collimated by lens 1 and a portion of this 
collimated beam is split off by a beam splitter to provide the 
local oscillator for heterodyne detection. While collimation 
between lenses 1 and 2 is not required, it provides for easy 
transport to the heterodyne detectors 3 and 4. Lens 2 focuses 
the beam into a two-window cell as a scattering light source 
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for particles passing through the cell. The scattered light 
is collected by two optical systems;, a high angle heterodyne 
system for particles below approximately 0 . 5 microns and a low 
angle non-coherent detector for 0.4 to 1.2 micron diameter 
particles . Each system has multiple detectors to measure 
scattering at multiple angles. Figure 1 shows a representative 
system, where the representative approximate mean scattering 
angles for detectors 1, 2, 3,. and 4 are 10, 20, ao, and 80 
degrees, respectively. However, other angles and numbers of 
detectors could be used. All four scattering intensity 
measurements are used for each particle passing through the 
intersection of the field of view of each of the two systems 
with the focused source beam. Detectors 1 and 2 use non- 
coherent detection because the signal levels for the larger 
particles measured by these two detectors are sufficiently 
large to avoid the complexity of a heterodyne system- Also the 
Doppler frequency for particles passing through the cell at 
meter per second speeds are too low to accumulate many cycles 
within the single particle pulse envelope at these low 
scattering angles. The Doppler frequencies may be much larger 
at larger scattering angles where the heterodyne detection is 
needed to measure the small scattering intensities from 
smaller particles . 

Lens 4 collects scattered light from particles in the flowing 
dispersant. Slit 1 is imaged by lens 4 into the cell, The 
intersection of the fays passing through that image and the 
incident source beam define the interrogation volume 1 where 
the particle must reside to be detected by detectors 1 and 2. 
Detectors 1 and 2 each intercept a different angular range of 
scattered light. Likewise for lens 3, slit 2 and detectors 3 
and 4. The intersection of the rays by back-projection image 
of slit 2 and the source beam define interrogation volume 2 
for the heterodyne system. The positions of slit 1 and slit 2 
are adjusted so that their interrogation volumes coincide on 
the source beam. In order to define the smallest interaction 
volume', the images of the two slits should coincide with the 
minimum beam waist in the sample cell. These slits could also 
be replaced by other apertures such as pinholes. A portion the 
source beam, which was split off by a beamsplitter (the source 
beamsplitter), is reflected by a mirror to be expanded by a 
negative lens 5. This expanded beam is focused by lens 6 to 
match the wavefront of the scattered beam through lens 3. This 
matching beam is folded through slit 2 by a second 
beamsplitter (the detector beamsplitter) to mix with the 
scattered light on detectors 3 and 4. The total of the optical 
pathlengths from the source beamsplitter to the particle in 
the sample cell and from the particle to detectors 3 and 4, 
must match the total optical pathlength of the local 
oscillator beam from the source beamsplitter through the 
mirror, lenses 5 and 6, the detector beamsplitter, and slit 2 
to detectors 3 and 4. The difference between these two total 
optical pathlengths must be less than the coherence length of 
the source to insure high- interferometric visibility in the 
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heterodyne signal. The scattered light is Doppler shifted by 
the flow velocity of the particles in the cell. By mixing this 
Doppler frequency shifted scattered light with unshifted light 
from the source on a quadratic detector (square of the 
combined E fields) , a Doppler beat frequency is generated in 
the currents of detectors 3 and 4. The current oscillation 
amplitude is proportional to the square-root of the product of 
the source intensity and the scattered intensity- Hence, by 
increasing the amount of source light in the mixing, the 
detection will reach the Shot noise limit, allowing detection 
of particles below 0.1 micron diameter. By using a sawtooth 
drive function to vibrate the mirror with a vibrational 
component perpendicular to the mirror's surface, introducing 
optical phase modulation, the frequency of the heterodyne 
carrier can be increased to produce more signal oscillations 
per particle pulse. During each rise of the sawtooth function 
and corresponding motion of the mirror, the optical frequency 
of the light reflected from the mirror is shifted, providing a 
heterodyne beat signal on detectors 3 and 4 equal to that 
frequency shift. Then the mirror vibration signal could be 
used with a phase sensitive detection, at the frequency and 
phase of the beat frequency, to improve signal to noise. This 
could also be accomplished with other types of optical phase 
modulators (electro-optic and acousto-optic) or frequency 
shifters (acousto-optic) . The reference signal for the phase 
sensitive detection could be provided by a separate detector 
which measures the mixture of light, which is reflected by the 
moving mirror (or frequency shifted by another device) , with 
the unshifted light from the source. 

For particles above approximately 0.4 microns, signals from 
all 4 detectors will have sufficient signal to noise to 
provide accurate particle size determination. The theoretical 
values for these 4 detectors vs. particle size may be placed 
in a lookup table. The 4 detector values from a measured 
unknown particle are compared against this table to find the 
two closest 4 detector signal groups, based upon the least 
squares minimization of the function: 

(S1-S1T)^2 + (S2-S2T) ,'^2 + (S3-S3T,) ^^2 + (S4-S4T) ^2 

where S1,S2,, S3,S4 are signals from the 4 detectors, 
SIT, S2T, S3T, S4T are the theoretical values of the four signals 
for a particular particle size, and ^2 is the power of 2 or 
square of the quantity preceding the . 

The true size is then determined by interpolation between 
these two best data sets based upon interpolation in 4 
dimensional space. The size could also be determined by using 
search algorithms, which would find the particle size which 
minimizes the least square error while searching over the 4 
dimensional space of the 4 detector signals. For particles of 
size below some empirically determined size (possibly around 
0.4 micron), detector 1 and 2 signals could be rejected for 
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insufficient signal to noise,- and only the ratio of the 
signals (or other function of two signals) from shot noise 
limited heterodyne detectors 3 and 4 would be used to size 
each particle. If the low angle signals from detectors 1 and 2 
are needed for small particles^ they could be heterodyned with 
the source light using the same optical design as used for 
detectors 3 and 4. In any case^ only signals with sufficient 
signal to noise should be used in the size determination, 
which may include only the use of detectors 1 and 2 when 
detector 3 and 4 signals are low. The look up table could also 
be replaced by an equation in all 4 detector signals, which 
would take the form of: particle size equals a function of the 
4 detector signals. These techniques, least squares or 
function, could be extended to more than 4 detectors. For 
example, 3 detectors could be used for each system, discarding 
the low angle non-coherent detection when the signal to noise 
reaches unacceptable levels. In this case^ a 6-dimensional 
space could be searched, interpolated, or parameterized as 
described above for the 4 detector system. This disclosure 
claims the use of any number of detectors to. determine the 
particle size, with the angles and parameterization functions 
chosen to minimize size sensitivity to particle composition. 

By tracing rays back from slits 1 and 2, the field of view for 
systems 1 and 2 are determined, as shown in Figure 2. The 
traced rays and source beam converge into the interrogation 
volTome, where they all intersect. Figtire 2 shows these rays 
and beam in the vicinity of this intersection volume, without 
detailed description of the converging nature of the beams. 
The intersection volume is the intersection of the source beam 
and the field of view of the detector. In this case, the /beam 
from slit 1 may be wider than that from slit 2, so that the 
source beam and slit 2 field of view fall well within the 
field of view of slit 1. And the source beam falls well within 
the field of view of slit 2 . By accepting only particle signal 
pulses which show coincideiice with pulses from detector 4 
(which has the smallest intersection with the source beam, 
shown by the crosshatched area) , the interrogation volume is 
matched for all 4 detectors. The source beam could also have a 
rectangular crossection, with major axis aligned with the long 
axis of the slits. This would reduce the edge effects for 
particles passing near to the edges of the beam. The slit 
images are designed to be much longer than the major axis of ^ 
the source beam, so that both slits only need to be aligned in 
the direction perpendicular to the source major axis. This 
provides for very easy alignment to assure that the 
intersection of images from detectors 3 and 4 and the source 
beam fall within the intersection of images from detectors 1 
and 2. The slit position could also be adjusted along the 
optical axis of the detection system to bring the crossover 
point of both detector fields of view to be coincident with 
the source beam. Another configuration is shown in Figure 2b, 
where slit 2 is wider than slit 1. Here detector 2 defines the 
smallest common volume as indicated by the cross hatched area. 
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And so only particles which are counted by detector 2 can be 
counted by the other detectors. All other particles detected 
by the other detectors, but not detected by detector 2, are 
rejected because they do not produce concurrent signals in 
every detector- This process can be extended to more than 4 
detectors. In some cases three or more detectors per optical 
system may be required to obtain accurate size measurement. In 
this case, the size could be determined by use of a look up 
table or search algorithm. 

The data for each particle would be compared to a group of 
theoretical data sets. Using some selection routine, such as 
total RMS difference, the two nearest size successive 
theoretical sets which bracket either side of the measured 
set would be chosen. Then the measured set would be used to 
interpolate the particle size between the two chosen 
theoretical sets to determine the size. The size determination 
is made very quickly (unlike an iterative algorithm) so as to 
keep up with the large number of data sets produced by 
thousands of particles passing through the sample cell. In 
this way each particle could be individually sized and counted 
according to its size to produce a number-vs . -size 
distribution which can be converted to any other distribution 
form. These theoretical data sets could be generated for 
various particle refractive indices and particle shapes. 

In general, a set of design rules may be created for the 
intersection of fields of view from multiple scattering 
detectors at various angles. Let us define a coordinate system 
for the incident light beam with the z axis along the 
direction of propagation and the x axis and y axis are both 
perpendicular to the z axis, with the x axis in the scattering 
plane and the y axis perpendicular to 'the scattering plane. 
The scattering plane is the plane which includes the source 
beam axis and the axis of the scattered light ray. In all 
cases the detector slits are oriented parallel to the y 
direction. Many configurations are possible, including three" 
different configurations which are listed below; 

1) The incident beam is smaller than the high scattering angle , 
detector field crossection, which is smaller than the low 
scattering angle detector field crossection. Only particle 
pulses that are coincident with the high angle detector 
pulses are accepted. The incident beam may be spatially 
filtered in the y direction, with the filter aperture imaged 
into the interaction volume. This aperture will cut off the 
Gaussian wings of the intensity profile in the y direction, 
providing a more abrupt drop in intensity. Then fewer small 
particles, which pass through the tail of the intensity 
distribution, will be lost in the detection noise and both 
large and small particles will see the same effective 
interaction volume . 
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2) The incident beam is larger than the low scattering angle 
detector field crossection, which is larger than the high 
scattering angle dete'ctor field crossection. Only particle 
pulses that are coincident with the high angle detector 
pulses are accepted. The correlation coefficient of the 
pulses or the delay (determined by cross correlation) 
between pulses is used to insure that only pulses from 
particles seen by every detector are counted. 

3) The incident beam width and all fields from individual 
detectors progress from small to large size. Then particles 
counted by the smallest entity will be sensed by all of the 

, rest of the detectors . Only particles sensed by the 
smallest entity will be counted. 

The two detector pairs, 1+2 and 3+4, could also be used 
independently to measure count vs. size distributions. The 
lower angle pair could only measure down to the size where the 
ratio of their angles is no longer sensitive to size and the 
scattering crossections are too small to maintain signal to " 
noise. Likewise for the high angle detectors, they can only 
measure up to sizes where their ratio is no longer monotonic 
with particle size. However, absolute scattered signal levels 
could be used to determine the particle size outside of this 
size region. Since extremes of these operational ranges 
overlap on the size scale, the two pairs could be aligned and 
operated independently. The small angle detectors would miss 
some small particles and the high angle detectors would miss 
some large particles. But the two independently acquired 
particle size distributions could be combined using their 
particle size distributions in the size region where they 
overlap. Scale one distribution to match the other in the 
overlap region and then use the distribution below the. overlap 
from the high angle detectors for below the overlap region and 
the distribution from the low angle detectors for the 
distribution above the overlap region. In the overlap region, 
thie distribution starts with the high angle result and blends 
towards the low angle result as you increase particle size. 
Detector triplets could also be used, where the largest angle 
of the low angle set and the lowest angle of the high angle 
set overlap so as to scale the scattering measurements to each 
other. 

in some cases, the angular range of each of the heterodyne 
detecbtors must be limited to maintain heterodyne signal 
visibility. 

The flat window surfaces could be replaced by spherical 
surfaces with centers of curvature which coincide with the 
center of the interrogation volume. Then the focal positions 
of all of the beams would remain in the same location for 
dispersing liquids with various refractive indices. These 
systems can also be designed using fiber optics, by replacing 
beamsplitters with fiber optic couplers. Then the vibrating 
mirror could be replaced by a fiber optic phase modulator. 
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Figure 3 shows an alternate optical configuration for Figure 
Ir where the low angle scattering system is placed on the 
opposite side of the cell from the high angle system. In some 
cases ^ this configuration will facilitate the mechanical 
design of the support structure for the cell and optical 
systems . 

The detector currents from the low angle system and the high 
angle system must be processed differently. Every particle 
passing through the interaction volume will produce a pulse in 
the detector current. Detectors 1 and 2 will show simple 
pulses, but detectors 3 and 4 will produce modulated pulses. 
The heterodyne detection measures the Doppler beat frequency 
as the particle passes through the beam. So each heterodyne 
pulse will consist of a train of oscillations which are ^ 
amplitude modulated by an envelope determined by the intensity 
profile of the incident beam, as shown in Figure 101 for a 
Gaussian beam profile. The heterodyne signal must pass through 
a high pass filter or bandpass filter (to remove the large 
local oscillator offset) and then an envelope detector (see 
Figure 4) to remove the heterodyne oscillations, producing the 
signal envelope for further processing- This preprocessing 
envelope detection is used in the process steps below. 

For small particles the heterodyne signals will be buried in ■ 
laser source noise. Figure 5 shows an additional detector 5 
which measures the intensity of the local oscillator laser 
noise. If we define a heterodyne detector current as II and 
the detector 5 laser monitor detector current as 12 we obtain 
the following equations which hold for each of the heterodyne 
detectors . 

II = sqrt( R*Io(t) * Is(t) ) * COS(F*t + A) + R*Io 
II = sqrt{ R*Io(t) *S(1-R) lo ) * COS(F*t + A) + R*Io 



12 = K*Io(t) 
where : 

COS (x) = cosine of x 

K is a constant which includes the product of the 
reflectivities of the beamsplitter 1 and beamsplitter 3 

R and (1~R) are the effective reflectivity and transmission of 
the beam splitters, respectively 
R = R2*R3* (1~R1) 
(1-R) == (1-R2) * (1-R3) 

R2 is the reflectivity of beamsplitter 2 
R3 is the reflectivity of beamsplitter 3 
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sqrt(x) = square root of x 

Io(t) is the source beam intensity as function of time t 

F is the heterodyne beat frequency at a heterodyne detector 
due to the motion of the scatterer in the sample cell. And A 
is an arbitrary phase angle for the particular particle - 
Is(t) is the scattered light intensity from the particle: 

Is(t) = S* (l-R) *Io (t) where S is the scattering efficiency 
or scattering crossection for the particle 

The light source intensity will consist of a constant portion 
loc and noise n(t): 

Io(t) = loc + n{t) 

We may then rewrite equations for II and 12: 

11 - sqrt (S* (1-R) *R) * (loc +n (t) ) *COS (F^t -H A) + R*(Ioc + 

n(t)) . 

12 = (loc + n (t) ) 

The heterodyne beat from a particle traveling with nearly 
constant velocity down the sample cell will cover a very 
narrow spectral range with high frequency F. For example, at 1 
meter per second flow rate, the beat frequency would be in the 
megahertz range. If we use narrow band filters to only accept 
the narrow range of beat frequencies we obtain the narrow band 
components for II and 12: 

Ilnb = sqrt,(,S*(l-R)*R)*Ioc*COS(F*t + A) + R*n(t) 
I2nb = K*n(t) 

where we have assumed that n(t) is much smaller than loc. And 
also n(t} is the portion of the laser noise that is within the 
electronic narrowband filter bandwidth (see below) . 

The laser noise can be removed to produce the pure heterodyne 
signal, Idiff, through the following relationship: 

Idiff = Ilnb - (R/K)*I2nb = Sqrt (R* ( 1-R) *S) *Ioc*COS (F*t + 
A) 

This relationship is realized by narrowband filtering of each 
of the II and 12 detector currents. One or both of these 
filtered signals are amplified by programmable amplifiers, 
whose gains and phase shifts are adjustable. The difference of 
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the two outputs' of these amplifiers is generated by a 
difference circuit or differential amplifier. With no 
particles in the beam, the gain and phase shift of at least 
one of the programmable amplifiers is adjusted, under computer 
or manual control, to minimize the output of the difference 
circuit (i.e. (gain for I2)*R/K = 1, assuming gain for 11=1) 
At this gain, the source intensity noise component in the 
detector 3 or detector 4 beat signal, with particles present, 
is eliminated in the difference signal, which is fed to an 
analog to digital converter (A/D) , through a third narrowband 
filter, for analysis to sense the beat signal buried in noise. 
This filtered difference signal could also be detected by a 
phase locked loop, which would lock in on the beat frequency 
of current from the heterodyne detector. 

The particle dispersion flow rate could also be adjusted to 
maximize the heterodyne signal, through the electronic 
narrowband filter, from flowing particle scattered light. 

This entire correction could also be accomplished in the 
computer by using a separate A/D for each filtered signal and 
generating the difference signal by digital computation inside 
the computer. The phase and gain adjustments mentioned above, 
without particles in the beam, could be adjusted digitally. 
Also these gain adjustments could also be determined from 
measurement of the signal offsets Ildc and I2dc (the average 
value of the signal due to the local oscillator) . If the 
scattering component of the heterodyne signal is negligible 
compared to the offset caused by the local oscillator, this 
adjustment could be determined from measurements taken with 
particles in the beam. In this case, the contribution from the 
source intensity noise should be proportional to the offset 
level because the noise is the same percentage of the average 
level of the intensity in both II and 12. Then the coefficient 
ratio R/K in the equation for Idiff can be calculated from: 

R/K = Ildc/I2dc 

Where Ildc and I2dc are the average of the unfiltered signals 
II and 12, respectively. And the gain (or digital multiplier) 
of 12 is then I2dc/Ildc (relative to a gain for signal II = 
1) . 

If both signals were digitized separately, other correlation 
techniques could be used to reduce the effects of source 
intensity noise. Beamsplitter 2 and 3 reflections are 
adjusted to obtain shot noise limited heterodyne detection, 
with excess laser noise removed by the difference circuit. 

The noise correction techniques described on the prior pages 
(and Figure 5) can be applied to any heterodyning system by 
simply adjusting the filtering of currents II and 12 to pass 
the signal of interest, while blocking the low frequency 
component (loc) of Io(t)- Excess laser noise and any other 
correlated noise component, which is present in both the 
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heterodyne signal and the light source, can be removed from 
the signal of interest through this procedure. One application 
is dynamic light scattering, where the heterodyne signal is 
contaminated by laser source noise in the optical mixing 
process. The filters on II and 12 would be designed to pass 
the important portion of the Doppler broadened spectrum (using 
a lower frequency broad band filter or high pass filter 
instead of the high frequency narrow band filter) and to 
remove the large signal offset due to the local oscillator . 
Then by using the subtraction equation described below 
(where the narrow band filter is replaced by said broad band 
filter in all equations) the effects of laser noise can be 
removed from the Doppler spectrum, improving the particle size 
accuracy. 

Idiff = Ilnb - (R/K)*I2nb = Sqrt (R* ( 1-R) *S) *Ioc*COS (F*t + 
A) 

In this case, the heterodyne signal is the sum of many COS 
functions with various frequencies and phases. The noise, 
common to both the heterodyne signal and incident light source 
intensity, will still be completely removed in Idiff. In the 
case of fiber optic heterodyning systems, the laser monitor 
current, 12, could be obtained at the exit of the unused 
output port of the fiber optic coupler which is used to 
transport the light to and from the particle sample, because 
this port carries light only from the optical source, without 
any scattered light. 12 can be measured with a light detector 
at any point in the optical system where the light source 
intensity vs. time is available. This subtraction shown in the 
equation above could be accomplished by the analog difference 
circuit or by digital subtraction after digitization of both 
the filtered contaminated signal and the filtered source 
monitor as outlined previously. This procedure could also be 
accomplished using the unfiltered signals, but with much 
poorer accuracy due to the large signal offsets . 



Figure 6 shows the system with some additional features . The 
sample cell windows contain spherical surfaces with center of 
curvature at the interaction volume. The light source beam and 
detector acceptance cones pass through these spherical 
surfaces in order to avoid focal shift of the source and 
detector beams when the refractive index of the dispersing 
fluid is changed. The heterodyne detector currents from 
detectors 3 and 4 are passed through a high pass filter to 
remove the large local oscillator current and then (after 
completing the noise removal described above) they are passed 
through an envelope detector to remove the heterodyne 
oscillation due to the Doppler shifted spectrum of the 
scattered light from the moving particles. As mentioned 
earlier, this Doppler frequency may be increased by vibrating 
the mirror so as to add phase modulation to the local 
oscillator. This will provide more signal oscillations per 
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signal pulse. After the high pass or narrowband filter, the 
signal will consist of a sinusoid which is amplitude modulated 
by the scattering pulse due to the particle's transit through 
the source beam (see Figure 101) . The envelope of this 
modulated sinusoid is measured by an envelope detector as 
shown in Figure 4. The resulting single pulse is digitized by 
an analog to digital converter (A/D) before analysis by a 
computer. This process is similar for each of detectors 3 and 
4 . Since lower angle scattering produces lower Doppler 
frequency, the lower scattering angle signals are usually 
measured without heterodyne detection when the signals are 
large. So for large signal levels. Detectors 1 and 2 do not 
require heterodyne detection; but a heterodyne optical system, 
as used for detectors 3 and 4, could be used for detectors 1 
and 2 if the signal levels were small. Then the vibrating 
mirror phase modulator, shown in figure below, could be used 
to increase the heterodyning frequency. If the signals are 
large, the scattered light current pulses from detectors 1 and 
2 can be digitized directly before computer analysis, without 
envelope detection. The analysis of these signals is described 
below. 

One other aspect of this invention is a means for auto- 
alignment of the optics. As shown in figure 2, the source beam 
and all four fields of view, from the four detectors, must 
intersect at the same point to all see scattering from the 
same particle. Images of slits 1 and 2 define the point where 
the view fields from each detector pair (1+2 or 3+4) 
intersect- The slit apertures usually only need alignment in 
one direction, perpendicular to the slit, but position 
adjustment may also be needed along the optical axis of the 
detector system to place the intersection between the fields 
of view from detectors 1 and 2 (or 3 and 4) on the source 
beam. Pinhole or rectangular apertures must be aligned in two 
orthogonal directions which are in the plane perpendicular to 
the optical axis of the scattering detection optical system. 
Either one or both slits may be adjusted to obtain alignment. 
Figure 5 shows an example where both slit positions are 
optimized by a computer controlled micro-positioner. For 
example, the digitized signals from detector 2 and detector 3 
could be digitally multiplied (after the envelope detector) 
and the resulting product integrated or low pass filtered to 
produce a correlation between the two detector signals. The 
position slit 1 is adjusted until this correlation signal is 
maximum with particles flowing through the interaction volume. 
If needed, both slits 1 and 2 may be moved to optimize this 
correlation signal. In general these should be small 
adjustments because the spherical window surfaces will prevent 
large beam refractions and focal shifts due to changing 
particle dispersant refractive index. In systems with large 
beam shifts, the slits may need to be moved perpendicular and 
parallel to the optical axis of each optical system to 
maximize the correlation between the detectors . This could be 
accomplished with dual axis micro-positioners, which could 
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also be used when tfie slits are replaced by pinholes or 

rectangular apertures, which require alignment in two 
orthogonal axes perpendicular to the optical axis of the 
scattering detection optical systern . 

Figure 2 shows larger fields of view for detectors 1 and 2 
than for detectors 3 and 4 . This is accomplished with slit 1 
wider than slit 2 or by larger magnification for lens 4 than 
for lens 3 (image of the slit in the 'interaction volume) . 
Hence the alignment of slit 1 is much less critical then slit 
2 because image of slit 1 in the interaction volume is wider 
and has larger depth of focus than for slit 2 . By placing 
computer controlled micro-positioners on slit 1 and slit 2, 
the system can be aligned by using the correlation between the 
signals. The micro-positioners move each slit perpendicular to 
the long axis of the slit opening and perpendicular to the 
optical axis of that lens. The alignment procedure is 
described below: 

1) With low concentration of particles in the flow stream, 
adjust the position of slit 2 to maximize the correlation 
(using an analog multiplier and RMS circuit) between the 
signals from detectors 3 and 4. At this point the 
intersection of the fields of views of both detectors cross 
at the incident beam and the signals are maximum. 

2) Then adjust the position of slit 1 until the correlation of 
detectors 1 and 2 with detectors 3 and, 4 is a maximum. 
After this adjustment both detectors 1 and 2 view the 
intersection defined by step 1. 

During particle counting and measurement, only particles seen 
by both detectors 3 and 4 are counted by all of the detectors, 
because they are a subset of the particles seen by detectors 1 
and 2 . By using different slit image sizes and using the 
smaller slit images to determine count acceptance, the system 
will accept only particles which are seen by all four 
detectors. If the slit images from detectors 3 and 4 are 
larger than the images from detectors 1 and 2, then detectors 
1 and 2 would be adjusted before adjusting 3 and 4; and 
detectors 1 and 2 would select which particles are counted. 
The general rule is that the detector images which have the 
smallest intersection with the incident beam are adjusted 
first and they determine which particles will be counted. The 
slit widths are chosen to create one slit image with a small 
intersection volume and the other with a larger intersection 
volxame so that when a particle is detected in the smaller 
volume, it is clearly within the larger volume. The smaller 
slit image only needs to cross the incident beam near to its 
image plane. Then the larger slit image only needs to cover 
the intersecting volume to insure that it sees all of the 
particles passing through the smaller slit image. Then by only 
counting particles detected by the smaller slit image, only 
particles which are seen by both detectors will be counted. If 
the slit images were comparable in size, very precise 
alignment of both slit images with each other would be 
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required and the correlation between the detector signals 
would be needed to choose which particles to count. This 
comparable sized slit ca3e is also claimed in this disclosure. 
Also the replacement of slits with rectangular apertures or 
pinholes is also claimed, but with the requirement for two 
axis alignment as indicated previously. 

Figure 6b shows an example of scattered detection pulses from 
the four detectors. These signals are measured after a high 
pass filter for each of detectors 1 and 2, and after a high 
pass filter and envelope detector for each of detectors 3 and 
4, The high pass filters could be replaced by narrow band 
filters. This data describes the case where the particle 
passes through a corner of the volume which is common to the 
source beam and the field of view from detector 4 (see figure 

2) . Detectors 1/ 2, and 3 show similar profiles as a function 
of time as the particle passes through the interaction volume. 
However the signal from detector four is truncated at the 
leading edge due to the edge of the detector field of view. 
Over the region where the particle is well within the detector 
fields of view, each detector signal will maintain the same 
ratio with another detector signal as the detector signal 
amplitudes follow the particle passing through the source 
crossection intensity distribution. This region of stable 
signal ratio must be determined in order to eliminate the 
effects of the variation in source intensity by ratioing pairs 
of detector signals. Each of the four detector signals is 
digitized and the ratio of signal from the detector with the 
minimum interaction volume with one of the other detectors is 
calculated at each A/D (analog to digital conversion) sampling 
point. The A/D may be only turned on by a comparator during 
the period where all the detector signals are above a noise 
threshold, between times Tl and T2 in Figure 6b. In this case 
the ratio between detectors 3 and 4 is used to determine the 
optimum portion the sampled data to use. The ratio of detector 
4 /detector 3 increases as the particle enters the field of 
view of detector 4. Once the particle is completely inside the 
field of view, the ratio between the two signals is nearly 
constant even though the individual signals are changing due 
to the source intensity distribution non-uniformity. 
Eventually, the signal levels drop and the signal ratio 
becomes very noisy. If we assume that there are 2 0 samples 
between Tl and T2, we could measure the variance of the ratio 
for samples 1' through 5 and then the variance of the ratio for 
samples 2 through 6, and so on up to samples 16 through 20. 
The 5 sample set with the lowest variance for the detector 3/ 
detector 4 ratio would be chosen to determine the detector to 
detector ratios for all detector combinations for that 
particle, either by choosing the sample in the middle (sample 

3) of that set or by averaging all 5 ratios to obtain an 
averaged ratio for the 5 samples in the set. The assumption of 
20 samples and 5 samples per set is an example. This invention 
claims any appropriate data set size and segmentation. 
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The pulses shown in Figure 6b are the result of some prior 
electronic filtering and envelope detection. The signals from 
detectors 1 and 2 will be simple pulses which may be cleaned 
up by a high pass filter before the A/D conversion. The 
signals from the heterodyne detectors 3 and 4 are the product 
of a pulse and a sinusoid. The pu'lse may consist of a 
megahertz sine wave, amplitude modulated by the intensity 
profile of the source beam over a period of about 100 
microseconds, depending upon the size of the interaction 
volume and the particle flow velocity. This oscillatory signal 
sits on top of a large offset due to the local oscillator 
intensity. This offset and other source noise components may 
be removed from the heterodyne signal by high pass or narrow- 
band electronic filtering. The power spectrum of these pulses 
will reside in a 100 kilohertz band which is centered at 1 
megahertz . Hence a narrow-band filter may provide optimal 
signal to noise for the heterodyne signals. After the 
filtering, the signals could be digitized directly for digital 
envelope detection or an analog envelope detector could be 
used to remove the 1 megahertz carrier, reducing the required 
sampling rate to only 10 to 20 samples per pulse instead of 
400 samples per pulse. By using a dual phase lock-in amplifier 
with reference oscillator set to the heterodyne frequency (1 
megahertz in this example), extremely high signal to noise 
could be obtained by measuring the filtered signal without the 
envelope detector. By using the zero degree and quadrature 
outputs of the dual phase lock-in amplifier, the phase 
sensitive signal would be recovered even though the reference 
and signal carriers are not in necessarily phase. 

The particle counting rate can also be increased by digitizing 
the peak scattered signal (directly from detectors 1 and 2 and 
after the envelope detector from detectors 3 and 4) from each 
particle instead of digitizing many points across the 
scattering pulse and finding the peak digitally. This is 
accomplished by using an analog peak detector whose output is 
digitized in sync with the positive portion of the signal 
pulse derivative and reset by the negative portion of the 
derivative. Then only one digitization is needed for each 
particle, as shown in figure 7. The negative comparator 
switches on when the input signal drops below the reference 
setting and the positive comparator switches on when the 
signal is greater than the reference setting. 

Another variation of this concept triggers on the actual 
signal instead of the derivative, as shown in Figure 8. When 
the signal rises above a preset threshold, the positive 
comparator takes the peak detector out of reset mode. As the 
signal rises, the output of the peak detector (see Figure 9) 
follows the input signal until the signal reaches a peak. 
After this point, the peak detector holds the peak value with 
a time constant given by the RC of the peak detector circuit. 
The input signal drops below this as it falls down the 
backside of the peak. When the signal reaches some percentage 
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of the peak value, the A/D is triggered to read and then reset 
the positive comparator. This percentage value is provided by 
a voltage divider (shown in the figure as 0 . 5X voltage 
divider, but other divider ratios would also be appropriate 
between approximately 0.2 to 0.8) which determines the 
reference level for the negative comparator. The A/D is only 
triggered once per signal pulse and measures the peak value of 
the pulse. Using this circuit, the detector with the smallest 
interaction volume generates the A/D trigger for all of the 
other detectors, so that only particles seen by all of the 
detectors are counted. 

The detector signals are either digitized directly, peak 
detected with the circuit in Figure 7 or Figure 8, or 
integrated and sampled at a lower rate. The signal can be 
continuously integrated (up to the saturation limit of the 
integrator) . Then the integrated signal needs to be sampled 
only at points of zero slope in the integrated signal, between 
each pulse. By subtracting the integrated values on either 
side of the pulse, the integral of each pulse is sampled 
separately without having to sample the pulse at a high 
sampling rate. Also each pulse could be sampled at a lower 
rate and then a function could be fit to these samples to 
determine the peak value of the pulse. This should work 
particularly well using Gaussian functions which model the 
intensity profile of the laser beam. The parameters of the 
best fit Gaussian solution directly provides the peak, half 
width, or integral of the pulse. In any case, the final signal 
from each pulse will be analyzed and counted. One problem 
associated with conventional particle counters is the incident 
beam intensity profile in the interaction volume. Identical 
particles passing through different portions of the beam will 
see different incident intensity and scatter light 
proportionally to that intensity. But the scattered intensity 
also depends upon the particle diameter, dropping as the sixth 
power of the diameter below 0.3 microns. So the effective 
interaction volume will depend upon particle diameter and 
detection noise, because particles will not be detected below 
this noise level. Therefore small particles will be lost in 
the noise when they pass through the tail of the intensity 
distribution. This means that larger particles have a larger 
effective interaction volume than smaller particles and 
therefore the number distribution is skewed in favor of large 
particles. This invention includes a method for creating a 
particle diameter-independent interaction volume, using signal 
analysis. The systems in Figures 1, 2, and 3 use at least 2 
detectors per scatter collection system to remove the incident 
intensity dependence by using the ratio of two scattering 
angles to determine the particle size. Also for any number of 
detectors, ratios between any pair of detectors could be used 
to determine the size of particles in the size range covered 
by that pair. Instead of detector pairs, detector triplets or 
quadruplets, etc. could also be used with appropriate 
equations or lookup tables to determine the size of each 
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particle independent of the incident intensity on the 
particle. In the case of detector pairs, both scattered 
signals, S(A1) and S{A2), are proportional to the scattering 
function, at that angle, times the incident light intensity: 

S (al) = K^IO^F{D,al) 

S(a2) = K^I0^F(D,a2) 

where 10 is the incident intensity, K is an instrumental 
constant in this case, and F{D,A) is the scattering per 
particle per unit incident light intensity for a particle of 
diameter D, at scattering angle a. The variable ^'a'' can refer 
to a single angle or a range of angles over which the signal 
is collected- Then al and a2 would refer to the al range of 
angles and a2 range of angles, respectively. For more than two 
detectors, there is a similar equation for each detector 
signal for angles al, a2, 'a3, a4, etc. The scattering signals 
S (a) may be the pulse peak value or pulse integral of the 
envelope of the heterodyne signal (detectors 3 and 4) or of 
the direct non-coherent signals (detectors 1 and 2) . So The 
ratio of the scattering at two angles is equal to 
F(D,al) /F(D,a2) , which is independent of incident intensity 
and relatively independent of position in the Gaussian beam 
profile of a laser. Figure 10 shows a conceptual plot of 
number of particles vs. S(al)/S(a2) and S(a2). A plot of 
number of particles vs. S(al)/S(a2) and S(al) could also be 
used. The scattering signal, for any diameter D, will show a 
,very narrow range of S(al)/S(a2) but a broad range of S(a2). 
Particles passing, through, the peak of the laser intensity 
profile will produce pulse peak amplitudes at the upper limit, 
maximum S(a2). The surface, describing this count 
distribution, is determined by fitting a surface function to, 
or by interpolation of, this count surface in figure 10, This 
surface function provides the parameters to determine an 
accurate particle count, because S(al)/S(a2) is a strong 
function of particle size, but a very weak function of 
particle path through the beam. By setting an acceptance 
threshold for S(a2) at a certain percentage of this maximum 
value, separately at each value of S(al)/S(a2), only particles 
passing through a certain voliame (independent of particle 
size) of the beam will be accepted and counted. Because the 
particles, which are counted by these detectors, are all much-, 
smaller than the source beam crossection, they all have the 
same probability functions for describing the percentage of 
particles passing through each segment of the beam. Therefore, 
at any value of S(al)/S(a2), the shape of the count vs. S(a2) 
function is nearly identical when you normalize the function 
to maximum S(a2). By setting a count-above threshold at a 
certain percentage of maximum S(a2) (but well above the noise 
level) at. each value of S(al)/S(a2), only particles passing 
through a certain portion of the interaction volume, with 
acceptable signal to noise, will be counted and sized, as 
shown in figure 10b. The noise threshold is chosen so that all 
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particles with signals above that level will be accurately 
sized based upon the scattering signals. 

This analysis is usually done for the detector with the 
smallest interaction volume, the heterodyne system in the case 
of Figurel. All four detectors are used to determine the 
particle diameter^^ but the acceptance criteria is determined 
by only detector 3 and 4 (S (al) =detector 3 and S (a2 ) =detector 
4) . This analysis can also be performed, individually, on any 
pairs of signals, as long as the noise threshold is always the 
same percentage of the maximum value of signal used for the 
horizontal axis in Figure 10b. These counts are accumulated 
into a set of size ranges, each range defines a different size 
channel. In many cases, each size range has a very narrow 
width in size and S{al)/S{a2). The optimum channel size width 
is the minimum width which still contains sufficient particle 
counts in that channel to avoid statistical errors. Hence, the 
distribution of S(a2) for the range of S(al)/S(a2) within a 
certain channel can determine the S(a2) acceptance limit for 
counts in that channel. This is accomplished by only counting 
particles with S(a2) above a certain percentage of the maximum 
S(a2) for that channel. If the theoretical scattering 
efficiency changes substantially across any channel, the S(a) 
for each count is divided by the theoretical scattering 
efficiency indicated by the size corresponding to the 
S(al)/S(a2) for that particle. This may be especially 
important below 0.4 microns where the scattering efficiency 
drops as the inverse of the sixth power of the particle 
diameter. 

The shape and width of the S(a2) profile is determined by how 
sharply the source crossectional intensity distribution drops 
off at the edges of the beam. If the beam profile was a step 
function, the effective interaction volume would be only 
weakly particle size dependent near the edges. This shape can 
be accomplished by spatially filtering the source, with the 
spatial filter aperture in a plane conjugate to the 
interaction volume. Then an image of the aperture, which is 
smeared by aberrations and diffraction limits, defines the ' 
sharpness of intensity drop at the edges of the beam. The 
intensity tails of the Gaussian beam are cut off by the 
aperture, which could be sized to cut off at any appropriate 
percentage of the peak intensity to limit the variation of 
scattering from a particle as it traverses the beam. The beam 
crossectional intensity distribution may also be shaped by use 
of appropriate apodization of the beam or by using diffractive 
beam shapers . 

Figure 10b describes this concept for data collected in one 
size channel. The pulse signal is collected and stored for 
each pulse above the count noise threshold. But only pulses 
which are within the range of S(al)/S(a2) for that channel are 
collected into that channel. The frequency distribution of 
counts at each pulse level is plotted for a beam with a 
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Gaussian intensity profile. The problem is that some particle 
pulses fall below the noise threshold and are not counted. The 
amount of missed particles depends upon the scattering 
efficiency of the particles. For smaller particles with lower 
scattering efficiency;, a higher percentage of particles will 
be lost below the noise threshold- So the count error will be 
particle size dependent. The source beam can be spatially 
filtered to cut off the low intensity wings of the source 
intensity distribution. Then the count distribution would be 
as shown in the ""'with aperture (ideal)" curve and no particles 
would be lost in the noise. This could be accomplished by 
using a rectangular spatial filter that cut the wings off in 
the Y direction, because the particle flows in the XZ plane 
and the tails in this plane are actually measured in each 
pulse shape. However, the image of the spatial filter aperture 
in the interaction volume will be aberrated and diffraction 
limited as shown by the '^^with aperture (aberrated)'' curve. In 
this case a few particles may still be lost in the noise and a 
count threshold must be, set above this level to reject all 
questionable particle pulses. The maximum S(a2) value changes 
in each channel due to the change in scattering efficiency for 
the particles in that channel. As long as the count threshold 
is set to be the same percentage of that maximum S(a2) for 
each channel, all channels will lose the same percentage of 
particles and the distribution will be correct. Without this 
channel specific threshold, the smaller particle channels will 
lose a larger percentage of particles than the larger particle 
channels and the distribution will be skewed towards larger 
particles. This assumes that the sample is a homogeneous 
mixture of all the particle sizes and that the sufficient 
count exists- in each channel to obtain an accurate estimate of 
the maximum S(a2) . 

The method described above handles the variations caused by 
the particle passing through various random paths in the 
interaction volume. This method can also correct for the 
variations due to random positions along the path where the 
digitization occurs . Therefore the peak detectors or 
integrators could be eliminated. The signal from the envelope 
detector (detectors 3+4) and direct signals from detectors 1+2 
could be digitized directly at approximately 3 points per 
pulse. The maximum signal data point from each pulse would be 
added to the data list for input to the analysis described 
above. Any deviation from the peak value would not be a 
problem because the ratio of the signals determines the size 
and all four detectors will be low by the same percentage if 
they are not sampled at the peak intensity position in the 
interaction volume . 

Also all the signals could be digitized directly after the 
high pass (or narrow band) filter on the detectors (detectors 
3 and 4 are high pass filtered to remove local oscillator 
current and detectors 1 and 2 could also be high pass filtered 
to remove low frequency noise) . Then all of the analog and 
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digital operations (phase sensitive detection, envelope 
detection, etc.) could be done digitally but at the cost of 
reducing data collection rates. Also the source could be 
modulated for detectors 1 and 2 to use phase sensitive 
detection (lock-in amplifier) when their signals are low. 

All of the optical design and algorithm techniques described 
in this disclosure may leave some residual size response 
broadening which may be particle size dependent. This 
instrument response broadening is determined by measuring a 
nearly mono-sized particle sample (such as polystyrene 
spheres) . For example, due to noise or position dependence in 
the beam, a certain size particle will produce a range of 
S(al)/S(a2) values as it repeatedly passes through different 
portions of the interaction volume. In any event, the 
broadening may be removed by solving the set of equations 
which describe the broadening phenomena. If the broadening is 
relatively the same for all size particles, the response 
broadening can be described by a convolution of the broadened 
niomber distribution response and the actual number vs. size 
distribution. Iterative deconvolution algorithms may be used 
to deconvolve the measured number vs. measured parameter 
distribution to obtain size resolution enhancement. This 
resolution enhancement will work for any ergodic stochastic 
process, where the broadening statistics are stable over time. 
This idea could be applied to (and is claimed for) any 
broadened counting phenomena with stable stochastic or 
deterministic broadening mechanisms. In particle counting 
measurements the amount of scatter from a particle may vary 
due to the random orientation and position of a particle as it 
passes through the exciting light beam, or by other structural 
and optical noise sources. The counting and classification of 
each of a group of identical particles will not produce a 
narrow peak when plotting count number vs. measured parameter. 
Here ^^measured parameter'' refers to the parameter which is 
measured from each particle to determine its size. Examples of 
measured parameters are scattering optical flux amplitude, 
ratio of flux from two scattering angles, a function of fluxes 
from multiple scattering angles, or the decrease in intensity 
due to particle scattering and absorption, as will be 
described later in this disclosure. The peak of the number- 
vs. -measured parameter function from a group of monosized 
particles will be broadened in a predictable way. This 
broadening can be determined experimentally with a calibrated 
group of particles (by measuring the response from monosized 
.particle samples) or it can be calculated theoretically based 
upon models for the random and deterministic broadening 
sources. Then the entire system can be modeled using a matrix 
equation, where each column in the matrix is the broadened 
measured parameter distribution from a certain sized particle. 
This broadening is reproducible as long as a large number of 
particles are counted for each trial. The matrix equation is 
described by the following relationship: 
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Where Nm is the vector of values of the measured (broadened) 
number-vs . -measured parameter distribution and N is the vector 
of values of the actual particle huinber-vs . -size distribution 
which would have been measured if the broadening mechanisms 
were not present. is a matrix vector multiply. The number 

distribution is the number of particles counted with parameter 
amplitudes within certain ranges. It is a differential 
distribution which describes counts in different channels or 
bins, each bin with a different range of parameter, which may 
be size, scattering ratio, etc.. M is a matrix of column 
vectors with values of the broadened number-vs . -measured 
parameter function for each particle size in N. For example, 
the nth column of M is a vector of values of the entire 
measured number- . vs . -measured parameter distribution obtained 
from a large ensemble of particles of the size which is 
represented by the nth element of vector N. This matrix 
equation can be solved for the particle number-vs . -size 
distribution, N, by matrix inversion of M or by iterative 
inversion of the matrix equation. This particle number-vs. - 
size distribution can be determined by using this matrix 
equation in many different forms. The term ^^measured 
parameter'' in this paragraph can refer to many size dependent 
parameters including: scattering signal amplitude (pulse peak 
or integral, etc. ) ^ the ratio (or other appropriate 
mathematical relationship) between scattered signals at two or 
more different angles, or even particle diameter (a broadened 
particle size distribution determined directly from a 
broadened process can also be ^^unbroadened'' by using broadened 
particle size distributions for each monosized sample column 
in matrix M) . So we solve for N, given Nm and M. 

If each column of M is simply a shifted version of the prior 
coliomn, then the instrument response is shift invariant and 
the relationship is a convolution of N with the system impulse 
response IMP: 

Nm = IMP **N 

where is the convolution operator 

For this case, deconvolution algorithms may be used tp solve 
for N, given Nm and IMP. 

The generalized matrix equation above may also include the 
effects of coincidence counting. As discussed earlier, over 
one million particles should be counted for a uniform volume 
distribution to be accurately determined in the large particle 
region. In order to insure low coincidence counts, the source 
spot size in the interaction region must be reduced to 
approximately 2 0 microns in width so that the particle 
concentration can be raised to count 1 million particles at 
flow rates of 1 meter per second in a reasonable time. For 
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example, the worst case is slit 1 being the largest slit, 
because then the largest interaction volume might be 
approximately 2 0 micron x 20 micron x 200 micron, for example 
If we require approximately 5 volumes per particle to avoid 
coincidence counts then the inter-particle spacing is 7 4 
microns. 1 million particles spaced by 74 microns (on average) 
moving at 1 meter per second will take 74 seconds to measure. 
This spot size would provide good count reproducibility for 
the worst case of uniform volume distribution. However, a 20 
micron spot and the corresponding detector fields of view may 
be difficult to align, requiring larger source spot size with 
a higher coincidence level. Even with a 20 micron spot, some 
coincidences will be seen at the 74 micron particle spacing. 
These coincidences can be corrected for by including their 
effects in the generalized matrix equation. If M were 
correcting for coincidences, a column in matrix M which 
corresponds to the large size end of vector .N will have 
negative values in the region corresponding to the small size 
end of vector Nm, because the larger particles will block the 
scattered light from smaller particles which are ahead or 
behind that larger particle in the source beam. Also a column 
in matrix M which corresponds to the small size end of vector 
N will have a tail of positive values in the region 
corresponding to the large size end of vector Nm, because some 
smaller particles will be counted coincidently with the larger 
particles and increase their measured size relative to their 
actual size. The effects of coincident counts can be mitigated 
by using a wedge shaped cell as shown in Figure 9b. The cell 
consists of two windows at an angle so as to produce regions 
of different optical path along the cell. This cell could 
replace the cells in Figure 11 and 12. The red rays define the 
edges of the source beam. Then at any point along the wedge 
direction, only particles smaller than a certain size may pass 
through that portion of the cell. The size distributions 
gathered at different points along the wedge, from the 2 
dimensional detector array in Figures 11 and 12, may be 
combined by correcting the count in the larger particle areas 
for coincidently counted smaller particles by using counts in ' 
the smaller particle regions of the wedge. This correction can 
be accomplished by solving a matrix equation of the form 
shown previously. 

The correction for coincidences may also be accomplished by an 
iterative procedure, which solves for N, given Nm, and then 
corrects each scattered signal for coincidences. Each 
scattered signal, SI and S2, consists of light scattered (or 
light lost due to absorption or scattering) from all the 
particles in the interaction volume. Ideally, the particle 
concentration is low and most of the time each scattering 
event is from a single particle. But for the general case, 
multiple coincident particles can be modeled by the following 
equation : 



Ai= SUMj ( G{Ni,Nj) A j ) 
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where SUMj means summation over the j index. Ai is the 
'^particle signal" for a particle of the ith size bin in the 
particle size distribution. Particle signals can include SI, 
S2, or the log of attenuation or obscuration (described later 
in this disclosure) due scattering and absorption of a 
particle. G{Ni,Nj) is a function which describes the most 
probable total particle signal from a combination of particles 
of ith and jth sizes based upon their particle numbers (or 
concentrations), Ni (for the ith size) and Nj (for the jth 
size) . Since the combination of particles in the interaction 
volume is a random process, G(Ni,Nj) represents the sum of all 
combinations (given Ni and N j ) , weighted by their probability 
functions . 

In the case of signals S(al) and S(a2), the procedure for 
determining the number vs. size distribution is the following: 

1) Use the surface plot of Figure 10 and 10b to determine the 
raw number distribution Nm. 

2) Solve the matrix equation Nm = M*N for the true number 
distribution N, 

3) Recalculate S(al) and S(a2) using the equation above and 
the distribution N: 

S(al)i= SUMj ( Gl(Ni,Nj) S(al)j) 
S(a2)i- SUMj( G2(Ni,Nj) S(a2)j)' 

4) Do steps 1 through 3 again 

5) repeat iteration loop of step 4 until the change in number 
distribution N between successive loops is below some 
threshold. 



For particles between approximately 1 and 10 microns, the 
ratio of scattered intensities at two angles below 
approximately 3 degrees scattering angle is optimal to provide 
highest size sensitivity and accuracy. A white light source or 
broad band LED should be used to reduce the Mie resonances for 
spherical particles. Above 10 microns, the measurement of 
total scatter from a white light source provides the best size 
sensitivity and depth of focus for a spatially filtered 
imaging system as shown in Figure 11. A white light or broad 
band LED source is spatially filtered by lens 1 and pinhole 1 
to provide a well collimated beam through lens 2. If a well 
collimated source beam is required to measure scattering at 
very low scattering angles (for large particles), a laser 
source might also be used. This collimated beam passes through 
a cell consisting of two windows, which confine the flowing 
particle dispersion. Lens 3 focuses this collimated beam 
through pinhole 2, which removes most of the scattered light 
from the beam. This transmitted beam is transferred to a 2 
dimensional detector array through lens 4, which images the 
center of the sample cell onto the array. This array will see 
dark images of each particle on a bright background due to the 
light lost through scattering or absorption by the particle. A 
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iDeamsplitter after lens 3 diverts a! portion of thd light to an 
aperture and lens 5. The aperture defines a narrow pencil of 
light through the cell and a small scattering volume, lowering 
the probability of coincidence counts tor detectors 1 and 2, 
which are near to the focal plane of lens '5. The aperture is 
optimally placed in the optical plane which is conjugate to 
the center of the sample cell, through lens 3. Conjugate 
planes are a pair of image and object planes of an optical 
system. Detectors 1 and 2 are nominally placed in the optical 
plane which is conjugate to the optical source, through lenses 
1, 2, 3 and 5. Detectors 1 and 2 measure scattered light at 
two angles, which are nominally below 3 degrees for larger 
particles but which can cover any angular range appropriate 
for the size range of the particle detector. Also more than 
two detectors could be used to increase the size range for 
this portion of the particle detection system. These detectors 
could also be annular ring detectors centered on the optical 
axis to reduce sensitivity to particle shape. For example, 
detector 1 could measure a scattering angular region around 1 
degree; and detector 2 could measure around 3 degrees. By 
combining the particle scattering pulse signals from these 
detectors, by ratio or polynomial, a relatively monotonic 
function of particle size is created without strong Mie 
resonances (due to the white light source and signal ratio) . 
Detectors 1 and 2 count and size particles in much the same 
way as the system in Figure 1 . The concept is to use two 
angles to remove the variations in scattering intensity due 
to particles passing through different portions of the 
incident beam and to reduce calculated size sensitivity to 
particle and dispersant composition. Particles of size between 
approximately 1 and 10 microns could be handled by detectors 1 
and 2 of Figure 11; and particles above approximately 5 
microns are handled by the 2-dimensional detector array. The 
two size distributions from these two measurements are 
combined with blending in the overlap region between 5 and 10 
microns . 

The 2-dimensional detector array is imaged into the center of 
the sample cell with magnification corresponding to 
approximately 10x10 micron per array pixel in the sample cell 
plane. A 10 micron particle will produce a single dark pixel 
if it is centered on one pixel or otherwise partially darkened 
adjacent pixels. By summing the total light lost in these 
adjacent pixels, the total light absorbed or scattered outside 
of pinhole 2 for each single particle in the view of the array 
is determined. The particle concentration is limited to 
prevent coincidence counting in each separate 10x10 micron 
projection through the sample cell. At low concentrations, any 
group of contiguous pixels with reduced light levels will 
represent a single particle. And the total percentage light 
lost by these contiguous pixels determines the particle size. 
All pixels below a certain percentage of their non-obscured 
values are accepted as particle pixels. All contiguous 
particle pixels are then combined as representing one 
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particle. This is accomplished by siommihg the pixel values of 
contiguous pixels and comparing that sum to the sum of those 
same pixels without the presence of a particle. This works 
well for smaller particles where the total scattered light is 
well outside of the pinhole 2 aperture, because then the total 
percentage drop represents the scattering and absorption 
extinction of the particle. For larger particles, a larger 
portion of the scattered light will pass through pinhole 2 and 
cause a deviation which will not agree with the total 
theoretical scattering extinction- This can either be 
corrected for in the theoretical model by calculating the 
actual percentage loss for larger particles by integrating the 
actual scattered light outside of the pinhole, or the particle 
can be sized directly by counting contiguous pixels, because 
for the larger particles the pixel size may be less than 0.1 % 
of the total crossectional area of the particle. The accuracy 
of this calculation is improved by adding partial pixels at 
the edge of the particle image based upon their attenuation as 
a fraction of the attenuation of nearby interior pixels. Hence 
if a pixel in the interior of the image is attenuated by 10% 
and an edge pixel is attenuated by 4%, that edge pixel should 
count as 40% of its actual area when added to the sum of all 
contiguous attenuated pixels to determine the total 
crossectional area and size for that particle. Otherwise the 
theoretical loss per particle could be used. 

This detector array system has an enormous particle size 
dynamic range. The particle will remove approximately the 
light captured by twice its crossectional area. So a 2 micron 
particle will reduce the total light flux on a 10x10 micron 
pixel by 8 percent. But the entire array of 1000x1000 pixels 
can cover a crossection of 10x10 millimeters. So the size 
range can cover 2 microns to 10000 microns. The size dynamic 
range is almost 4 orders of magnitude. The smallest particles 
are detected by their total light scattering and absorption. 
For very large particles, the angular extent of the scattering 
pattern may fall within the aperture of pinhole 2 . Then the 
summed light from all the contiguous pixels may not indicate 
accurate size. For the larger particles, the actual imaged 
size is determined by counting contiguous pixels. Pixels at 
the outer boundary are counted as partial pixels based upon 
the amount of light lost as a fraction of the amount lost from 
pixels in the interior of the contiguous set. The light loss 
in each pixel is determined by storing the light value for 
each pixel without particles in the sample cell and 
subtracting the particle present values these stored values. 
The source intensity can also be monitored to normalize each 
pixel measurement for light source intensity fluctuations . 

In order to avoid smeared images, the detector array must 
integrate the current from each pixel over a short time to 
reduce the distance traveled by the flow during the exposure. 
This may also be accomplished by pulsing the light source to 
reduce the exposure time. Smearing in the image can be 
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corrected for using deconvolution techniques. But the 
scattering extinction measurements will be accurate as long as 
the contiguous pixel groups do not smear into each other. 
Simply add up all of contiguous pixel signals (from the 
smeared image of the particle) after presence, of the particle 
to determine the particle scattering attenuation and size. If 
the particle image is smaller than one pixel, then the 
attenuation of that pixel is the scattering extinction for 
that particle. Essentially, you are measuring nearly the total 
amount of light scattered or absorbed by the particle during 
the exposure. Using this total lost optical flux divided by 
the incident intensity provides the scattering crossection for 
the particle, even if the particle is not resolved by the 
optical system or if that loss is distributed over more pixels 
than expected from perfect imaging of the particle. This is 
the power of this technique, A 10mm by 10mm detector array, 
with 10x10 micron pixels, can measure particle diameters from 
a few microns up tolO mm, with thousands of particles in the 
source beam at one time. The 10 mm particles will be sized 
directly by adding up pixels and multiplying the interior 
pixels by 1 and the edge pixels by their fractional 
attenuation and adding all of the pixels up to get the total 
crossectional, area and size. A 5 micron particle, centered on 
one 10x10 micron pixel, will attenuate that pixel by 50% (the 
total scattering extinction crossection is approximately twice 
the actual particle area, outside of the Mie resonance 
region) . In both cases the particles are easily measured. You 
are simply adding up all of the signal differences (signal 
without particle - signal with particle) of contiguous changed 
pixels to get the total light lost due the particle. Pinhole 
2 blocks all of the light scattered outside of the angular 
range of the pinhole 2, whose maximum scattering angle is 
equal to the inverse tangent of the pinhole 2 radius divided 
by the focal length of lens 3. So the signal difference 
(signal without particle - signal with particle) is the amount 
of light scattered by the particle at scattering angles above 
this maximum angle of the pinhole, including any light 
absorbed by the particle. The particle size is determined 
using scattering theory and the ratio of signal change (signal 
without particle - signal with particle) to the signal without 
a particle. 

Image smearing could also be reduced by using pulsed flow. The 
particle sample flow would stop during the period when the 
light source is pulsed or when the detector array is 
integrating. Then a flow pulse would push the next slug of 
sample into the detector array field of view before the next 
signal collection period. The sample would be approximately 
stationary during the signal collection on the detector array. 
This pulsing could be accomplished by pressurizing the 
particle dispersion chamber and using a pulsed valve to leave 
short segments of the sample dispersion through the source 
beam interaction volume . 
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The nearly parallel window cell could also fee replaced by a 
wedge shaped cell which would control tHe particle count in 
different size regions, as discussfed above (see Figure 9b) . 

Non-spherical particles present anSfeh'er problem for single 
particle sizing: non-syitrnietridaii scattering patterns. Assume 
that the incident light beam is propagating along the Z 
direction and the XY plane is perpendicular to the Z 
direction, with origin at the particle. The XZ plane is the 
center scattering plane of the group of scattering planes 
which are intercepted by detectors 1,2,3, and 4. Each detector 
subtends a certain range of scattering angles, both parallel 
and perpendicular to the center scattering plane. For 
spherical particles, the scattering pattern is symmetrical 
about the Z axis and the scattering function could be 
described in cylindrical coordinates as a function of Z and of 
radius R from the Z axis, at some distance ZO from the 
scattering particle. However, for npn-spherical particles the 
scattering pattern is not symmetrical about the Z axis at ZO. 
The 2 -dimensional array in. Figure 11 measures approximately 
the total light lost to scattering or absorption at all 
scattering angles, in all scattering planes. Hence it will 
produce particle size estimates which are related to the total 
crossectional area of the particle, for both spherical and 
non-spherical particles, without sensitivity to particle 
orientation. But detectors 1, 2, 3, and 4 in Figures 1, 3, or 
11 measure the scattering only over scattering planes close to 
the XZ plane (or a limited range of scattering planes) . If the 
pattern is not symmetrical, the particle size estimate will 
depend upon the orientation of the particle. So a group of 
particles with identical crossectional areas, but random 
orientations, would be reported over a wide range of particle 
crossectional area and size. This particle size distribution 
width could be corrected by deconvolution of the number vs. 
size distribution, . as described by the matrix equation shown 
previously, where matrix M would describe the broadening of a 
count distribution from a group of particles, each with the 
same particle volume, but with all possible orientations . But 
the theoretical model would change with the particle shape. 
Another way to reduce spread is to use two sets of detector 
systems, one centered on the scattering plane which is +45 
degrees with respect to the XZ plane and the other at -45 
degrees from that plane, to sample two perpendicular particle 
orientations and maintain the optimum orientation for 
heterodyne detection. The average of the size distributions 
from these two systems would reduce the spread of the 
distribution. Another more effective method is to collect all 
of the scattered planes at a certain scattering angle, using 
the system shown in Figure 12. A light source is focused into 
the sample dispersion. This focused spot is imaged onto a 
pinhole which removes unwanted background light. The light 
passed by the pinhole contains the incident light beam and the 
scattered light from the particles . This light is collected by 
lens 3 which projects the light onto two masks, using a 
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beamsplitter. Each mask contains ah annular aperture which 
defines the range of scattering angle accepted by the 
collection optics. Lens 4 collects high angle scattered light 
passing through mask 1 and focuses it onto detector 1. 
Likewise the low angle scatter is measured by mask 2^ lens 5 
and detector 2 . 

Figure 12 shows the annular aperture for mask 1, defining 
equal scattering collection in all the scattering planes. The 
ratio of signals from detector 1 and 2 would precisely 
determine the average radius of a non-spherical particle, 
without size broadening of the system response due to random 
particle orientation. The beam splitter and dual mask concept 
could also be applied to the system in Figure 11. Lens 5 and 
detectors 1 and 2 (all of Figure 11) would be replaced by lens 
3, the beam splitter, the dual mask system, and detectors 1 
and 2 (all of Figure 12), with the masks in the same optical 
plane as detectors 1 and 2. Masks 1 and 2 act as angular 
filters which only pass scattered light in a certain range of 
scattering angles. The 2-dimensional array in Figure 11 is 
already insensitive to particle orientation and needs no 
modification . 

The particle concentration must be optimized to provide 
largest count levels while still insuring single particle 
counting. The concentration may be optimized by computer 
control of particle injection into the flow loop which 
contains the sample cell, as shown in Figure 13. Concentrated 
sample is introduced into flow loop 2 through sample vessel 2, 
The sample vessel may also contain a stirring means for 
maintaining a homogenous dispersion in the vessel. Pump 2 
pumps the dispersion around the loop to provide a homogenous 
dispersion in the loop and to prevent loss of larger particles 
through settling. A second flow system, flow loop 1, is 
attached to flow loop 2 through a computer controlled valve 
with minimal dead space . The computer opens the valve for a 
predetermined period to inject a small volume of concentrated 
dispersion into loop 1. The optical system counts the 
particles and determines the probability of coincidence 
counting based upon Poisson statistics of the counting 
process . The computer then calculates the amount of additional 
particles needed to optimize the concentration and meters out 
another injection of concentrated sample into loop 1, through 
the valve. Actually, both the concentration and pump speed for 
loop 1 may be controlled by computer to optimize counting 
statistics. When the particle concentration is low, higher 
pump speed will maintain a sufficient particle count rate for 
good count statistics. The optimum concentration may be 
different for different detectors and detection systems. 
Therefore the computer valve may adjust the concentration to 
various levels in succession. At each concentration level data 
is taken with the appropriate detector (s) or detector array 
for a sufficient period and flow rate to accumulate enough 
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counts to reduce the count uncertainty (due to Pdis3ori 
statistics ) to an acceptable level. 

Another consideration for Figure 6b is the determination of 
signal baseline. The baseline for the scattered signals must 
be determined for each detector. Digitized values; measured 
before and after the scattered signal pulse, determine the 
signal baseline to be subtracted from the pulse signal, by 
interpolation of those values through the pulse region. These 
regions before and after each detector pulse should be chosen 
to be before and after the widest pulse of the group (in some 
rare cases, the pulse with the largest amplitude should be 
used If the signals are lost in noise) . Then the baseline will 
certainly be determined from values in a region where no 
particle scattering has occurred in each of the detectors. 

The system shown in Figure 11 can also be modified to look at 
only scattered light over a certain angular region, instead of 
the total light removed from the beam by absorption and/or 
scattering. Figure 14 shows such an optical system where the 
light source is spatially filtered by lens 1 and pinhole 1. 
Lens 2 collimates and projects the source beam through the 
particle sample, which is imaged onto the 2 dimensional- 
detector array by lens 3. An annular spatial mask (or spatial 
filter) IS placed in the back focal plane of lens 3 to only 
pass scattered light over a certain range of scatter angle as 
defined by the inner and outer radii of the annular spatial 
mask, which is similar to mask 1 shown in figure 12. The very 
low angle scattering and incident beam are blocked by central 
stop of the annular aperture in the back focal plane of lens 
3. Hence the detector array 1 sees an image of the particles, 
and the sum of the contiguous pixels associated with each 
particle's image is equal to the scattered light from that 
particle over the angular range defined by the aperture (or 
spatial mask) in the back focal plane of lens 3. A beam 
splitter splits off a portion of the light to a second annular 
filter (in the back focal plane of lens 3) and detector array 
2. The angular ranges of the two annular filters are chosen to 
produce scattered values which are combined by an algorithm 
which determines the size of each particle. The sum of signals 
from contiguous pixels which view the same particle are 
analyzed to produce the particle size. One such algorithm 
would be a simple ratio of the corresponding sums (the sum of 
contiguous pixels from the image of each particle) from the 
same particle detected by both arrays. The key advantage is 
that when the particle size becomes too small to size 
accurately by dimensional measurements on the image 
(resolution is limited by pixel size) then the total scattered 
light from each particle may be used to determine the size. 
And if the total scattered light is sensitive to particle 
composition, then the ratio of the two scattering signals can 
be used to determine the particle size more accurately. In 
figure 14, scattered light is only present when a particle is 
present. In figure 11 the particle image creates a decrease in 



30 



wo 2005/091970 



PCT/US2005/007308 



light, from a bright background level > oh the 2-Dimensional 
array xn the corresponding pixels, i^hiie in the system of 
figure 14 the particle image creates aih increase from a dark 
background level. If the particle is smaller than a single 
pixel, then the amount of scattered light measured by that 
pixel will indicate the total light scattered from that 
particle in the angular range defined by the focal plane 
aperture, providing that particle's size. If more than one 
pixel is associated with a particle, those pixel values are 
summed together to obtain the scattered signal from that 
particle in a similar fashion as described before for figure 
11. The only difference is that the increase in pixel signals, 
relative to the signal without particles, are summed to 
produce the total light scattered from that particle in the 
angular range of the annular aperture in figure 14. In Figure 
11, the decrease in pixel signals, relative to the signal 
without particles, are summed to produce the light lost due 
to scattering outside of pinhole 2 or, absorption by that 
particle. Signal to background should be better for figure 14, 
but with higher sensitivity to particle composition and 
position in the sample cell. The depth, of focus and signal to 
noise should be better for figure 11 than for figure 14, 
because the pixel values drop by the total light scattered and 
absorbed by the particle in figure 14 as opposed to the light 
increasing by. only the amount scattered through a narrow range 
of scattering angles defined- by the aperture. As with all 
other systems described in- this disclosure, these ideas can be 
extended to more than two detector arrays or more than two 
scattering angles, sim:ply by adding more annular spatial masks 
and detectors by using beamsplitters. Also the 2~dimensional 
array optics in Figure 11 could be combined with those in 
Figure 14 to provide total scatter information for determining 
size and to provide the unscattered intensity for each pixel 
to normalize the pixel scatter data of Figure 14 detectors for 
the incident light intensity which may vary across the beam. 
In this way, each pixel in the detector array creates a small 
independent interaction volume with low coincidence 
probability . But yet contiguous pixels can be combined to 
measure particles of sizes approaching the dimensions of the 
entire detector array's image in the sample cell. The size 
dynamic range is enormous. Figure 14 co.uld also be used with a 
source beam which is focused into the sample cell to reduce 
the interaction volume and increase the beam intensity and 
scattered signal. In this case the center portion of the 
annular mask must be increased in size to block the diverging 
light from the source so that each detector array only sees 
scattered light. 

The optical source used with the detector arrays in Figures 11 
and 14 could be a pulsed broad band source such as a xenon 
flash lamp which produces broadband light to wash out the Mie 
resonances, and a short light pulse to freeze the motion of 
particles flowing through the cell. 
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One problem with the techniques described above is coincidence 
counting- The cell path must be large to pass the largest 
anticipated particle (eicc^pt for thfe wedge cell shown in 
figure 9b, where the pathlength chainges across the cell) . 
Hence for these collimated systems, many small particles may 
be in any sample volume seen by a single pixfel- These 
coincidences could possibly be eliminated by measuring at 
various particle concentrations;, but in order to count 
sufficient large particles to obtain reasonable count 
accuracy, the concentration must be raised to a level where 
more than one small particle is present in the sample cell 
volume, seen by each pixel. The scattered signals from these 
multiple particles can be separated, to be counted 
individually, by measuring their settling velocities. This is 
accomplished by the optical system shown in Figure 15. A light 
source is collimated and spatially; filtered by lenses 1 and 2 
and pinhole 1. A modulation transfer target or mask with a 
spatially periodic transmission function is placed in the 
collimated beam to create a, sinusoidal (or other periodic 
function) intensity pattern in the collimated beam. Examples 
of the sinusoidal (or other periodic) patterns are shown in 
Figure 16. Each line in the patterns represents the peak of 
the sinusoidal transmission function which oscillates along 
the particle settling direction but is constant along the^ 
direction perpendicular to the settling. The mask consists of 
multiple regions with different spatial modulation 
frequencies. The projection of each region into the sample 
cell is imaged onto a separate detector by lens 3. The. light 
from lens 3 is split into one or more directions, each having 
a different annular spatial mask which defines a different 
range of scattering angles. Each image plane for each spatial 
mask has multiple detectors, each of which intercept light 
from only one of mask regions in the sample cell. As a 
particle settles through the sinusoidal intensity pattern, the 
scattered light on the detector is modulated because the 
scattered light is proportional to the light intensity 
incident on the particle and the mask provides a spatially 
modulated illumination field. When a particle passes through a 
region where the spatial modulation wavelength is greater than 
the particle size, the scattered light from that particle will 
attain a large modulation visibility (the ratio of peak to 
trough values will be large) . The scattered signal from the 
"largest particles will have lowest modulation visibility in 
the high spatial frequency region because the particle will 
span over multiple cycles of the spatial modulation. Larger 
particles settle faster and produce higher frequency detector 
signals, because they have a higher terminal velocity. 
Therefore, a lower spatial modulation frequency can be used 
with larger particles to increase modulation visibility while 
still maintaining high signal oscillation frequency, because 
the scattering signal frequency is equal to the product of 
settling velocity and the mask spatial frequency. The size 
range is increased by using multiple regions with different 
spatial modulation frequencies, with higher frequencies for 
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smaller slower settling particles. The area of the higher 
frequency portions of the mask are smaller to reduce the 
number of particles measured at one time by each detector, 
because typically there are much higher small particle counts 
per unit volume than for larger particles. Figure 15b shows a 
similar system with a single pinhole filter. Detector signals 
will show the same oscillation characteristics as Figure 15, 
but with a large offset due to the light transmitted by the 
pinhole. The power spectrum of the detector currents for the 
systems shown in both Figures 15 and 15b will be similar 
except that figure 15b will contain a large component at (and 
near to ) zero frequency. Figure 17 shows the power spectra 
from the low and high angle detector signals in figure 15, for 
the two detector elements, Bl and B2, which view the same 
portion, B, of the mask, but pass through different annular 
filters. Since each particle settles at a different velocity, 
each particle will produce a separate narrow peak at the same 
frequency in both of the power spectra of detectors Bl and B2 
(from signals for scattering angle 1 and scattering angle 2, 
respectively) . This is due to the fact that the detector 
signal power at any certain frequency, measured by each of 
the corresponding low and high angle detectors, will originate 
from the same particle or group of particles. Since the 
smaller particles will create a. continuum at lower 
frequencies, they can be removed from the spectrum of the 
larger particles. The corresponding single peak values from 
power spectrum of current from detector Bl at frequency fl and 
the power spectrum from current of detector B2 at frequency fl 
for example (from each scattering angle) can be ratioed (or 
analyzed by other . algorithms) to determine the size of the 
particle which created that peak in each spectra. In this way, 
multiple particles in the sampling volume can be counted 
individually. When particle size is close to the line spacing 
of the modulation target, the modulation of the scattered 
light will decrease because the signal is the convolution of 
the particle with the modulation target. However, the 
amplitude of the scattered signals at both angles will both 
decrease by the same percentage so that their ratio will still 
accurately indicate the size. Any peaks with amplitudes that 
are higher than that which would be expected from one particle 
are expected to originate from more than one particle. The 
expected amplitude for a single particle can be determined 
from the minimum value of other peaks in that frequency region 
for prior digitization sets. These multiple particle peaks can 
be either corrected for the second particle's contribution or 
eliminated from the particle count. If the particle density, 
liquid density and viscosity, are known, each individual 
particle size can also be determined by the frequency of the 
corresponding peak, by calculating the corresponding settling 
velocity and using the Stokes equation for settling to solve 
for the particle size. 

The signal frequency for each particle signal pulse could also 
be determined individually by either the timing of zero 
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crossings or by using ai phase locked loop, avoiding the power 
spectrum calculation. Each particle pulse will consist of a 
train of oscillations which are modulated by the intensity 
profile through that particular mask region. The oscillation 
amplitude and frequency provide the scattering amplitude and 
settling velocity, respectively, for that particle. The size 
can be determined from the settling velocity, if the particle 
density and fluid viscosity are known, or the size can be 
determined from the ratio of amplitudes from two different 
scattering angles (or angular ranges), or the amplitude at one 
scattering angle (or angular range) (but with possible higher 
sensitivity to particle composition) . 

The particle density or fluid viscosity can.be determined by 
combination of the scattering amplitudes and the signal 
oscillation frequency - 

Figure 18 shows another system for measuring the settling of 
particles, using crossed laser beams. An interference pattern 
is formed at the intersection of the beams. As the particles 
pass through this pattern, the scattered light is modulated, 
producing a power spectrum as described above. 

All slit apertures in this disclosure (for example, slit 1 and 
slit 2 in Figure 1) can be changed to pinholes or rectangular 
apertures, whose images at the source beam may be or may not 
be smaller than the source beam. Unlike slits, the pinholes or 
rectangular apertures may require alignment in the both of the 
mutually perpendicular X and Y directions, which are both 
approximately perpendicular to the optical axis of the 
detection system. 

Figure 19 shows more detail of the actual beam shapes in 
Figure 1 for the angular ranges specified for detectors 1 
through 4 . The scattering angle range for each detector is 
controlled by the detector size or by an aperture on front of 
a larger detector. Figure 20 shows the detail of the 
intersection of the each detector field of view with source 
beam. 

If the signal to noise is sufficient for non-coherent 
detection with any detector in Figure 1, or in any other 
variation of Figure 1 shown in this document, the local 
oscillator optics for that detector can be removed and non- 
coherent detection can be used. 

The progression of crossectional size, in the interaction 
voliome, from smallest to largest is: light source, fields of 
view from detectors 3 and 4, and the fields of view from 
detectors 1 and 2, as shown in Figure 20. The progression 
could also progress from the smallest to largest as light 
source, field of views of detector 4, detector 3, detector 2, 
and detectorl, respectively. However, this would require 
different sized apertures for each detector. This would 
require a separate lens and aperture for each detector, but 
would insure that any particle passing through the 
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intersection of the source beam and the detector 4 field of 
view, will be seen by all of the other detectors during the 
entire pulse period from detector 4. Then if all of the 
detector signal ratios are measured during a period near to 
the peak of the detector 4 signal (after the envelope 
detector) , valid scattering ratios will be recorded for all of 
the detectors . 

The particles being measured by system in Figure 19 are all 
much smaller than the crossection of the source beam. 
Therefore particles of different sizes produce the same 
count-vs . -parameter distribution for the following parameters: 

1) scattered light amplitude normalized to the maximum 
scattered light amplitude measured over all the particles of 
the same size at that scattering angle. 

2) Pulse width at certain fraction of pulse peak level 

3) delay between pulses from two different detectors 

4) correlation between pulses from two different detectors 

Any of these parameters can be used to define a threshold for 
counting particles as shown in Figures 10 and lOb^ by 
replacing the S(a2) axis with one of the parameters from 
above. If the particle count is large, the statistics of the 
above parameters will be stationary for all particle sizes. 
Then the strategy outlined for Figure 10 can be used to 
properly threshold particles of all sizes , in an equivalent 
manner, by rejecting the same percentage of particles in each 
size bin in the count-vs . -size distribution. The 3 dimensional 
surface which corresponds to the one in Figure 10, can be 
interpolated or fit to a surface function in order to 
determine the rejection threshold. Based upon the function or 
interpolated values, a rejection criteria can be determined 
which eliminates the particles with poor signal to noise and 
also removes the same percentage of particles from each size 
range so as to maintain a true particle size distribution- The 
rejection threshold is chosen to maintain a sufficiently high 
signal-to-noise for any particles which are accepted into the 
total count distribution. In fact, this process will 
computationally define an interaction volume for the source 
beam and all detector fields of view, for all particle sizes 
being detected, where all scattering signals have sufficient 
signal-to-noise, to produce accurate sizes based upon their 
amplitudes, ratios of amplitudes, or other multi-parameter 
functions of scattering amplitudes. This selection process is 
required to reduce the effects of the tails in the intensity 
distribution of the source and the spatial response tails at 
the edges of the detector fields of view, where they intersect 
the source beam. If these tails are sharpened (or cut off) by 
spatial filtering the source or by using slits, pinholes, or 
other apertures with low aberration optics for the detectors, 
the errors due to these tails are further reduced, as shown in 
Figure 10b. Also diffractive optics can be used to produce a 
^^flat top'' intensity distribution from a Gaussian laser 
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intensity profile. Then the corresponding flat top shape 
should be used to produce the functions in Figure It) and 10b 
This would improve the accuracy of the particle rejection 
threshold and the resulting particle count distribution, m 
any event, there is always some parameter which is 
statistically well described by the millions of particles 
which are detected. And by eliminating particles from the 
count based upon this parameter, you can define a group of 
particles which are sorted by the same criteria at all 
particle sizes, thereby creating an accurate size 
distribution, while removing count events which have poor 
signal to noise. 

In the sample cell with flat windows, many of the incident 
source beams and scattered light rays are at high angles of 
incidence on the sample cell windows. The interior surface of 
the window is in contact with a liquid which reduces the 
Fresnel reflection at that surface. However, the exterior 
surface is in air which can cause an enormous Fresnel 
reflection at these high incident angles. This reflection can 
be reduced by anti-reflection coating the exterior surface, 
but with high cost. A better solution is to attach prisms (see 
Figure 21) to the exterior surface with index matching optical 
adhesive. The prism .surfaces present low angles of incidence 
for the source beams and the scattered light. Even simple 
antiref lection coatings on the external prism surfaces will 
reduce the Fresnel reflections to negligible levels. A 
spherical plano-convex lens, with center of curvature near to 
the center of the cell could also be used instead of each 
, prism, with piano side attached to the window. 

Another configuration for the sample cell is a cylindrical 
tube. The particle dispersion would flow through the tube and 
the scattering plane would be nearly perpendicular to the tube 
axis and flow direction. In this case, the beam focus and 
detector fields of view would remain coincident in the 
scattering plane for various dispersant refractive indices and 
only inexpensive antiref lection coatings are needed. However, 
since the flow is perpendicular to the scattering plane, the 
heterodyne oscillations cannot be produced by the particle 
motion. The optical phase modulation mirror in the local 
oscillator arm (called ^'mirror") in Figures 1 and 19 (and 
other figures) could be oscillated to provide a heterodyne 
signal on detectors 3 and 4 as described before. This could 
also be accomplished with other types of optical phase 
modulators (electro-optic and acousto-optic) or frequency 
shifters (acousto-optic) . 

Any of the measurement techniques described can be used 
individually or in combination to cover various particle size 
ranges. Examples of possible combinations are listed below: 

For particle diameter 0.05-0.5 microns use Figure 19 with 
detectors 3 and 4 at 30 and 80 degrees, respectively. 
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Use heterodyne detection (if heeded)' . TaSke ratio of the 
detector signals. 

For particle diameter 0.4-1.2 microns use Figure 19 with 
detectors 1 and 2 at 10 and 20 degrees, respectively. 
Use heterodyne detection (if needed) . Take ratio of the 
detector signals. 

All 4 signals can also be used together for the range 0.05 to 
1.2 microns, using a 4 parameter function or lookup tables. 

Use the detector with the smallest interaction volume to 
trigger data collection. 

These angles are only representative of general ranges. For 
example, instead of 10, 20, 30, 80 degree angles, any group of 
angles with one widely spaced pair below approximately 30 
degrees and another widely spaced pair above approximately 30 
degrees would work. Each detector sees an angular range 
centered about the average angle specified above. But each 
detector angular range could be somewhat less than the angular 
spacing between members of a detector pair. In some cases, 
without optical phase modulation, the angular ranges of each 
heterodyne detector should be limited to maintain heterodyne 
inter ferometric visibility. The particle size distribution, 
below 2 microns, from this system is combined with the size 
distribution above 1 micron as determined from various other 
systems described in this disclosure, including the following 
systems : 

For particle diameter 1-10 microns use detectors 1 and 2, in 
Figures 11 or 14, at approximately 1 degree and 3 degrees ,• 
respectively. Use white light or LED source- Use ratio of 
detector signals to deteirmine individual particle size. 

For particle diameter greater than 1 microns, use the 2 
dimensional array in Figure 11 with a pulsed white light 
source (such as a pulsed xenon source) to freeze the motion of 
the flowing particles . The two dimensional array could also be 
used alone to measure greater than 1 micron, without detectors 
1 and 2. The system in Figure 14 could also measure particles 
greater than 1 micron in diameter, with one (absolute) or both 
(ratio) detector arrays. 

Another configuration is to use the scattering flux ratio of 
scattering at 4 degrees and 1 degree, in white light, for 0.5 
to 3.5 microns. And use absolute flux at 1 degree (white 
light, same system) for 3 to 15 microns- And use the 2 
dimensional array in Figures 11 or 14 for particles above 10 
microns. As noted above, each array in Figure 14 can 
separately be used to size particles by measuring the absolute 
scattered light from each particle, or light lost due to the 
particle, over one angular range. However, the system size 
response will be more composition (particle and dispersant 
refractive indices) dependent using data from a single array 
than the ratio of the corresponding measurements from both 
arrays . 
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In all of these systeiris, the white ligtit feohrce can be 
replaced by a laser. However, the particle size response will 
become more sensitive to particle and dispersant composition. 
And also the response vs. size may not ibe mono'toriic due to 
interference effects between the particle scattered light and 
light transmitted by the particle (Mie resonances), producing 
large size errors. If lasers or LEDs are required for 
collimation or cost requirements, scattering measurements can 
be made at more than one wavelength, using multiple sources, 
to reduce the composition dependence - Particle size of each 
object would be determined from all of these multi-wavelength 
measurements by using a multi-parameter function (size = 
function of multiple parameters) , by interpolation in a 
lookup table as described above, or. by a search algorithm. And 
in all cases the angles are nominal. Many different 
combinations of average angles, and ranges of angles about 
those average angles,- can be used. Each combination ha.s a 
different useable particle size range based upon size 
sensitivity, composition sensitivity, and monotonicity . All of 
these different possible combinations are claimed in this 
description. Also note that the ratio of S(:al)/S{a2) in Figure 
10 can be replaced by any parameter or function of parameters 
which have nearly exclusive sensitivity , to particle size, with 
low sensitivity to incident intensity on the particle. This 
may include functions "of scattering measurements at more than 
2 angles. 

These cases are only examples of system combinations which 
could be combined to provide a larger particle size range. 
Many other combinations are possible and claimed by this 
inventor - 

One problem associated with measuring large particles is 
settling. The system flow should be maintained at a 
sufficiently high level such that the larger particles remain 
entrained in the dispersant. This is required so that the 
scattered light measurements represent the original size 
distribution of the sample. For dense large particles, 
impracticable flow speeds may be required. This problem may be 
avoided by measuring all of the particles in one single pass, 
so that the total sample is counted even though the larger 
particles may pass through the light beam as a group (due to 
their higher settling velocities) before the smaller 
particles . 

A small open tank is placed above the sample cell region, 
connected to the cell through a tube. The tube contains a 
valve which can be shut during introduction of the particle 
sample into the tank to prevent the sample from passing 
through the cell until the appropriate time. The liquid in the 
tank is continuously stirred during the introduction of sample 
to maximize the homogeneity in the tank. A light beam may be 
passed through the mixing vessel via two windows to measure 
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scattered light or extinction to assist in determining the 
optimum amount of sample to add to obtain the largest counts 
without a high coincidence count level. The optical detectors 
are turned on and the valve is opened to allow the particle 
mixture to pass through the cell with gravitational force. 
This can also be accomplished by a v^lve below the sample cell 
or by tilting the tank up to allow the mixture to flow over a 
lip and down through the cell, as shown in Figure 22. The 
light beam could be wider than the width of the flow stream 
through the cell so that all of the particles passing through 
the cell are counted. Single particle counting is assured by 
only introducing a sufficiently small amount of sample into 
the tank. Since all of the larger particles are counted in one 
pass, the count distribution is independent of the sample 
inhomogeneity . After the entire sample has passed through the 
cell, the conventional flow system for ensemble scattering 
measurements is then turned on to circulate the sample through 
the cell. This flow rate must only be sufficient to suspend 
the particles which were too small, or too numerous, to be 
counted during the single pass. During this flow period, 
settling of the larger particles which have been counted does 
not matter because their count distribution (from the single 
pass) will be combined with the count~vs . -size distribution of 
the smaller particles (obtained during flow) to produce a 
single volume distribution over the entire size range, using 
the size distribution blending and combination methods 
described previously. 

The system in Figure 22 could also be used alone to provide 
significant cost savings by eliminating the pump and 
associated hardware. Figure 22 shows the concept for 
dispersing and measuring the particles in a single pass 
through the optical system. The particles and dispersant are 
mixed continuously in a mixing vessel which is connected to 
the optical system through a flexible tube. The mixing vessel 
is tilted while being filled, so that no sample enters the 
optical system. Then once the dispersion is well mixed, the 
mixing chamber is sealed with a gas tight cover and the mixing 
vessel is moved into an upright position to allow the 
dispersion to either fall through the sample cell under 
gravity (without gas tight cover) or to be pushed through the 
cell using gas pressure. If gas pressure is not used, the 
flexible tube could be eliminated and the mixing vessel could 
pour into a funnel on top of the sample cell. 

The systems based upon Figure 19 and similar systems, where 
scattered light at multiple -angles are measured from a single 
particle, will have a lower particle size limit due to the 
lower scattering intensity of small particles and 
insensitivity of the scattering ratio to particle size. If the 
particle refractive index and dispersant refractive index are 
known, various scattering theories can be used to calculate 
the scattering signal ratios and absolute scattering signals 
vs. particle size to provide the look up table or function to 
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calculate the size of each particle tfbm fehes^ signals. Since 
these refractive index values are hot aiways available, the 
scattering model (the effective refractive indices J may need 
to be determined empirically from scattering data measured in 
a region where the sca:tterihg iratio is independent of 
refractive index, but Where fchd absolute scattering amplitudes 
are dependent upon refractive index . in this region the 
signal ratio will determine the particle size. Rnd then the 
absolute scattering signals (relative to the incident source 
intensity) and this size value can be used to determine the 
refractive indices by using global search algorithms to search 
for the refractive indices which give the best fit to the 
particle size and absolute signal values. This process assumes 
that the entire particle ensemble is homogeneous in individual 
particle composition (however, a method is proposed for 
dealing with inhomogeneous samples ^ in Figure. 27b) . This method 
can be used for any particle which shows a ratio in the size 
sensitive region of the response, can be used to determine the 
effective refractive index of the particle by using the ratio 
and the absolute values of the scattering signals, because 
these are unique for particle and dispersant refractive 
index. The optical model for this effective refractive index 
could then be used to extend the size response range of any 
set of detectors to a size range outside of the ratio 
determined size region (where the scattering signal ratio is 
sensitive to size and the ratio cannot be used) . This process 
can extend the size to both smaller and larger particles by 
using the absolute scattering intensity in regions where the 
scattering signal ratio no longer works. Theoretically, very 
small particles are Rayleigh scatterers, where the shape of 
the angular scattering distribution is not size dependent - 
However for very small particles, the peak of the scattering 
intensity distribution scales as the 6^^ power of the particle 
diameter and the heterodyne signal scales as the 3^^ power of 
the diameter. So as the particle size decreases, the ratio of 
the intensity at two angles becomes constant, but the actual 
intensities continue to drop as the particle size decreases. 
So when the intensity ratio approaches this constant, the 
particle size algorithm should use absolute scattering 
intensity to determine the size. The absolute scattered 
intensity is proportional to a constant (which is a function 
of scattering angle, particle and dispersant refractive 
indices, and light wavelength) divided by particle diameter to 
the 6^^ power (3""^ power for the heterodyne signal) . This 
constant is determined from the absolute scattered intensities 
of particles, in the distribution, whose intensity ratios 
still provide accurate particle size. Also the functional form 
for the absolute intensities can be calculated using various 
scattering theories (Mie theory for example) . This process can 
also be used to extend the useful range to larger particles. 
As the particle becomes larger, the scattering signal ratio 
from the detector pair will become more dependent upon the 
refractive index of the particle. The absolute intensity data 
from particles, in the region where the ratio is independent 
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of particle composition, can be usfed to deterinine the 
effective composition of the particles and determine which 
theoretical scattering model to use for absolute scattering 
intensity from particles outside of that size region. Then the 
larger particles are measured using this model and the 
absolute values from each scattering detector, instead of the 
ratio of scattering signals. Absolute intensity data from 
particles, in the region where the ratio is independent of 
particle composition is preferred, but absolute scatter data 
from multiple scattering angular ranges in any size region can 
be used to determine the particle and dispersant refractive 
indices by using the global search algorithm described above. 
The scattering values depend more strongly on the ratio of 
particle refractive index to dispersant refractive index, than 
their absolute values. So in many cases, only the refractive 
index ratio must be determined by this search algorithm. 

One other remaining problem with absolute scattering 
intensity measurements is the sensitivity of pulse intensity 
to the position where the particle passes through the beam. 
Measuring the distribution of pulse amplitudes from a nearly 
mono-sized calibration particle dispersion (with a low 
coefficient of variation for the size distribution) provides 
the response of the counter for a group of particles of nearly 
identical size. This count distribution, which is the same for 
any particle whose size is much smaller than the light beam 
crossection, provides the impulse response for a deconvolution 
procedure like the one described previously. The scattering 
pulses can be selected based upon their pulse length by only 
choosing pulses with intensity normalized lengths above some 
threshold or by using the various pulse selection criteria 
listed below. This selection process will help to narrow the 
impulse response and improve the accuracy of the 
deconvolution. This process is also improved by controlling 
the intensity profile to be nearly "'flat top'' as described 
previously. 

The scattered signal from any particle is proportional to the 
intensity of the light incident on the particle. Hence as the 
particle passes through different portions of the incident 
light beam, each scattered signal will vary, but the ratio of 
any two signals (at two different scattering angles) will 
theoretically be constant as long as the field of view of each 
detector can see the particle at the same time. This can be 
insured by eliminating the signals from long intensity tails, 
of the Gaussian intensity profile of the laser beam, which 
may not be seen by all detectors. This is accomplished by 
placing an aperture, which cuts off the tails (which may be 
Gaussian) of the incident light intensity distribution, in an 
image plane which is conjugate to the interaction volume. This 
aperture will produce a tail-less illumination distribution in 
the interaction volume, providing a narrower size range 
response to mono-sized particle samples (the impulse 
response) . In the case of an elliptical Gaussian from a laser 
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diode, the aperture size could £>e chosen feo cut tHe 
distribution at approximately thd 50% points in both the x and 
y directions (which are perpendicuiat to the propagation 
direction) . Such an aperture would dause higher an^le 
diffractive lobes in the far field of the interaction volume, 
which could cause large scatter background for low scattering 
angle detectors. Since this aperture should only be used for 
measurements at high scattering angles where the background 
scatter can be avoided, the low angle detector set and high 
angle detector set may need to view separate light beams. The 
apertured beam size should be much larger than the particles 
which are being measured in that beam. Hence, to cover a large 
size range, apertured beams of various sizes could be 
implemented. The particle size distributions from these 
independent systems (different source beams or different 
detector groups) could be combined to produce one continuous 
distribution- These apertures could also have soft edges to 
apodize the beam, using known methods to flatten the beam 
intensity profile while controlling the scattering by the 
aperture. This could also be accomplished by using diffractive 
optics for producing flat top distributions from ' Gaussian beam 
profiles as mentioned earlier. Also apertures can be oriented 
to only cut into the beam in the appropriate direction such 
that the diffracted light from that beam obstruction will be 
in the plane other than the scattering plane of the detectors. 
This is accomplished by orienting the aperture edges so that 
they are not perpendicular to the scattering plane. The 
aperture edges which cut into the beam at higher levels of the 
intensity profile should be nearly parallel to the scattering 
plane to avoid high scattering background. 

The apertured beams will help to reduce the size width of the 
system response to a mono-sized particle ensemble, because the 
intensity variation of the portion of the beam which is passed 
by the aperture is reduced. Other analysis methods are also 
effective to reduce the mono-sized response width for absolute 
scattering and scattering ratio measurements. Methods which 
accept only scattering signal pulses, or portion of pulses, 
which meet certain criteria can be very effective in narrowing 
the size width of the system response to mono-sized particles. 
Some examples of these acceptance criteria are listed below. 
Any of these criteria can be used to determine which peaks or 
which portion of the peaks to be used for either using the 
scattering signal ratios or absolute values to determine the 
particle size. 

1. Choose only the time portion of both pulses where the pulse 
from the detector which sees the smaller interaction volume, 
or has the shorter duration, is above some threshold. The 
threshold could be chosen to be just above the noise level 
or at some higher level to eliminate any possibility of 
measuring one signal while the second signal is not present. 
Then either take the ratio of the signals (or ratio the 
peaks of the signals with peak detector) over that time 
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portion or the ratio of the integrals of the signals over 
that time portion. The absolute integrals or peak values 
during this time portion could ^Iso be used to determine 
size^ as described before. 
2- Only accept pulses where the separaition (or time delay 
between peaks or rising edges) between pulses from the 
multiple detectors is below some limit 

3. Only accept pulses where the width of a normalized pulse or 
width of a pulse at some threshold level is above some limit 

4. Only accept the portion of the pulses where the running 
product S1.*S2 (a vector containing the products of SI and' 
32 for every point during the pulses) of the two signals is 
above some limit. Then either take the ratio of the signals 
over that portion or the ratio of the integrals of the 
signals over that portion. 

5. Only use the portion of pulses where 

sum(Sl.*S2) / (sum(Sl) *sum(S2) ) is greater than some limit 
(sum(x) = summation of the data points in vector x) 

6. Use only the portion of the pulses where (SI . *S2) / (S1+S2) 
is greater than some limit 

7 . integrate each pulse and normalize each integral to the, 
pulse length or sample length 

8. Use only the portion of the pulses where the value of S1*S2 
is be greater than some fraction of the peak value of the 
running product Sl.*S2 

9. Integrate both signals SI and S2 only while the signal from 
the smaller interaction volume is above a threshold or while 
any of the above criteria are met . 

10. fit a function to the selected portion (based upon various 
criteria described above) of each pulse. The fitting 
function form can be measured from the signal of a particle 
passing through the center of the beam or can be based upon 
the beam intensity profile 

11. When both SI and S2 have risen above some threshold, start 
•integrating (or sample the integrators from) both signals. 
If the integrators for Si and S2 are integrating continually 
(with resets whenever they approach saturation) then these 
integrators could be sampled at various times and the 
differences would be used to determine the integrals in 
between two sample times. Otherwise the integrators could be 
started and stopped over the period of interest. These 
sampled integrals are ITIO and IT20 for SI and S2 
respectively, when each of them rises above the threshold. 
When the first signal to drop falls back down below the 
threshold, sample the integrator on each of SI and S2 
(integrals ITla and IT2a) . When the second signal (signal 
number *) to drop falls below the threshold, sample the 
integral IT^b for that signal. Use the ratio of the integral 
differences, (ITla-ITlO)/ (IT2a-IT20) , during the period 
when both signals are above the threshold to determine size. 
Accept and count only pulses where a second ratio (IT*a~ 
IT-^^^O) / (IT^b-IT*0) is above some limit. This second ratio 
indicates the fraction of the longer pulse which occurs 
during the shorter pulse. As the particle passes through the 
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light beam further away from the center of the interaction 
volume, this ratio will decrealse. Only particles which pass 
through the beam close to the center of the interaction 
volume will be chosen by only accepting pulses where the 
shorter pulse length is a large fraction of the longer pulse 
length. These pulse lengths could also be determined by 
measuring the difference in the length of time between the 
above trigger points for each pulse. Pulses with a shorter 
difference in time length are accepted into the count by 
ratioing their integrals during the period when both of them 
are above the threshold. 

These criteria can be easily implemented by digitizing SI and 
S2 and then doing the above comparisons digitally. However, 
full waveform digitization and digital analysis of 1 million 
particles may require too much time. Figures 23a and 23b show 
configurations for implementing some of these criteria using 
analog circuits . The signal digitization and computational 
load is greatly reduced by using analog equivalents to 
preprocess data before digitization. This concept is 
particularly effective when the thresholding or comparative 
functions, of the criteria described above, are replaced by 
analog equivalents; but the actual signal analysis used for 
size deteinnination is done digitally to avoid the poorer 
linearity and accuracy of the analog equivalents. An example 
of this is shown in Figure 23b, where both integrator outputs 
(one integration per signal pulse) are separately digitized by 
the A/D converters to do the amplitude or signal ratio 
calculations digitally instead of using analog ratio circuits; 
but the criteria related functions, analog multiply and 
comparator, are done analog to reduce the digitization load. 
This overall concept of using analog circuitry specifically 
for only the criteria related functions to reduce the 
digitization load is claimed by this disclosure, along with 
applications to other systems . 

All of these variations will not be perfect. Many of them rely 
upon approximations which can lead t,o variation in calculated 
size for a particle that passes through different portions of 
the beam. The important advantage is that the broadening of 
the mono-sized particle response is the same for all size 
particles which are much smaller than the source beam. 
Therefore this broadened response, which is calculated by 
measuring the count distribution from a mono-sized 
distribution or by theoretical modeling, can be used as the 
impulse response to deconvolve the count distribution of any 
size distribution . 

The intensity ratio is sensitive to size and mildly sensitive 
to particle and dispersant refractive index. Size accuracy is 
improved by using Mie theory, for the actual refractive index 
values, to calculate the scattering ratio vs. particle 
diameter function. However, sometimes these refractive indices 
are not easily determined. Three scattering angles could be 
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measured to generate a function which has reduced sensitivity 
to refractive index. 

D = A1*(S2/S1) + A2* (S2/S1) -^2 + A3*(S2/S1)^3 + 

B1*(S3/S1) + B2*(S3/S1)^2 -h B3* ( S3/S1 ) ^^3 + CI 

D = particle diameter 

Al, A2, A3, Bl, B2, B3 are constants 

51 = scattering signal at the first scattering angle (over the 
first scattering angle range) 

52 = scattering signal at the second scattering angle (over 
the second scattering angle range) 

53 = scattering signal at the third scattering single (over the 
third scattering angle range) 

Solve the set of equations: 

Di = Al*(S2/Sl)ij + A2^( (S2/Sl)ij)^2 + A3^ ( (82/81 ) i j ) ^3 4- 

Bl^(S3/Sl)ij + B2*( (S3/Sl)ij)^2 + B3* ( (S3/S1) ij ) ^3 + 

CI 

where i = diameter index 

j = index of refraction index 

and 

Slij = theoretical scattering signal over scattering angular 
range #1, for particle diameter D=Di and the jth index of 
refraction 

S2ij = theoretical scattering signal over scattering angular 
range #2, for particle diameter D=Di and the jth index of 

refraction 

( ) i j indicates that all the variables inside the 

parentheses have index i j . , 



A set of simultaneous equations are created for various 
diameters Di using signal ratios calculated from the 
appropriate scattering theory (Mie theory or non-spherical 
scattering theory) for various particle and dispersant 
refractive indices. These equations are then solved for the 
constants Al, A2, A3, Bl, B2, B3,C1. Of course this process 
can be extended to more than 3 angles and for polynomial order 
greater than 3 . 

Particles which are two small for single particle counting 
may be measured by stopping the flow and using the heterodyne 
signal of the scattered light to measure the size distribution 
from the Brownian motion of the particles. The particle size 
distribution is determined by inverting either the power 
spectrum or the auto-correlation function of the Doppler 
broadened scattered light from the moving particles- The 
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particle size distribution from Browniaii mofcion can also be 
used to determine the effective particle/dispersaht refractive 
index (scattering model) by measuring the hydrodynamic size of 
a particle along with the scattering signal amplitudes. The 
scattering model can be determined from the scattering 
intensity at each angle, and the true size for a 
representative single particle or a group of particles. The 
true size can be determined from the power spectrum or 
autocorrelation function of the heterodyne signal via Brownian 
motion, from the ratio of intensities of light scattered at 
two angles in the size region where the ratio is an accurate 
indicator of size, or by other size measurement techniques. 
This scattering model could then be used for computations of 
particle size in a counting process, which does not use 
Brownian motion. 

Other methods of generating the heterodyne local oscillator 
are also claimed in this disclosure for systems like in Figure 
19. For example a small reflecting sphere or scattering object 
could be placed in the interaction volume to scatter light 
into the heterodyne detectors along with the light scattered 
by the moving particle. Since this sphere or object is 
stationary, the optical phase difference between the scattered 
light from the moving particle and light scattered (or 
reflected) from the sphere or object would increase as the 
particle passed through the beam,; creating an .oscillating beat 
scatter signal on the detectors, at high frequencies. Then the 
local oscillator beam, which passes through lens 6, could be 
eliminated. 

In Figure 19, the effective scattering angles seen by each 
detector can depend upon the position of the particle in the 
beam. The addition of lenses 7 and 8, as shown in Figure 24 
(which only shows the detection portion of Figure 19) , will 
lower the scatter angle sensitivity to particle position. Each 
of these lenses place the detectors in a plane which is 
conjugate to the back focal plane of either lens 3 or lens 4. 
Essentially the back focal plane of lens 3 is imaged by lens 7 
onto detectors 3 and 4; and the back focal plane of lens 4 is 
imaged by lens 8. onto detectors 1 and 2. Also the detectors 
could be placed in the back focal planes of lens 7 and lens 8, 
where each point in the focal plane corresponds to the same 
scattering angle from any point in the interaction volume. 
This configuration nearly eliminates dependence of detector 
scattering angles on the position of particles in the beam. 

By powering the source at various intensity levels, the 
scattered light from particles which span a large range of 
scattering intensities can be measured with one analog to 
digital converter. Even though the dynamic range of the 
scattered light may be larger than the range of the A/D, 
particles in different size ranges can be digitized at 
different source intensity levels . The resulting signals can 
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be normalized to their corresponding source intensity and then 
used to determine the size of each particle. 

In Figure 11, the 2-Dimensional detector array could also be 
moved farther away from lens 4 to the plane which is nearly 
optically conjugate to the center of the sample cell. This may 
provide better imaging resolution of the particles on the 
array - 

Also it is recognized that many of the ideas in this 
disclosure have application outside of particle counting 
applications . Any other applications for these ideas are also 
claimed. In particular, the ideas put forth in Figures 21 and 
22 would also have application in ensemble particle sizing 
systems . 

. Many drawings of optical systems in this disclosure show 
small sources with high divergence which are spatially 
filtered by a lens and pinhole and then collimated by a second 
lens. In all cases, a low divergence laser beam could replace 
this collimated source^ as long as the spectral properties of 
the laser are appropriate for smoothing of Mie resonances if 
needed. 

Another issue is interf erometric visibility in the heterodyne 
signals described before. Misalignment of beamsplitter or 
lenses 5 or 6 in Figure 24 can lower the visibility of the 
heterodyne signals. Since this loss may be different on 
different detectors , the ratio of two signals may not be 
preserved. However, the ratio of the visibilities for two 
detectors will be the same for all particles. Therefore a 
correction for the effects of low visibility, for both 
absolute signals and signal ratios, can be determined by 
measuring scattered signals from one or more nearly mono-sized 
particles of known size and comparing the results with 
theoretical values to determine the effective visibility for 
each channel or visibility ratio for pairs of channels. This 
is most easily accomplished by measuring larger particles with 
scattering signals of very high signal to noise and looking at 
the actual heterodyne signals to determine the interf erometric 
visibility for each detector. This could be determined by 
blocking the local oscillator light and measuring the 
scattered signal pulse with and without the local oscillator 
to calculate the theoretical heterodyne signal from the 
measurement of the local oscillator power and the scattered 
pulse amplitude. Also simply comparing the ratio of two 
scattering heterodyne signals to the theoretical value for 
that particle size would also provide a correction factor for 
the ratio, directly. 

For particles which are much smaller than the size of the 
laser spot, the scattered signal for particles passing through 
various portions of the laser spot will be distributed over a 
range of peak amplitudes. For a group of monosized particles, 
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the probability that a peak amplitude Will be between value S- 
deltaS/2 and S+deltaS/2 is Ph{S)deltaS, where Ph(S) is the 
probability density function for scattering amplitude in 
linear S space. ^MeltaQ'' means the difference in Q between the 
end points of the interval in where Q may be S or Log(S) 
for example- For a group of monosized particles of a second 
size (diameter D2 in Figure 25), the scattering amplitudes, S, 
for particles that have passed through the same region of 
laser beam, as particles of the first size diameter Dl, are 
changed by a multiplier R and the probability density 
amplitude is changed by a multiplier of 1/R, as shown in first 
graph of Figure 25 for two particle diameters, Dl and D2 and 
the following equation: 

Pn{S) deltas = Pn (RS) deltaS/R 

If we switch to logarithmic space for S, we find that the 
probability density becomes shift invariant to a change in 
particle size (Pg(Log(S)) only shifts along the Log(s) axes as 
R or particle size changes . 

Using deltas = R^deltaLog ( S ) 

Pg(Log(S) )deltaLog(S) = Pg ( Log(R) +Log(S) )deltaLog(S) 

Where Pg(Log(S)) is the probability density function in Log(S) 
space. This shift invariance means that the differential 
count-vs . -Log (S) distribution, Cg, in logarithmic space is a 
convolution of the probability function Pg shown in Figure 25,' 
with the number-vs . -size Ng, where all are functions of 
Log (S) . 

Cg = Ng e Pg in convolution form where Pg is the response 
(impulse response) from a monosized particle enseinble 

Cg = Ng * Pgm in matrix form, where each column in matrix Pgm 
is the probability function for the size corresponding to the 
element of Ng which multiplies it. This equation can also be 
used when Pg is not a convolution form. 

These equations can be inverted to solve for Ng, given Cg and, 
Pg, by using deconvolution techniques or matrix equation 
solutions. Pg is determined theoretically from the laser beam 
intensity profile or empirically from the Cg measured for one 
or more monosized particle samples. If Pg has some sensitivity 
to particle size, the matrix equation is preferable. 

These relationships also hold for the above functions, when 
they are functions of more than one variable. For example, 
consider the case where Cg is a function of scattering values 
SI and S2 from two scattering detectors at different 
scattering angles . Then 

Cg(Log(Sl) ,Log(S2) )deltaLog(Sl)deltaLog(S2) is the number of 
events counted with log signals between Log (SI) - 
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delta (Log (SI) ) /2 and Log (SI) + delta (Log (sl) ) /2 , and between 
Log(S2) - delta (Log (S2) ) /2 and Log(S2) + delta (Log ( S2 )) /2 . Then 
Cg (Log (SI) , Log (82 ) ) could be plotted as a surface on the 
(Log (SI) , Log (S2 ) ) plane as shown in Figure 26. This surface is 
determined from the event density of the ^^scatter plot'' or 
^Mot plot'' of all of the particles on the Log (SI) , Log (S2) 
plane (each particle is represented by its values of Log (SI) 
on the X axis and Log(S2) on the Y axis on the scatter plot) . 
So an event is the dot or point in Log (SI) , Log (S2 ) space (or 
S1,S2 space) which represents a counted object. The 
distribution functions Pg and Cg are calculated from the 
number of events for each small area in this space. A group of 
monosized particles will theoretically produce a group of 
points in S1,S2 space which follow the function Log(S2) = 
Log (SI) + Log(R) . Hence the data points will line up along a 
line of slope =1 and with an offset of Log(R) . The 
distribution of points along the length of the line for the 
particle group is determined by range of SI and S2 for that 
group due to the intensity distribution of the source beam. If 
particles pass through the beam at random locations, the 
distribution of data points along the line will follow the 
intensity distribution along each of the SI and S2 axes. R, 
which is the ratio between S2 and SI, changes with particle 
size. As the particle size decreases below the wavelength of 
the source, R becomes a constant for all sizes, as determined 
from Rayleigh scattering theory. However, real measurements do 
not follow theory exactly due to structural imperfections in 
the optical system. These imperfections will cause broadening 
of the line. This broadening is illustrated by an elliptical 
shape (however the actual shape may not be elliptical) in 
Figure 26. Each ellipse represents the perimeter around a 
group of counted data points on the S1,S2 plane from particles 
all of one size. 

Notice, if the only cause of response broadening is due to the 
intensity distribution of the source spot, then the ellipses 
in Figure 2 6 will become nearly circles, because both signals 
will come from similar intensity distributions. 

A group of monosized particles will produce a differential 
count distribution in Sl,S2 or Log ( SI ) , Log ( S2 ) space. In each 
case, a differential count distribution from a polysized 
sample will be the sum of the monosized distributions, each 
weighted by the percentage of particles of that size in the 
total distribution. Hence the particle number-vs . -size 
distribution can be determined by inverting this total 
differential count distribution, as a function of SI and S2, 
or Log(Sl) and Log(S2), using algorithms which may include 
those already developed for image restoration. In 
Log (SI) , Log (S2) space the monosized response functions will be 
similar in shape over a large size range, because the 
functions are approximately shift invariant to size over the 
Log (SI) , Log (S2) space. In this logarithmic space, 
deconvolution can be used to invert the count distribution in 
either one or more dimensions. The signal pulse from each 
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event may pass through ah analysis or sorting as described 
before, to sharpen the monosized response for higher size 
resolution. These pre-processed pulses are counted vs. a 
parameter (SI, S2, etc) such as peak value, total area, total 
correlated signal, etc., using the methods outlined 
previously. Each counted event is placed into the S1,S2 or 
Log(Sl) ,Log(S2) space, where SI and S2 may be the pulse peak 
value, pulse area, or any of the other size related parameters 
which can be calculated from the scattering signals. Then this 
space is broken up into very small regions, and the events in 
each region are summed to give sampled values of the 
differential count-vs - -parameter distribution in the 2 
dimensional space. The known monosized response (in the 2 
dimensional space), which may be size dependent, is used to 
invert this differential distribution to produce the particle 
number vs. size distribution. This monosized response may be 
calculated from scattering theory and the optical design 
parameters, or it may be measured empirically by, recording the 
differential event count distribution in the space from 
monosized particle groups of known sizes. The known monosized 
response defines a region in the space, where scattering from 
single particles can produce counts as shown in Figure 27. 
Events which are outside. of this region, can be rejected as 
non-particle events (i.e. non-single particle events or 
noise) , which may be due to multiple particles in the 
interaction volume or noise. This is particularly important 
for small sized particles with low scattering signals, where 
detection noise can cause many non-particle counted events as 
shown in Figure 28. These noise events can also be included 
into the monosized response functions. For example, when a 
group of monosized particles are measured to produce an 
empirical monosized response function, many noise events, will 
be measured. These noise events can be included in the 
monosized response function so that they are removed as part 
of the inversion process, when that function is used as one of 
response set in the inversion algorithm- The complete 
monosized response function set can be generated from 
scattering data from only a few well chosen, monosized particle 
groups. The intervening response functions are interpolated 
from the trend of the theoretical scattering/optical system 
model. The empirical data from monosized samples may only be 
needed to locate the theoretical model in the space . The power 
of this process is that both the absolute signal data, which 
is needed to size particles in the Rayleigh scattering regime, 
and the scattering signal ratio information are combined into 
one space, where non-particle events are easily identified. 
This process can be applied to data taken in any number of 
dimensions, from one scattering angle to any number of 
scattering angles . Also any dimension of this process can be 
represented by a scattering signal related parameter (peak, 
integral, etc.) or combination of scattering signals (ratio of 
S2/S1, correlation between S2/S1, etc.) - Higher niomber of 
dimensions provides better discrimination against non-particle 
events, but with added cost of more detectors and computer 
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processing time. For example, more than two detectors could be 

placed behind each of slit 1 or slit 2 (aperture 1 or aperture 
2), in the previous figures, to provide additional dimensions 
to the problem. For 3 detectors you could plot each event on 
the S3/S1, S2/S1 plane. Then the efffects of source spot 
intensity variations would be reduced and the impulse 
responses (the ellipses shown previously) would become very 
localized, perhaps eliminating the need for deconvolution . 
All four detectors from these figures (detectors 1, 2, 3, 4) 
could also be combined into one four dimensional space as 
described in this section or into two 2-dimensional spaces, 
which are first solved separately and then the results are 
combined into one final size distribution by blending methods 
described previously. 

If all of the particles in a particular sample are in a size 
region where the signal ratio is not sensitive to particle 
size, such as the Rayleigh size regime, the scattering model 
could be determined empirically from dynamic scattering 
measurements. If the particle flow is stopped, the 
heterodyning detection system can measure the Doppler spectral 
broadening due to Brownian motion (dynamic light scattering) . 
The particle size distribution from this measurement may be 
used directly, or the optical scattering model may be 
determined from the dynamic scattering size distribution to 
invert the absolute scattering signal amplitudes from the 
count-vs . -scattering signal distribution. In this way, the low 
size resolution distribution from dynamic light scattering 
will provide scattering model selection for the higher size 
resolution counting method. This technique can be used over 
the entire size range of the dynamic light scattering to 
select the scattering model for counting particles inside or 
outside of the size range of dynamic light scattering. The 
scattering model may also be determined by inverting the count 
distribution in S1,S2 or Log (SI) , Log (S2) space. This inversion 
will create a line function in the space. The shape of this 
line function in the size transition from Rayleigh scattering 
(where the ratio between SI and S2 is constant) to larger 
particles will indicate the scattering model and refractive 
index of the particles. 

This multi-parameter analysis also provides for separation of 
mixtures of particles of different compositions such as 
polymer particles mixed with metal particles or polymer^ 
particles mixed with air bubbles. Hence, the count of air 
bubbles could be eliminated from the count distribution. 
Figure 27B shows the methodology. Particles of different 
composition will have different response profiles in the 
multi-dimensional space. And the count events will be grouped 
to follow different profiles for each particle composition. So 
the data points (events) for particles of different 
composition will occupy different response profiles as shown 
in Figure 27B. Individual particle size distributions and 
particle concentrations, for each particle type, could be 
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determined from analysis of this data lifeihg the techniques 
described in this disclosure, individually for ekch response 
profile. Hence, a different particle refractive index and 
optical scattering model would be determined and used to 
calculate the size distributions for particles in each 
composition group, separately. These techniques could also be 
extended to larger numbers of dimensions by measuring more 
signals . 

This process could also be used by replacing signals, in these 
multiple parameter plots, with ratios of signals. Any of these 
multiple angle configurations may be extended to many more 
angles simply by adding more scatter detectors which view the 
same interaction volume. For example, consider 4 such 
detector signals. Si, S2, S3, S4. Each of these signals could 
be pulse peak, pulse area, or correlated peak, etc. The ratios 
S4/S1, S3/S1, and S2/S1 are plotted in 3-dimensional space, 
one point for each particle counted. These ratios could also 
be S4/S2, S3/S2, S2/S1, etc. The point here is that the Mie 
resonances cause the scattering signals at various angles to 
oscillate together vs. particle size. For any size range, 
there is always a region of scattering angles where the ratio 
of scatter from two different ranges of angle are nearly 
independent of Mie resonances and particle composition. The 
path of these ratios in 2~dimensional space (e.g. S3/S1,S2/S1) 
or 3-dimensional space (e.g. S4/S1, S3/S1, S2/S1) , are only 
weakly dependent on particle composition. The strongest 
particle composition dependence is for spherical particles in 
the size region of the Mie resonances- When thousands of 
particles are measured, their points will follow a multi- 
dimensional curve or line, in this multi-dimensional space, 
which indicates the sphericity or composition of the 
particles. This multi-dimensional line is formed by the 
highest concentration of points in this multidimensional 
space. Only points which are within a certain distance of this 
line are accepted as true particles. The outliers represent 
particles which passed through the edge of the source beam or 
whose signals are contaminated by noise. Also non-single 
particle events, such as noise pulses or multiple particles, 
would also be rejected because their combination of 
coordinates in this space would not agree with the possible 
coordinates of a particle. The signal ratios could also be 
replaced by the signal values to take advantage of absolute 
signal information, which is particularly advantageous in the 
Rayleigh region where signal ratios are weakly dependent on 
particle size, but absolute signal levels are strongly 
dependent on particle size. For particles of size above the 
Rayleigh region, the signal ratios may be preferred because 
the spread of particle events around the multi-dimensional 
line is very small for signal ratios which remove the 
dependence upon the particle position in the beam, removing a 
portion of the monosized response broadening shown in Figures 
27 and 28. Both of the techniques, described here and in the 
description of Figures 27 and 28, will be needed to cover the 
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entire size range, because the signal ratid's arei not strongly 
size dependent in the Rayleigh region (for particles below 0.1 
micron in visible light) . Another op)tidn would be to use a 
multi-dimensional space where some 'dim:ehsi6hs wiere signal 
ratios and other dimensions where absolute signals. Then the 
monosized response would only be strongly broadened in the 
absolute signal dimensions. For small particles, the absolute 
signal dimension would be processed with deconvolution and 
noise event rejection as shown in figures 27 and 28; and for 
larger particles, signal ratio dimensions would be used to 
, determine size with outlier rejection and minimal 
deconvolution. The entire path of the line constructed by the 
particle events determines the optical scattering model. This 
is particularly important for the Rayleigh region where the 
particle refractive index must be determined to calculate the 
dependence of absolute scattering amplitude on particle size. 
In any case, there will be a line or curve, in multi- 
dimensional space, which follows the path defined by the 
highest concentration of counted events. This line will define 
the optical scattering model and the particle shape or 
composition, if the user cannot provide that information. Size 
accuracy may be improved by using any apriori information 
about the particles to determine the scattering model from 
theoretical models . For example, the path of the multi- 
dimensional line could be calculated from scattering theory, 
given the particle refractive index and shape, but sometimes 
this information is not well known. If the particle 
composition or shape is unknown, then this empirically 
determined line in multidimensional space is compared to the 
theoretical lines for particles of various compositions and 
shapes. The theoretical line, which most closely matches the 
measured line from the unknown particle dispersion, is assumed 
to represent the composition and shape of the unknown 
particles . 

The accuracy of the process described above improves as more 
scattering angles are measured. For example, the measured 
values of the scattered light for each of three scattering 
angles could be measured for each particle. These data points 
are then analyzed in a three dimensional scatter or dot plot. 
A line could be generated in 3 dimensional space by 
determining the path where the maximum concentration of 
particles (or dots in the plot) reside. In any one axis, this 
line may be multi-valued vs. particle diameter, especially in 
the region of Mie resonances. However, the line will not be 
multi-valued in 3 dimensional space. The spread of points 
about this line will be determined by the intensity 
distribution of the source beam in the interaction region. 
This group of points could be deconvolved in 3 dimensional 
space to produce a more sharply defined set of points, with 
less spread from the line, providing better size resolution 
along the line. But a better solution is to measure 4 
scattering values at 4 different scattering angles for each 
particle . And then take ratio of each of any 3 values with the 
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fourth value (or any other value) to remove the effect of 
intensity variation for particles which pass through different 
portions of the beam. Produce a scatter plot of these 3 ratios 
in three dimensions, where each point in 3 dimensional space 
is placed in Xm, Ym, and Zm values corresponding to the three 
ratios for each particle. Since the intensity distribution 
broadening is reduced, most of the points will tightly follow 
a line in three dimensional space. Outliers which are not 
close to the line passing through the highest concentration of 
data points may be eliminated as not being real single 
particles. The remaining data points (Xm,Ym, Zm) are then 
compared to different theoretical models to determine the 
composition and/or shape of the particles. The 3 dimensional 
function which describes the theoretical scattering is Zt 
where Zt is a function of Xt and Yt : 

Zt = Zt (Xt, Yt) 

Let (Xm,Ym, Zm) be the set of data points measured from the 
counted particles. Where the values in the X,Y,Z coordinates 
represent either absolute scattering signals SI, S2, S3 or 
signal ratios S4/S1, S3/S1, S2/S1 (or any other combination of 
ratios) . Then define an error function Et for a certain 
theoretical model as : 

'Et(Xm,Ym) = (Yt (Xm) -Ym) -^2 for Xm in the region Xmy where 
Yt (Xm) is single-valued 

Et(Xm, Zm) = ( Zt (Xm) -Zm) ^"2 for Xm in the region Ziriy where 
Zt(Xm) is single-valued 



Where Yt(Xm) is the theoretical value of Yt at Xm and Zt (Xm) 
is the theoretical value of Zt at Xm. Then find the 
theoretical model which produces the minimum sum of Esum over 
all values of Xm in the data set. 

Esum = SUM(Et (Xm, Ym) ) /Ny + SUM (Et (Xm, Zm) ) /Nz 

Where Ny is the number of points in region Xmy and Nz is the 
number of points in Xmz . 

And SUM is the sum of Et over its valid region of Xmy or Zmy. 
Esum is calculated for various theoretical scattering models, 
for spherical and non-spherical particles, and the model with 
the lowest Esum is chosen as the model for the sample. The sum 
of Esum values from multiple particle samples can also be 
compared for different theoretical models. The model with the 
lowest sum of Esum values is used to analyze all of those 
samples of that type. This calculation may be computationally 
intensive, but it only needs to be done once for each type of 
sample. Once the optimal theoretical model is determined for 
each particle sample type, the appropriate stored model can be 
retrieved whenever that sample type is measured. The chosen 
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theotetical model will provide the particle diameter as a 
function of Xm, Ym, and Zm for each detected object. 

Signal ratios show reduced sensitivity to' the position of the 
particle in the source beam because each scattering signal is 
proportional to the optical irradiance on the particle. 
Usually to obtain optimal signal to noise, a laser source will 
be used to provide high irradiance but with lower irradiance 
uniformity due to the Gaussian intensity profile. The 
broadening in the monosized response, as shown in Figures 27 
and 2 8 for example, may be reduced by insuring that only 
particles which pass near to the peak of the beam intensity 
profile are counted. Many methods have been described above to 
accomplish this selection process. Other methods could include 
the use of small capillaries or sheath flow to force all of 
the particles to go through the center of the source beam. But 
these methods are sometimes prone to particle clogging. In 
sheath flow, the particle dispersion is restricted to flow 
through a narrow jet, which is surrounded by a flow sheath of 
clean dispersant. If the particle concentration is low, the 
particles in this narrow stream will pass through the laser in 
single file and in locations close to the center of the beam. 
This method could be used with the ideas in this disclosure, 
but the wide range of particle sizes would require many 
different sized jets to handle the entire size range with the 
constant .danger of clogging. The methods described in this 
disclosure can be used within a flow system of much larger 
dimensions, because the optical system only views and counts 
particles within a small interaction volume of that much 
larger volume. Particles which pass through that volxome and 
which are outside of the size range for that measurement 
system will produce data points in the multidimensional space 
which are far from the multi-dimensional line of the optical 
model. They may be rejected based upon this criteria or simply 
based upon the length of the scattering pulse. The small 
particles are counted and sized by the higher angle system. 
The larger particles are sized by the 2-dimensional array or 
lower angle scattering systems. These independent particle 
size distributions are then combined to produce one size 
distribution over the full size range of the instrument. 

The measurement of particle shape has become more important in 
many processes. Usually the shape can be described by length 
and width dimensions of the particle. If the length and width 
of each particle were measured, a scatter plot of the counted 
particles may be plotted on the length and width space to 
provide useful information to particle manufacturers and 
users, and this type scatter plot is claimed in this 
invention. If the particles are oriented in a flow stream, the 
angular scattering could be measured in two nearly orthogonal 
scattering planes, one parallel and one perpendicular to the 
flow direction. Each of these scatter detection systems would 
measure the corresponding dimension of the particle in the 
scattering plane for that detection system. If the flow of 
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particle dispersion flows through a irestrictioh, so as to 
create an accelerating flow field, elongated particles will 
orient themselves in the flow direction. Figure 2 9 shows one 
of these scatter, detection systems where the scattering plane 
is parallel to Ys and measures the projected particle 
dimension in the Ys direction, which is parallel to the 
projection of the flow direction of the particles in the Ys/Xs 
plane. A second scattering detection system could be placed in 
a scattering plane which includes the Z axis and Yp as shown 
in Figure 30. This detection system would measure the particle 
dimension in the plane perpendicular to the flow. Each 
particle is counted with two dimensions, one parallel to and 
the other perpendicular to the flow, as measured concurrently 
by these two detection systems- In some cases, the particles 
cannot be oriented in the flow and they pass through the beam 
in random orientations. The detection configuration in Figure 
31 shows three scattering systems. Each system is in a 
scattering plane which is approximately 12 0 degrees from the 
next one. If the particle shapes are assumed to be of a 
certain type: rectangular, ellipsoid, etc . , three size 
measurements in various scattering planes can be used to solve 
for the length, width, (or major and minor axis, etc.) and 
orientation of each particle. These planes can be separated by 
any angles, but 120 degrees would be optimal to properly 
condition the 3 simultaneous equations , formed from these three 
size measurements. 

When measuring larger particles / which require smaller 
scattering angles, the scatter collection lens' may be centered 
on the Z axis, with scattering detectors in the back focal 
plane of the collection lens, as shown in Figure 32. As shown 
before, lens 1', pinhole 1, and lens 2 are not needed if a 
spatially clean collimated beam, such as a clean laser beam, 
is incident on the particle dispersion in the sample cell. 
Lens 3 collects scattered light from the particles and focuses, 
it onto a group of detectors in the back focal plane. As 
before, the length and width of each randomly oriented 
particle is determined by 3 independent size measurements or, 
in the case where the size measurements are not independent, 
you must solve a set of simultaneous equations as described 
below. If the particles are oriented in the flow, only the Ys 
direction (parallel to the flow) and a set of detectors in the 
direction perpendicular to Ys are needed. As before, these two 
directions can be at any angle, but parallel and perpendicular 
to the flow are optimal. In the random orientation case, each 
measurement is made by a separate arm of the detector set in 
the three directions Ys, Yl, and Y2 . These directions can be 
separated by any angles, but 12 0 degrees (see Figure 33) would 
be optimal to properly condition the 3 simultaneous equations 
formed .from these three size measurements. The scatter 
detector signals in each direction (or scattering plane) are 
combined by ratio of signals or other algorithms to determine 
the effective size in that direction. Then three simultaneous 
equations are formed from these size measurements to solve for 
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the width, length, and orientation of each particle. Figure 34 
shows how this detector configuration is used in the system 
from Figure 11. And as indicated above, the scatter in various 
scattering angular ranges can be measured each scattering 
plane by a separate optical system, as shown in figure 29 for 
example, in each scattering plane. The scattering plane is the 
plane which includes the center axis of the source beam and 
the center axes of the scattered light beams which are 
captured by the detectors . 

The accuracy of the methods outlined above is improved by 
solving another type of problem. The sizes calculated from 
angular scattering data in each of two or more directions are 
not usually independent. In order to accurately determine the 
shape parameters of a particle^ the simultaneous equations 
must be formed in all of the scattering signals. The form of 
the equations is shown below: 

Si = Fi (W,L,0) 

Where Si is the scattering signal from the ith detector. In 
the case of three directions (or scattering planes) and three 
detectors per direction, we have 9 total detectors and i = 1, 
2, , 9) 

W is the ^Vidth" parameter and L is the ^^length" parameter of 
the particle. In the case of a rectangular shape model, W is 
width and L is length. In the case of an ellipsoidal model, W 
is the minor axis and L is the major axis, etc. O is the 
orientation of the particle which could be the angle of the 
particle's major axis relative to Ys, for example. The 
functions Fi are calculated from non-spherical scattering 
algorithms and the form of Fi changes for different particle 
shapes (rectangles, ellipsoids, etc.). These equations. Si = 
Fi, form a set of simultaneous equations which are solved for 
W, L and O for each particle. If the Fi functions do not have, 
a closed form, iterative methods may be employed where the 
Jacobian or Hessian are determined by numerical, rather than 
symbolic, derivatives. Also the closed form functions for Fi 
could be provided by fitting functions to Fi(W,L,0) calculated 
from the non-spherical scattering algorithms. 

If we had two detection angles per each of three scattering 
planes, we would have 6 equations with 3 unknowns. With three 
detectors per scattering plane the size range may be extended 
and we will have 9 equations with 3 unknowns. For particles 
with more complicated shapes, such as polygonal, more 
scattering planes may be required to determine the particle 
shape parameters. In any case, a shape model is assumed for 
the particles and the set of equations Si == Fi are created for 
that model where Fi is a function of the unknown size 
parameters and Si is the scattered signal on detector i. This 
method can be applied to any of the shape measuring 
configurations shown before. This technique can also be 
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applied to ensemble size measuring systefrts when: the particles 
all have the same orientation as in accelerating flow. This 
invention claims scattering measurements from any number of 
angular ranges, in any number of scattering planes. 

Low scattering signals from small particles may be difficult 
to detect. Figure 35 shows another variation where the source 
beam is passed through a patterned target which is conjugate 
to the interaction volume . The image of the target occurs in 
the interaction volume which is defined by aperture 1 or 
aperture 2 . This target could consist of a sinusoidal 
transmission- pattern or a Barker code pattern. As the 
particles pass through the image of this pattern, the 
scattered light is modulated by the modulated source intensity 
distribution in the interaction volume and so the. scattered 
signal-vs . -time distribution is equivalent to the spatial 
intensity distribution. For a sinusoidal pattern, a phase 
sensitive detector with zero degree and quadrature outputs 
could be used to detect the sinusoidal signal of arbitrary 
phase. For a given particle velocity, the scattered signal 
could be filtered by a bandpass filter which is centered on 
the frequency equal to the particle velocity divided by the 
spatial wavelength of the sinusoidal intensity distribution in 
the interaction volume. The phase sensitive detector reference 
signal would also match this frequency. Better signal to noise 
may be achieved with other types of patterns- A Barker code 
target pattern will produce a single peak with very small side 
lobes when the scattering signal is correlated with a matching 
Barker code signal using a SAW or CCD correlator. When two 
scattering signals are multiplied and integrated, the ^ zero 
delay (tau) value of the correlation function is obtained. 
This value will have the lowest fluctuation when the two 
signals have strong correlation as when both signals are from 
the same particle, instead of uncorrelated noise. The 
integrated product of the two signals will show less noise 
than the separated integrated signals. So the product of 
signals from two different angular ranges or the integral of 
this product over the particle pulse period will provide a 
signal parameter which is less sensitive to noise and which 
can be substituted for Si in any of the analyses described 
above. Figure 35 also illustrates an additional scattering 
detector on aperture 2 for detection of three scattering 
angles. This can also be extended to a larger number of 
detectors . 

As shown before, ratios of scattering signals can be analyzed 
as a multi-dimensional function. Another method is to look at 
the individual signal ratios vs. particle diameter as shown in 
Figure 36 for the case of three signal ratios. Any real 
particle event should produce a point on each curve which 
align vertically at the same diameter- Each curve indicates 
the particle size, but the most accurate size is determined 
from the curve where its point is in a region of high slope 
and monotonicity- For any particle, the 3 measured ratios 
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would determine the approximate particle size region and allow 
selection of the one ratio which is in the region of hi-ghest 
slope vs. particle size and is also not in a multi-valued 
region (caused by Mie resonances) . This selected ratio would 
then be used to determine the precise particle size for that 
particle . 

The ratio of scattering signals from different scattering 
angles reduces the dependence of the particle size 
determination on the particle path through the light beam. 
Particles with signals below some threshold are eliminated 
from the count to prevent counting objects with low signal to 
noise. The accuracy of counts in each size bin will depend 
upon how 

uniform this elimination criteria is over the entire size 
range. Many methods have been described in this disclosure for 
reducing this problem. These methods are improved by having a 
source beam with a ^^flat top'' intensity , distribution and very 
sharply defined edges. This flat top intensity distribution 
can be provided by placing an aperture in an optical plane 
which is conjugate to the interaction voliome or by using 
diffractive optic or absorption mask beam shapers . Another 
technique which will accurately define an interaction volume 
is shown in Figure 37. No selection criteria is required for 
the direction which is parallel to the particle flow 
direction, because in this direction each particle passes 
through a similar intensity distribution and digitized signal 
values may be analyzed to find maximum or the integral for 
each particle signal. The primary criteria for eliminating 
particles from the count is based upon the position of the ■ 
particle along the axis which is perpendicular to particle 
flow direction. The position of the particle along the 
direction perpendicular to the particle flow and the 
scattering plane {y direction) can be determined by using a 3 
element detector which is in an optical plane which is 
conjugate to the interaction volume as shown in Figure 37. 
This figure shows the position and orientation of the 3 
element detector in the optical system and an enlarged view of 
the detector elements showing the path of various particles 
passing through the interaction volume as seen by the detector 
elements. The beamsplitter, following lens 3, splits off some 
scattered light to the 3 element detector. By measuring the 
signal ratio between elements 1 and 2 or elements 3 and 2, the 
y position of the each particle is determined and only 
particles within a certain y distance from the center of the 
light beam are accepted. This signal ratio criteria is 
extremely accurate and uniform among all particle sizes so 
that proper mass balance is maintained over the entire size 
range. The ratio is also insensitive to how well the particle 
is optically resolved because the fraction the particle image 
on each the two detectors spanning the image is not strongly 
dependent on the size or sharpness of the image, but is 
strongly dependent on the y position of the particle. Figure 
37 shows the 3 element detector in a heterodyne arrangement. 



59 



wo 2005/091970 



PCT/US2005/007308 



with a portion of the source light being mixed With the 
scattered light. However, this idea is also applicable to non- 
heterodyne configurations by just removing the beamsplitter 
, between lens 1 and lens 2, 

Many figures (Figure 19 for example) show the heterodyne 
system with a negative and positive lens pair (lenses 5 and 6) 
which provide a local oscillator beam which matches the 
wavefront of scattered light from the particles. Figure 38 
shows an alternative design where all beams are nearly 
collimated in the regions of the beamsplitters. This 
configuration may be easier to align and focus. Lens 1 
collimates the source light which is sampled by beamsplitter 1 
and directed to beamsplitter 2 by the mirror. Lens 3 
collimates the light scattered by the particle and this light 
is combined with the source light by beamsplitter 2 and 
focused through aperture 2 by lens 7. As before, aperture 2 is 
conjugate to the particle interaction volume and defines the 
interaction volume, in the sample cell, which is viewed by the 
detectors 3, 4 and 5. Usually, the focal length of lens 2 is 
long to provide a source beam of low divergence and the focal 
length of lens 3 is short to span a large range of scattering 
angles. If the light source is a laser diode, without 
anamorphic optics, the major axis of the intensity 
distribution ellipse at the interaction volume should be in 
the plane of the particle flow and scattering plane to provide 
a long train of heterodyne oscillations for signal detection 
and to provide the lowest beam divergence in the scattering 
plane. A circular source beam may require anamorphic optics to 
create an elliptical beam in the interaction volume to provide 
the advantages mentioned above. However, the advantages of the 
ideas described in this disclosure can be applied to a source 
beam with any intensity distribution. 

The matching of light wavefronts between the source beam and 
scattered light at the heterodyne detectors is important to 
-maintain optimum interf erometric visibility and maximum 
modulation of the heterodyne signal on each detector. Since 
perfect wavefront matching is not achievable, the 
interf erometric visibility must be determined for each 
detector to correct the signals for deviation from theoretical 
heterodyne modulation amplitude. The visibility is determined 
by measuring particles of known size and comparing the 
heterodyne signals to the signals expected from theory. To 
first order, the interf erometric visibility should be 
independent of particle size for particles much smaller than 
the source beam in the interaction volume. The visibility 
could be measured for particles of various sizes to measure 
any second order effects which would create visibility 
dependence on particle size. If only signal ratios are used 
for determining size determination, only the ratios of 
interf erometric visibility need to be calibrated by measuring 
scattering from particles of known size. 
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The number of cycles in the heterodyne modulated pulise is 
determined by the length of the trajectory of the particle 
through the source beam. The frequency of the heterodyne 
modulation is determined by the velocity of the particle 
through the beam. In general the power spectrum of the signal 
will consist of the spectrum of the pulse (which may be 10 KHz 
wide) centered on the heterodyne frequency (which may be 1 
MHz) . Both of these frequencies are proportional to the 
particle velocity. Actually the best frequency region for the 
signal will be determined by the power spectral density of the 
detector system noise and/or the gain-bandwidth product of the 
detector electronics. For this reason^ in some cases the 
particle flow velocity should be lowered to shift the signal 
spectrum to lower frequencies . The particle concentration is 
then adjusted to minimize the time required to count a 
sufficient number of particles to reduce Poisson statistic 
errors. This is easily accomplished for small particles which 
usually have higher count per unit volume and require lower 
noise to maintain high signal to noise - 

In cases where optical heterodyne detection is not used, the 
signal to noise may be improved by phase sensitive detection 
of the scattered light- Modulation ' of the optical source may 
provide for phase sensitive detection of the scattering 
signal. The source is modulated at a frequency which is much 
larger than the bandwidth of the signal. For example, consider 
a source modulated at 1 megahertz with a scattering pulse 
length of 0.1 millisecond. Then the Fourier spectrum of 
scattered signal pulse would cover a region of approximately 
10 KHz width centered at 1 megahertz.. If this signal is 
multiplied by the source drive signal at 1 megahertz, the 
product of these two signals will contain a high frequency 
component at approximately 2 megahertz and a difference 
frequency component which spans 0 to approximately 10 KHz. In 
order to eliminate the most noise but preserve the signal, 
this product signal could be filtered to transmit only the • 
frequencies "contained in the scattering pulse, without 
modulation (perhaps between 5000 and 15000 Hz) . This filtered 
signal product will have higher signal to noise than the raw 
signal of scattered pulses. This signal product can be 
provided by an analog multiplier or. by digital multiplication 
after both of these signals (the scattering signal and the 
source drive signal) are digitized. This product is more 
easily realized with a photon multiplier tube (PMT) whose gain 
can be modulated by modulating the anode voltage of the PMT. 
Since the PMT gain is a nonlinear function of the anode 
voltage, an arbitrary function generator may be used to create 
PMT gain modulation which follows the modulation of the 
source. The voltage amplitude will be a nonlinear function of 
the source modulation amplitude, such that the gain modulation 
amplitude is a linear function of the source modulation 
amplitude. An arbitrary function generator can generate such a 
nonlinear modulation which is phase locked to the source 
modulator. 
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As described at the beginning of this disclosure, multiple 
sized beams can also be used to control the effects of seeing 
more than one particle in the viewing voliime at one time. The 
key is to choose the proper scattering configuration to 
provide a very strong decrease of scattering signal with 
decreasing particle size. Then the scattering signals from 
smaller particles do not effect the pulses from larger ones. 
For example, by measuring scattered light at very small 
scattering angles, the scattered light will drop off as the 
fourth power of the diameter in the Fraunhofer regime and as 
the sixth power of diameter in the Rayleigh regime- In 
addition, for typically uniform particle volume vs. size 
distributions, there are many more smaller particles than 
larger ones . The Poisson statistics of the counting process 
will reduce the signal fluctuations for the smaller particles 
because individual particles pulses will overlap each other 
. producing a uniform baseline for the larger , particles which 
•pass through as individual pulses. This baseline can be 
subtracted from the larger particle pulse signals to produce 
accurate large particle pulses. A typical optical 
configuration is shown in Figure 41. 

Figure 41 shows such an optical system where the light source 
is spatially filtered by lens 1 and pinhole 1. Lens 2 
collimates and projects the source beam through the particle 
sample and an aperture mask, which is imaged onto the detector 
array by lens 3. Lens 2 could also focus the beam into the 
sample cell to increase light intensity. An annular spatial 
mask is placed in the back focal plane of lens 3 to only pass 
scattered light over a certain range of scatter angle as 
defined by the inner and outer radii of the annular filter, 
which is similar to mask 1 shown in figure 12. The very low 
angle scattering and incident beam are blocked by central stop 
of the annular aperture in the back focal plane of lens 3. 
Hence the detector array 1 sees an image of mask apertures and 
each detector element measures the scattered light only from 
particles in it's corresponding mask aperture over the angular 
range defined by the aperture (or spatial mask) in the back 
focal plane of lens 3. Each detector element is equal to or 
slightly larger than the image of the corresponding mask 
aperture at the detector plane, as long as each detector only 
sees the light from only it's corresponding mask aperture. A 
' beam splitter splits off a portion of the light to a second 
annular filter (in the back focal plane of lens 3) and 
detector array 2 . The angular ranges of the two annular 
filters are chosen to produce scattered values which are 
combined by an algorithm which determines the size of each 
particle. One such algorithm would be a simple ratio of the 
corresponding pulses from both arrays. And if the total 
scattered light is sensitive to particle composition, then the 
ratio of the two scattering signals can be used to determine 
the particle size more accurately. As with all other systems 
described in this disclosure, these ideas can be extended to 
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more than two detector arrays or more thain two' scattering 
angles^ simply by adding more annular spatial masks and 
detectors by using beamsplitters. And the signals can be 
processed and analyzed, using the methods described 
previously. 

This configuration allows each detector element to see scatter 
from only a certain aperture in the mask and over a certain 
scattering angle range determined by it's spatial mask. If 
the spatial mask defines a range of low scattering angles, the 
total scatter for the detectors viewing through that spatial 
mask will show a strong decrease with decreasing particle 
size. The signal will decrease at least at a rate of the 
fourth power of the particle diameter or up to greater than 
the sixth power of the diameter. Assuming the weakest case of 
fourth power, we can obtain a drop by a factor of 16 in signal 
for a factor of 2 change in diameter. This means that you need 
to control the particle concentration such that no multiple 
particles are measured for the smallest particle size measured 
in each aperture. The largest particle size which has 
significant probability of multiple particles in the aperture 
at one time should produce scattering signals which are small 
compared to the scattering from lower size measurement limit 
set for that aperture. However, this particle concentration 
constraint is relaxed if multiple pulses are deconvolved 
within a signal segment, as shown in Figures 71, 72 and 73. 

One annular filter aperture could also be replaced by a 
pinhole, which only passes the light from the source (the red 
rays) . Then the signals on each detector element would 
decrease as a particle passes through it's corresponding 
aperture at the sample cell. This signal drop pulse amplitude 
would directly indicate the particle size, or it could be used 
in conjunction with the other annular signals. No limits on 
the number of apertures in the sample cell mask or of annular 
filter/detector sets are assumed. More annular 

filter/detector sets can be added by using more beamsplitters. 

Figure 42 shows another configuration for determining particle 
size and shape. The light source light is focused through 
pinhole 1 by lens 1 and then focused into the sample cell by 
lens 2. The beam divergence and spot size in the sample cell 
are determined by the range of scattering angles to be 
measured and the size range of the particles. Essentially the 
spot size increases and the divergence decreases for larger 
particles. The scattered light is collected by lens 3, which 
focuses it onto many multi-element detector assemblies, which 
are in the back focal plane of lens 3. Each multi-element 
detector has multiple detector elements which measure a 
certain range of scattering angles along various scattering 
planes. Figure 42 shows an example with three scattering 
planes separated by approximately 120 degrees between adjacent 
planes. Howe^^-er, any number of scattering planes with any 
angular separation is claimed in this disclosure. Each multi- 
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element detector contains a central region which either 
captures of passes the source light so that it does not 
contaminate the measurement of the scattered light. The beam 
divergence will determine the size of the source light capture 
region on the multi-element detectors 

Each detector element has a shape which determines, how much 
of the scattered light at -each scattering angle is collected 
by the detector element. For example detector 1 has wedge 
shaped detectors which weights scattering angles 
progressively. Detector 2 has a higher order weighting, the 
larger scattering angles are gradually weighted more in the 
total signal for each detector element. These detector element 
shapes can take on many forms: rectangular, wedged, and higher 
order. Any shape, will work as long as the progression of 
collection width of the detector is different between the two 
multi-element assemblies so that when the particle pulse 
signals from the corresponding detector elements of the two 
multi-element detectors are ratioed, you obtain a ratio which 
is particle size dependent. The progression of the weighting 
function can also be defined by placing a variable absorbing 
plate, over each detector element, which varies absorption vs. 
radius r from the center of the detector assembly. This 
absorption plate can provide a weighting similar to that 
obtained by varying the width of the detector element vs. r. 
And since these size measurements are made in different 
scattering planes, multiple dimensions of each particle are 
determined separately. In general each detector element 
produces a signal Sab, where '^'^a'' is the multi-element detector 
assembly number and ■^b''' is the element number within that 
assembly. Then we can define Sab as: 

Sab = / wa(r) f(r,d) dr for the bth element in the ath 
assembly 

Here d is the dimension of the particle in the direction of 
the corresponding scattering plane. The scattered intensity at 
radius r from the center of the detector assembly 
(corresponding to zero scattering angle) for dimension of d is 
f{r,d). And wa{r) is the angular width (or weighting function) 
of the. detector element at radius r in assembly ^^a'', in other 
words the angle which would be subtended by rotating the r 
vector from one side of the element to the other side at 
radius r. For the simple 3 element assemblies shown in figure 
42, we obtain 6 measured values: 

Sll, S12, 813 for assembly 1 

S21, S22, S23 for assembly 2 

The signal on each detector element will consist of pulses as 
each particle passes through the beam. The Sab values above 
can be the peak value of the pulse or the integral of the 
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pulse or other signal values mentioned in this disclosure. For 
example, one possible case would be: 

Angular width for assembly 1 wl(r) == Ar 

Angular width for assembly 1 w2 (r) = Br*r 

For this case S11/S21, S12/S22, and S13/S23 are almost linear 
functions of the particle dimension in the direction of the 
corresponding scattering plane. These 3 dimensions can also be 
determined from an algorithm which uses all 6 S values by 
solving simultaneous equations which include the 
interdependencies of these values on each other, as described 
previously. In any event, the actual dimensions of the 
particle can be determined by assumption of a certain particle 
form such as rectangular, ellipsoidal, hexagonal, etc. More 
detector elements in each assembly will produce more accurate 
dimensions for randomly oriented particles . The true power of 
this technique is that the shape of each particle can be 
determined over a large size range by measuring only a few 
signals. Each element of each detector assembly could be 
broken up into sub-segments along the ^^r" direction to provide 
better size information by measuring the angular scattering 
distribution in each of the scattering planes. However, this 
may reduce the particle count rate because more digitizations 
and data analysis may be required per particle. 

These S values can also be analyzed using methods shown in 
Figures 2 6 through 28, by creating two dimensional plots of 
absolute measurements of SI* vs. S2* (Sll vs. S21, S12 vs. 
S22, and S13 vs. 323). Also a virtual 6 dimensional plot of 
all 6 values can be created (but not plotted) . Then the same 
methods can be used to eliminate particles which do not meet 
criteria for having passed through the central portion of the 
beam or which are caused by noise pulses (the absolute S 
values are not consistent with the size determined from the S 
ratios) 

The actual particle size system may consist of systems, each 
which is similar to the one shown in Figure 42. Each system 
would have a different source beam divergence and spot size in 
the sample cell to accommodate different size ranges. The 
count distributions from the systems are then concatenated (or 
blended as shown previously) into one total distribution over 
the entire size range of the product. For example, for 
rectangular or ellipsoidal particles, the width and length 
dimensions of each particle could be plotted on a ^^scatter^' 
plot to display the information in a useful format. 

For smaller particles, the source beam will be more focused 
(higher divergence and smaller spot size in the sample cell 
region) into the sample cell. This will help to define a 
smaller interaction volume for the smaller particles which 
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usually have higher number concentration than the larger 
particles . 

Figure 43 shows another version of this concept where the 
interaction volume for particle scattering is controlled by 
appropriately positioned apertures and by correlation 
measurements between signals from different scattering angles. 
The system is similar to that shown in Figure 42. But in this 
case additional apertures and lenses are added in the 
detection system.. Aperture 2 and aperture 3 are placed in 
optical planes which are conjugate to the source focused spot 
in the sample cell. These apertures are sized and oriented to 
only allow the image of the focused spot to pass on to the 
multi-element detector. The source beam in aperture 1 may have 
significant intensity variation so as to produce a large 
variation in scatter signals when particles pass through 
different portions of the source spot in the sample cell. In 
this case, the size of apertures 2 and 3 may be reduced such 
that their images, at the sample cell, only pass the uniform 
portion of the source beam intensity profile in' the sample 
cell. These apertures and lens 3 limit the volume, in the 
sample cell, from which scattered light can be detected by the 
multi-element detectors . Lenses 4 and 5 image the back focal 
plane of lens 3 on to the multi-element detector so that the 
detector sees the angular scattering distribution from the 
particles . The multi-element detectors 1 and 2 can also 
measure scattered light in the back focal plane of lens 4 and 
lens 5, respectively. Detector 3 collects all of the light 
that is scattered in a range of scattering angles which are 
defined by the annular aperture 4 . This detector provides a 
equivalent diameter based on equivalent spherical particle 
crossectional area, without shape dependence. In some cases 
the size as determined by the total scatter thorough aperture 
4 will be more accurate than the particle dimensions from the 
multi-element detectors. For example, detector 3 could be used 
to determine the particle area and the multi-element detectors 
could determine the aspect ratio of the particle using the 
ratio of the determined dimensions . Using the area and the 
aspect ratio, the actual dimensions could be determined- This 
may be more accurate than simply determining the dimensions 
separately using the multi-element detectors for particles 
whose major or minor axis may not line up with a scattering 
plane and which require the solution of simultaneous equations 
to determine the particle shape. 

Another feature of this design is the ability to use 
correlation or pulse alignment to determine which particle 
pulses are accurately measured and which pulses may be 
vignetted in the optical system. Figure 45 shows a crossection 
of the source beam focus in the sample cell. The outline of 
the beam is shown in red and the outline of the optical limits 
for scattered rays is shown in black. These optical limits are 
defined by the angular size of the detector elements and the 
size of aperture 2 or aperture 3. Two extreme scattered rays 
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are drawn in blue for scattered light at a particular 
scattering angle. The intersection (crosshatched area) of 
volume between those scattered rays and the source beam is the 
interaction volume in which a detector can detect scattered 
light at that scattering angle. For examp^le, consider the 
highest angle detector elements^ ic and 4C (see Figure 44) 
which are both in a scattering plane which is parallel to flow 
direction of the particles. The bottom portion of Figure 45 
shows an approximation to the interaction volumes for each of 
these detector elements. Notice that as particle A passes 
through these interaction volumes, both scattering signals 
from IC and 4C will be highly correlated, they will rise and' 
fall together with a large amount of overlap in time. However, 
particle B, which is farther from best focus, will show very 
poor correlation between these two detectors. In fact the 
pulses will be completely separated in time. This pulse 
separation between IC and 4C can be used to determine where 
the particle has passed through the interaction volume; and 
particles that are too far from best focus can be eliminated 
from the particle count. This correlation or pulse separation 
can be measured between any two detector elements, within a 
group (i.e. lA and IC) or between groups (i.e. IC and 4C) . 
Typically the scattering plane for detector groups 1 and 4 
would be parallel to the particle flow to obtain maximum 
delay. The correlation or pulse separation can be determined 
from the digitized signals using algorithms. However, this may 
.require very high speed analog to digital converters and 
enormous computational load to obtain a high particle count 
and size accuracy . Another solution is to use analog 
electronics to measure the correlation or the pulse separation 
as shown in Figure 46, where the P boxes are processing 
electronics which measure the pulse peak (peak detector) or 
pulse integral. The X box is an analog multiplier. And SIA, ' 
SIB, and SIC are the analog signals from detector elements lA, 
IB, and IC, respectively. The following equations will provide 
an estimate to the correlation between the pulses: 

R12 = P12/(P1*P2) 

R13 = P13/ (P1*P3) 

When R12 or R13 are small, the pulses have poor correlation 
and they should be eliminated from the count. 

The delay between any two pulses from separate detector 
elements can also be used to select valid pulses for counting. 
As the particle passes through the beam farther from best 
focus of the source beam, the delay between the pulses will 
increase. Some threshold can be defined for the delay. All 
pulse pairs with delays greater than the threshold are not 
included in the count. One example is shown in Figure 47, 
where the delay between pulses from detector elements IC and 
4C (see Figure 44) is measured to reject particles which pass 
through the beam too far from the source beam best focus. 
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Another criteria for pulse rejection is pulse width. As shown 
in Figure 45, particle B will produce a shorter pulse than 
particle A, because the detector element will only see 
scattered light from the particle while it is in the 
interaction volume (the crosshatched area) for that particular 
detector element. The pulses could be digitized and the pulse 
width would then be computed as the width at some percentage 
of the pulse peak height to avoid errors caused by measuring 
pulses of different heights. Any of these techniques discussed 
above can be implemented using digitization of the detector 
element signals and computation of parameters of interest from 
that digitized data or using analog modules which directly 
produce the parameter of interest (pulse delay, pulse width, 
correlation, etc.). While the analog modules may have poorer 
accuracy, they can be much faster than digitization and 
computation, allowing a larger particle count and better count 
accuracy. 

The pulse rejection criteria described above is used to reduce 
the number of coincidence counts by using apertures to limit 
the volume which is seen by the detectors. The interaction 
volume can also be limited by providing a short path where the 
particles have access to the beam as shown in Figure 48. Two 
transparent cones are bonded to the inner walls of the sample 
cell windows using index matching adhesive. The tips of each 
cone is cut off and polished to either a flat or a concave 
optical surface. The optical windows and transparent cones 
could also be replaced with solid cell walls with holes which 
are aligned to hollow truncated cones with optical windows on 
the truncated tip of each cone. This way the light travels 
through air, except for a thin layer of particle dispersion 
between the two windowed cone tips . The gap between the cone 
tip surfaces provides the only volume where the flowing 
particles can pass through the source beam and scatter light 
to the detectors. A dispersion with a large range of particle 
sizes will not clog this gap because the larger particles will 
flow around the gap and the particle concentration is very 
low. The cone tip surfaces can be tilted slightly so that the 
spacing between them is smaller on the side where the 
particles enter the gap and larger where they leave the. gap. 
In this way, particles larger than the minimum width of gap, 
but smaller than the maximum gap width, are prevented from 
jamming inside the gap. 

The beam focus may shift with different dispersant refractive 
indices due to refraction at the flat surface on the end of 
each cone- This shift in focus and angular refraction can be 
corrected for in software by calculating the actual refracted 
rays which intercept the ends of each detector element to 
define the scattering angular range of that element for the 
dispersant refractive. This correction is not needed for 
concave surfaces, on each cone tip, whose centers of curvature 
are coincident with the best focal plane of the source between 
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the two tips. Then all of the beam rays and scattered rays 

pass through the concave surface nearly normal to the surface 
with very little refraction and low sensitivity to dispersant 
refractive index. 

Another problem that can be solved by particle counting is the 
problem of background drift in ensemble scattering systems 
which measure large particles at low scattering angles. An 
ensemble scattering system measures the angular distribution 
of scattered light from a group of particles instead of a 
single particle at one time. This angular scattering 
distribution is inverted by an algorithm' to produce the 
particle size distribution. The optical system measures 
scattered light in certain angular ranges which are defined by 
a set of detector elements (see Figure 53) . Each detector 
element is usually connected to it's own separate electronic 
integrator, which is connected to a multiplexing circuit which 
sequentially samples each of the integrators which may 
integrate while many particles pass through the beam (see a 
portion of Figure 54) . So particle pulses cannot be measured 
in the ensemble system. 

The detector elements which measure the low angle scatter 
usually see a very large scattering background when particles 
are not in the sample cell. This background is due to debris 
on optical surfaces or poor laser beam quality. Mechanical 
drift of the optics can cause this background light to vary 
with time. Usually the detector array is scanned with only 
clean dispersant in the sample cell to produce background 
scatter signals which are then subtracted from the scatter 
signals from the actual particle dispersion. So first the 
detector integrators are scanned without any particles in the 
sample cell and then particles are added to the dispersion and 
the detector integrators are scanned a second time. The 
background scan data is subtracted from this second scan for 
each detector element in the array. However, if the background 
drifts between the two scans, a true particle scattering 
distribution will not be produced by the difference between 
these two scans. A third scan could be made after the second 
scan to use for interpolation of the background during the 
second scan, but this would require the sample cell to be 
flushed out with clean dispersant after the particles are 
present . 

A much better solution is to connect each of the detector 
elements, for the lowest angles of scatter, to individual 
analog to digital converters or peak detectors as shown before 
in this disclosure (see Figure 54) . Then these signals could 
be analyzed by many of the counting methods which are 
described in this disclosure. This would essentially produce 
an assemble/counting hybrid instrument which would produce 
counting distributions for the large particles at low 
scattering angles and deconvolved particle size distributions 
from the long time integrated detector elements at higher 
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scattering angles for the smaller particles. These 
distributions can be converted to a common format (such as 
particle volume vs , size or particle count vs. size) and 
combined into one distribution. The advantage is that the 
frequency range for the particle pulses is so much higher than 
the frequencies of the background drift. And so these pulses 
can be measured accurately by subtracting the local signal 
baseline on either side of each pulse. At very low scattering 
angles, the scattering signal drops off by at least the fourth 
power of particle diameter. Therefore larger particle pulses 
will stand out from the signals from many smaller particles 
which may be in the beam at any instant of time. Also the 
number concentration of larger particle will be low and 
provide for true single particle counting. 

As mentioned before in this disclosure, the particle shape can 
be determined by measuring the angular distribution of 
scattered light in multiple scattering planes, including any 
number of scattering planes. The particle shape and size is 
more accurately determined by measuring the angular scattering 
distribution in a large number of scattering planes, requiring 
many detector elements in the arrays shown in figures 33 and 
44. As the number of detector elements becomes large, the use 
of less expensive 2-dimensional detector arrays, such as CCD 
arrays, becomes more attractive to take advantage of the 
economies of scale for production of commercial CCD cameras.. 
However, the use of these 2-dimensional detector arrays 
presents some problems, which are not associated with custom 
detector arrays with optimally designed elements as shown 
previously. These arrays usually have poor dynamic range, poor 
sensitivity, poor A/D resolution, slow digitization rates, and 
high levels of crosstalk between pixels (blooming for 
example) . Methods for mitigation of these problems are 
described below. 

The detector array could be scanned at a frame rate, where 
during the period between successive frame downloads (and 
digitizations) each pixel will integrate the scattered light 
flux on its surface during an entire passage of only one 
particle through the source beam. Each pixel current is 
electronically integrated for a certain period and then its 
accumulated charge is digitized and stored; and then this 
cycle is repeated many times. During each integration period 
the pixel detector current from scattered light from any 
particle which passes through the beam will be integrated 
during the particle's total passage through the light source 
beam. Therefore the angular scattering distribution for that 
particle will be recorded over a large number of scattering 
planes by all of the detector elements in the array. This 2- 
dimensional scattering distribution could be analyzed as 
described previously, using a large number of simultaneous 
equations and more shape parameters, by assuming a certain 
model for the particle shape (ellipsoidal, rectangular, 
hexagonal, etc.). As shown before, the particle shape and 
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random orientation can be determined from these equations . 
Also, conventional image processing algorithms for shape and 
orientation can be used on the digitized scattering pattern to 
find the orientation (major and minor axes, etc.) and 
dimensions of the scatter pattern. The particle size and shape 
can be determined from these dimensions. Also the particle 
size and shape can be determined from the inverse 2- 
dimensional Fourier transform of the scattering distribution 
for particles in the Fraunhofer size range, but with a large 
computation time for each particle. The inverse Fourier 
transform of the scattering distribution, which is measured by 
the 2 dimensional detector array, will produce an image, of 
the particle, from which various dimensions can be determined 
directly. 

For example, consider an absorbing rectangular particle of 
width and length dimensions A and B, with both dimensions in 
the Fraunhofer size range and minor and major axes along the X 
and Y directions , The irradiance in the scattering pattern on 
the 2-dimensional detector array will be given by: 

I(a,b) = lo * (SINC(na) *SINC(nb) ) ^2 

Where SlNC(x) = SIN(x)/x 



Jo is the irradiance in the forward direction at zero 
scattering angle relative to the incident light beam direction 

a = A^sin (anga) /wl 
b = B^sin (angb) /wl 
where : 

wl = wavelength of the optical source 

anga - the scattering angle relative to the incident source 
beam direction in the scattering plane parallel to the A 
dimension of the particle 

angb = the scattering angle relative to the incident source 
beam direction in the scattering plane parallel to the B 
dimension of the particle 

The corresponding x and y coordinates on the 2 dimensional 

detector array will be: 

X = F*tan(anga) and y = F*tan(angb) 

The scattering pattern crossections in the major and minor 
axes consist of two SINC functions with first zeros located 
at: 

xo = F*tan (arcsin (wl/A) ) 
yo = F-^tan (arcsin (wl/B) ) 

where F is the focal length of the lens 3 in Figure 32 or F = 
M*F3 in figure 49 where F3 is the focal length of lens 3 in 
figure 49 and M is the magnification of lens 4 between the 



71 



wo 2005/091970 



PCT/US2005/007308 



jDacK tocal plane of lens 3 (source block) and the detector 
array, which is in an image plan4 of said back focal plane • 
By inspection of these equations, the dimension of the 
scattering distribution is inversely proportional to the 
particle dimension along the direction parallel to direction 
of the dimension measurement. So for a rectangular particle, 
use known image processing methods to determine the major and 
minor axes of the scattering pattern and then measure the 
width and length of the pattern at the first zeros (xo and yo) 
in the scattering distribution in the directions of the major 
and minor axes. Then the particle dimensions are given by: 

A = wl/ (sin (arctan (xo/F) ) 
B = wl/ (sin (arctan (yo/F) ) 

These equations describe the process for determining particle 
shape for a randomly oriented rectangular particle where we 
have assumed that the particle is much smaller than the 
uniform intensity portion of the source beam. Other parameters 

(such as the point in the scatter distribution which is 50% 
down from the peak) which describe the width and length of the 
scattering pattern can be used instead of xo and yo, but with 
different equations for A and B. In general, the corresponding 
particle dimensions can be determined from these parameters, 
using appropriate scattering models which describe the 
scattering pattern based upon the effects of particle size, 
shape, particle composition and the fact that the scattering 
pattern was integrated while the particle passed through a 
light source spot of varying intensity and phase. This 
analysis for rectangular particles is one example for 
rectangular particles. The model for each particle shape 

(polygon, ellipsoid, cube, etc.) must be computed from 
scattering theory for nonspherical particles using algorithms 
such as T-matrix method . 

The hardware concept is shown in figure 49. This i^ystem is 
very similar to those shown previously in this disclosure in 
Figure 32. Pinhole 1 removes high angle background from the 
light source and lens 2 collimates the source for passage 
through the sample cell through which the particle dispersion 
flows. The light source, lenses 1 and 2, and pinhole 1 could 
be replaced by a nearly collimated light beam such as a laser 
beam. Two optical systems view the particles. The 2- 
dimensional array #1 measures the scattering distribution from 
each particle and 2-dimensional array #2 measures the image of 
each particle. Array #1 is used to measure the dimensions of a 
smaller particle and array #2 measures the larger particles 
where the array pixel size can provide sufficient size 
resolution as a percentage of particle dimension for accurate 
dimension measurement. The scattering pattern from the 
particle is formed in the back focal plane of lens 3 and this 
scatter pattern is imaged onto array #1 by lens 4. A small 
block is placed in the back focal plane of lens 3 to block the 
unscattered focused light from the light source so that it 
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Will not reach array #1. The source light would saturate some 
pixels on that array and these pixels may bloom or crosstalk 
into adjacent pixels where very low level scattered light is 
being measured. The source block could also be replaced by an 
annular spatial mask (as shown previously) to measure scatter 
only over a certain range of scattering angles. Aperture 1 
(also see figure 50) is placed in an image plane of the sample 
cell to define a restricted region of the beam where particles 
will be counted. This region is confined to where the 
intensity profile of the source beam has sufficient 
uniformity. If this region of confinement is not required and 
access to the surface of the CCD is available (windowless CCD 
array) then lens 4 and aperture 1 could be removed and the CCD 
array could be placed in the back focal plane of lens 3, 
directly behind the source block. The beam diameter and lens 
diameters are not drawn to scale. They are drawn to display 
the details of beam divergence and conjugate planes. Lenses 
2B, 3 and 3B might actually be much larger than the beam 
diameter to collect scatter over a large range of scattering 
angles. 

A second similar optical system (system B which contains 
lenses IB, 2B, 3B, 4B etc.) is placed upstream of the particle 
flow from the system (the main system which contains lenses 1, 
2, 3, 4, etc.) described above (see figure 49). This system B 
reduces many of the problems associated with CCD arrays which 
are mentioned above. System B measures the scattered light 
from each particle before it passes through the main system 
described above. This scattered light level determines the 
particle size and predicts the signal levels which will be 
seen from that particle .when it passes though the main system. 
These predicted levels provide the ability for the system to 
either adjust the intensity of light source 1 or the gain of 
array #1 to nearly fill the range of the analog to digital 
converter which digitizes the scattering pattern data from 
array #1 and to improve signal to noise. The analog to digital 
converter and array pixel dynamic signal range is not 
sufficient to measure scattered light levels from particles 
over a large range of sizes. For example, the lowest scatter 
angle signal will change over 8 orders of magnitude for 
particles between 1 and 100 microns. However, the dynamic 
range of most CCD arrays is between 200 and 1000. Therefore, 
by adjusting the source intensity so that the maximum pixel 
value on array #1 will be just below saturation for each 
particle, the optimum signal to noise will be obtained. The 
time of the pulse from the upper system B will predict when 
each particle will pass through the main system, using the 
flow, velocity of the dispersion through the cell. So the array 
only needs to integrate during the particle's passage through 
the source beam. This minimizes the . integration time and shot 
noise of array #1. This timing could also be used to pulse the 
laser when the particle is in the center of the source beam 
for imaging by array #2 to freeze the particle motion during 
the exposure. Also the predicted size from system B could be 
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used to choose only particles in a selected size range for 
shape measurement. Or some smaller particles could be passed 
without dimensional measurement to increase the statistical 
count of larger particles relative to. the smaller particles to 
improve the counting statistics for the larger particles which 
are usually at lower number concentration than the smaller 
particles. But the size distribution could then be corrected 
by the total count distribution from the upper system B, while 
the particle shape count distribution is determined from the 
fewer particles counted by the main system. The size of the 
scattering pattern could also be predicted by system B so that 
an appropriate sub-array of array #1 would be digitized and 
analyzed to save digitization time. The size prediction can 
also determine which array (#1 , #2 or both) will be digitized 
to determine the particle dimensions. The upstream system B 
particle counts could also be used to determine coincidence 
counting while the, particle concentration is being adjusted. 
Many of these techniques are used to reduce the digitization 
load on the analog to digital converter, if required. 

Lens 4B acts as a field lens to collect scattered light and 
place the scatter detector in the image plane of the sample 
cell. This detector could be a single element detector which 
simply measures the all of the scattered light over a large 
range of scattering angles. However, this single detector 
measurement could be complicated by the variation of light 
intensity across the source beam. The use of a three detector 
multi-element detector (see Figure 50) could be used in this 
image plane of the sample cell. Then only particles which 
produce signal primarily on the center detector (of the three 
detector set) would be accepted for counting. This particle 
selection could be based upon the ratios between the signal 
from the center detector element with the corresponding signal 
from either of the outer elements, as described previously in 
Figure 37. The particle will be counted only if these two 
ratios are both above some threshold. If a large single 
detector is used instead of the multi-element detector, lens 
4B could be removed and the detector could be placed directly 
behind the source block if it is large enough to collect all 
of the scattered flux. Otherwise a lens should be used to 
collect the scattered light and focus it onto the detector. 

If the upstream system B is not used, the CCD array scans of 
each scatter pattern should be made over multiple long periods 
(many individual particles counted per period with one array 
scan per particle) where the light source intensity or 
detector pixel gain is chosen to be different during each 
period. In this way particles in different size and scattering 
efficiency ranges will be counted at the appropriate source 
irradiance or detector pixel gain to provide optimal signal to 
noise. So during each period, some particles may saturate the 
detectors and other particles may not be measured due to low 
scattering signals. Only particles whose scattering efficiency 
can, produce signals within the dynamic range of the array for 
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that chosen light source level or gain will be measured during 
that period. So by using a different source level or gain 
during each period, different size ranges are measured 
separately, but with optimal signal to noise for each size 
range. The counting distributions from each period are then 
combined to create the entire size and shape distribution. 
This method will require longer total measurement time to 
accumulate sufficient particle counts to obtain good accuracy 
because some particles will be passed without counting. The 
use of system B to predict the optimal light source level or 
pixel gain provides the optimum result and highest counts per 
second. 

The main system counting capability, as shown in figure 4 9, 
could be added to any diffraction ensemble system by using the 
scatter collection lens (lens 3 in Figure 53) in the ensemble 
system to act as lens 3 in the counting system in figure 49. 
The light path after the ensemble system scatter collection 
lens (the lens forming the scatter pattern) would be partially 
diverted, by a beamsplitter, to a detection system as shown in 
figure 49 after lens 3. This detection system could be any 
appropriate variation of the detection system (with or without 
array#2 and its beamsplitter as shown in figure 49 for 
example) , Also if source region confinement (as discussed 
above) is not required and access to the surface of the CCD is 
available (windowless CCD array) then lens 4 and aperture 1 
could be removed and the CCD array could be placed in the back 
focal plane of lens 3, behind the source block. System B could 
also be added to the ensemble scattering system, but with 
significant added expense. 

Linear CCD arrays do not have sufficient dynamic spatial range 
to accurately measure scatter pattern profiles from particles 
oyer a large range of particle size. For example, for a 
million pixel array, the dimensions are 1000 by 1000 pixels. 
If at least 10 pixel values are needed to be measured across 
the scatter profile to determine the dimension in each 
direction, then 1000 pixels will only cover 2 orders of 
magnitude in size. This size range can be increased to 4 
orders of magnitude by using two arrays with different angular 
scales. Figure 51 shows a system similar to that shown in 
figure 49, but with an additional scatter detection array (2- 
dimensional array #2) . Arrays #1 and #2 are in the back focal 
planes of lens 4 and lens 5, respectively. Lens 5 has a much 
longer focal length than lens 4, so that each pixel in array 
#1 covers a proportionately larger scattering angle interval. 
As each particle passes through system B, the particle size is 
estimated to determine which array (#1 of #2) should be 
scanned and digitized. Array #1 should be used for small 
particles which scatter over large angles and Array #2 should 
be used for larger particles. The use of two arrays with 
different angular scales provides much higher particle count 
rates. For example, for 2 orders of size magnitude with a 
single array, 1 million pixels must be digitized (1000 by 1000 
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with minimum of 10 pixels for the largest particle) . However, 
if two smaller 100 by 100 pixel arrays were used for array #1 
and array #2 and the focal length for lens #5 was 10 times 
longer than the focal length of lens #5. Then these two 10,000 
pixel arrays could cover 2 orders of magnitude in size, 
equivalent to that of the single 1 million pixel array; but 
only a maximum of 10,000 pixels must be digitized for each 
particle by using the size estimate from system B to determine 
which array to digitize. This design provides a factor of 100 
increase in the particle count rate. This rate could be 
further increased by only digitizing the minimum subarray 
needed to measure each particle, based upon the size 
prediction provided by system B. Also, figure 51 shows the use 
of separate apertures (aperture #1 and aperture #2) which have 
different size openings. A smaller opening is used for the 
smaller particle detector array to reduce the scatter voliime 
and reduce the probability of coincidence counting. 

Figure 49B shows an analog version of the laser power control 
by system B. The peak detector receives the total signal, 
through port E, from the detector elements of the multi-- 
element detector behind lens B. This detector could also be a 
single element detector if particle position detection is not 
used to define a small interaction voliome, as described 
previously. When the peak detector breaks a threshold, it 
starts (by port C) the integration of the detector arrays in 
the main system after an appropriate delay, accounting for the 
distance between the systems and the flow velocity. When the 
integration is finished, the array (through port D) resets the 
peak detector to start to look for the next particle. The peak 
value held by the peak detector is input (input B) to an 
analog ratiometer with an adjustable reference voltage input 
A, which can be set to adjust the laser power, and hence the 
scatter signal, to nearly fill the analog to digital converter 
of the detector array in the main system. In this way, the 
light source intensity is rapidly changed to always nearly 
fill the range of the A/D converter for particles over a large 
range of size and scattered light signal. The value A/B could 
also be used to set the gain on the detector arrays, but this 
will probably not have sufficient speed. This entire process 
could also be replaced by its digital equivalents, but with 
much slower response and lower count rate. 

One note must be made about diagrams in this disclosure. The 
size of the scatter collection lens, (i.e. lens 3 in Figure 49 
and 50) is not shown in proper size relationship to the source 
beam in order to show more detail of the source beam and 
different focal planes in the design. This is true for all 
scatter collection lenses shown in this disclosure. In all 
cases we assume that the scatter collection lens is of 
sufficient diameter to collect scattered light from the 
particles over all of the scattering angles being measured. In 
some cases this may require the lens diameter to be much 
larger than the diameter of the source beam. 
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This disclosure also describes concepts for combining three 
different particle size measurement modalities: particle 
counting, ensemble scattering measurements, and dynamic light 
scattering. In this case, particle counting is used for the 
largest particles (>100 microns) which have the largest 
scattering signals and lowest particle concentration and least 
coincidence counts. The angular scatter distribution from a 
particle ensemble is used to determine particle size in the 
mid-sized range (0.5 to 100 microns). And dynamic light 
scattering is used to measure particles below 0.5 micron 
diameter. These defined size range break points, 0.5 and 100 
microns, are approximate. These methods will work' over a large 
range of particle size break points because the useful size 
ranges of these three techniques have substantial overlap: 

Single beam large particle counting (depends on the source 
beam size) 10 to 3000 microns 
Particle ensemble 0.1 to 1000 microns 
Dynamic light scattering 0.001 to 2 microns 

One problem that can be solved by particle counting is the 
problem of background drift in ensemble scattering systems 
which measure large particles at low scattering angles. An 
ensemble scattering system measures the angular distribution 
of scattered light from a group of particles instead of a 
single particle at one time. Figure 53 shows an ensemble 
scattering system (except for detectors Bl and B2 which 
illustrate additional detectors) which illuminates the 
particles with a nearly collimated light beam and collects 
light scattered from many particles in the dispersion which 
flow through the sample cell. The light source is focused 
through pinhole 1, which removes high angle defects in the 
beam intensity profile. Lens 2 collimates the beam through the 
sample cell. Lens 3 collects scattered light from the 
particles in the sample cell and focuses that light onto a 
detector array in the back focal plane of lens 3 . An example 
of the detector array design is shown in Figure 54. The 
optical system measures scattered light in certain angular 
ranges which are defined by the set of detector elements. The 
elements can have different shapes, but in general the 
scattering angle range for each element is determined by the 
radius from the optical axis in the back focal plane of lens 
3. In some cases, the detector array will have a central 
detector, DO, which captures the light from the source beam. 
Detectors Dl, D2, etc. collect various angular ranges of 
scattered light. Each detector element is connected to its own 
separate electronic integrator, which is connected to a 
multiplexing circuit and analog to digital converter (ADC) as 
shown in figure 54 for detectors D3, D4, D5, and D6 . This 
multiplexer sequentially samples each of the integrators which 
may integrate while many particles pass through the beam. So 
particle pulses cannot be measured in the ensemble system. All 
detector elements are connected to the multiplexer through 
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mregrators in a particle ensemble measuring system. Figure 54 
shows the modification, for detector elements DO, Dl and D2 
proposed by this invention. ' 

The detector elements which measure the low angle scatter (for 
example Dl and D2 ) usually see a very large scattering 
background without particles in the sample cell This 
background is due to debris on optical surfaces "or poor laser 
beam quality. Mechanical drift of the optics can cause this 
background light to vary with time. Usually the detector array 
xs scanned with only clean dispersant in the sample cell to 
produce background scatter readings which are then subtracted 
from the subsequent readings of the actual particle 
dispersion. So first the detector integrators are integrated 
and scanned without any particles in the sample cell and then 
particles are added to the dispersion and the detector 
integrators are integrated and scanned a second time. The 
first background scan data is subtracted from this second scan 
for each detector element in the array. However, if the actual 
background drifts between the two scans, a true particle 
scattering distribution will not be produced by the difference 
between these two scans. 

A much better solution is to connect each of the detector 
elements, for the lowest angle scatter, to individual analog 
to digital converters, or peak detectors as disclosed before 
by this inventor. Then these signals could be analyzed by 
many of the counting methods which were disclosed by this 
inventor. This would essentially produce an ensemble/counting 
hybrid instrument which would produce counting distributions 
for the large particles at low scattering angles and 
deconvolved particle size distributions from the long time 
. integrated detector elements (ensemble measurement) at higher 
scattering angles for the smaller particles. These 
distributions can be converted to a common format (such as 
particle volume vs. size or particle count vs. size) and 
combined into one distribution. The advantage is that the 
■frequency range for the particle pulses is much higher than 
the frequencies of the background drift. And so these pulses 
can be measured accurately by subtracting the local signal 
baseline (under the pulse), determined from interpolation of 
the signals on the leading and trailing edge of each pulse, 
using the digitized signal samples. At very low scattering 
angles, the scattering signal drops off by at least the fourth 
power of particle diameter. Therefore larger particle pulses 
will stand out from the signals from many smaller particles 
which may be in the beam at any instant of time. Also the 
number concentration of larger particle will be low and 
provide for true single particle counting. 

The smallest particles are measured using dynamic light 
scattering as shown in Figure 55. A fiber optic dynamic light 
scattering system, as described previously by the inventor, is 
inserted into the tubing through which the particle dispersion 
flows. The counting and ensemble scattering measurements are 
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made with dispersion flowing through the system. This flow 
would be turned off during th6 collection of dynamic light 
scattering signals to avoid Doppler shifts in the scattering 
spectriom due to particle motion. 

The particle counting uses the lowest angle zones (Dl^ D2, 
etc.) and the beam measuring zone (DO) of the detector array 
(an example of a detector array is shown in figure 54) • Each 
of these detector elements are connected to a separate ADC to 
measure the scattering pulse in Dl, D2, etc. and the signal 
drop on detector DO as each particle passes through the 
interaction volume where the be am illuminates and from which 
the scattering detectors can receive scattered light from the 
particles. One problem is that the amount of scattered light 
is nearly proportional to the illumination intensity of the 
source on the particle- Therefore as particles pass through 
different regions in the beam they may produce different pulse 
heights. Figure 56 shows a Gaussian intensity distribution 
which might be characteristic of the cross-section of a laser 
beam. Since the probability of particle passing through this 
beam at any position is approximately the same, we can 
generate the count vs . pulse amplitude distribution in Figure 
57, which shows the count distribution for large number of 
identical particles (also shown in Figures 10 and 10b) . Notice 
that most of the particles pass through the low intensity 
portions of the intensity distribution and many particles also 
pass through the uniform intensity region which is close to 
the peak of the intensity distribution, with a region of 
lower count level in between. This broad count response to a 
group of mono-sized particles will prevent accurate 
determination of complicated particle size distributions, 
because the pulse heights may be ambiguous for various sized 
particles. For example, a large particle passing through the 
lower intensity region can produce a pulse which is very 
similar to that from a smaller particle passing through the 
higher intensity region. The region of the intensity 
distribution which can produce scattered light into the 
detectors must be truncated by apertures in the source optics 
(aperture 1 in figure 52) or in the detection optics (aperture 
2 in Figure 52) . Either of these apertures can create a 
^'region passed by aperture''' as indicated in Figure 56 and 57. 
By using either or both apertures, only the upper region of 
the count vs. pulse amplitude distribution will be seen for 
many particles of a single particle size. This truncation by 
aperture can be used in any of the systems described in this 
document to reduce the broadening of the particle count peaks 
due to intensity variations of the source. Any residual 
broadening is then removed by algorithms such as 
deconvolution . 

Another method for eliminating this intensity distribution 
effect is to use ratios of detector signals. This works 
particularly well when many of the detectors have scatter 
signals. However, for very large particles, only scattering 
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detector Dl will see a high scatter signal with high signal to 
noise- So for very large particles, the apertures described 
previously may be required to use the absolute scatter from 
Dl . Another solution is to use the ratio of the drop in DO 
(signal SO) and the increase in Dl (signal SI) due to scatter 
as shown in figure 58A. As a particle passes through the beam 
DO will decrease by approximately the total scattered light 
and Dl will increase by only the amount of light scattered 
into the angular range defined by that detector. The drop in 
DO can be determined by subtracting the minimum of drop in SO 
from the baseline AO to produce a positive pulse AO-SO as 
shown in figure 58B. As shown in Figure 59, the ratio of 
either the integral or the peak value of the corresponding 
pulses from these two signals can be used to determine the 
size of the counted particle for the largest size particles 
which have insufficient scatter signal in D2 to. produce a 
ratio between SI and S2 . As long as D2 has sufficient scatter 
' signal and D2 captures a portion of the primary lobe of the 
angular scatter distribution, the ratio between SI and S2 will 
produce more accurate indication of size, than a ratio between 
SO and SI. The primary lobe of the scattering distribution is 
the portion of the distribution from zero scattering angle up 
to the scattering angle where the size information becomes 
more ambiguous and particle composition dependent. Usually 
this happens when the scatter function first drops below 20% 
of the zero angle (maximum) value of the function. For a 
certain range of smaller particles, the ratio between SI and 
S3 (if D3 were connected' to an A/D as D2) may have higher 
sensitivity to particle size than the ratio of SI to S2. For 
smaller particle diameters, ratios to larger angle scatter 
signals will provide better sensitivity- AO-SO, SI, etc. could 
be also be analyzed using the other methods described in this 
document. 

The signal ratio technique is needed when the ^'^region passed 
by aperture" in figures 5 6 and 57 is too large such that mono- 
sized particles produce pulse peaks over a large amplitude 
range. For example, if no aperture were used, then mono-sized^ 
particles will produce the entire count distribution shown in 
Figure 57, with ambiguity between small particles passing 
through the center of the Gaussian intensity distribution and 
large particles passing through the tail of the distribution. 
In cases where the ^^region passed by aperture'' is too large, 
the use of signal ratios (as described previously) is required 
to reduce the effect of the intensity variation (because the 
intensity variation drops out of the ratio, approximately) . If 
the source intensity distribution can be made more uniform by 
use of an aperture (aperture 1 of Figure 52) or by use of a 
non-coherent source, or if the viewing aperture (aperture 2 in 
Figure 52) of the detector only views a restricted region 
where the source intensity is more uniform, then scattering 
amplitude can be used directly to determine size as shown in 
figures 60 and 61- This may have some advantages when only one 
detector has sufficient signal so that two signals are not 
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avaixaoie to create a ratio. Also the absolute signal 
amplitude information, which is lost in the ratio calculation, 
can be useful in determining the particle composition and in 
eliminating pulses which are due to noise, as will be 
described in Figure 61. Figure 57 shows a count vs. pulse 
amplitude response with a ^"region passed by aperture''. This 
count distribution in the ^'region passed by aperture" is 
plotted on a logarithmic scale of S (or pulse peak or 
integral) for two different particle sizes, in Figure 60. Each 
function has an upper and lower limit in log(S). Notice that, 
in logarithmic S space, the two functions are shift invariant 
to particle diameter. The upper limit is due to particles 
which pass through the peak of the source intensity 
distribution and the lower limit is from the edge of the 
truncated source intensity profile. So that the count vs. 
log(S), Ns(log(S)), distribution from particles of the count 
vs. particle diameter distribution Nd(d) is the convolution 
between the shift invariant function in Figure 60, H{log{S), 
and the count vs. particle size distribution, Nd(d): 

Ns(log(S)) = Nd(d) 9 Hdog(S)) 

This equation is easily inverted by using iterative 
deconvolution to determine Nd(d) by using H{log(S)) to 
deconvolve Ns(log(S)). In some cases, for example when S=AO- 

50, the form of this equation may not be a convolution and a 
more generalized matrix equation must be solved. 

Ns{log(S)) =H(log(S)) ^Nd(d) 

Where H is the matrix and Nd(d) is a vector of the actual 
counts per unit size interval. Each column of matrix H is the 
measured count vs. log(S) response to a particle of size 
corresponding to the element in Nd(d) which multiplies times 
it in the matrix multiply . This matrix equation can be 
solved for Nd(d), given Ns(log(S)) and H(log{S)). This 
equation will also hold for the case where the functions of 
log{S) are replaced by other functions of S . 

Figure 61 shows a scatter plot of the counted data points in 
the two dimensional space, where the two dimensions are the 
logarithm of pulse amplitude or pulse integral for two 
different signals A or B. For example SA and SB could be 31 
and 82, or AO-SO and 

51. Two squares are shown which encompass the approximate 
region where counts from particles could occur for each of two 
particle diameters. SAU and SAL refer to the upper and lower 
limits of the signal, respectively, as shown in Figure 60. The 
lower limit, SAL, is determined by the cutoff of the aperture 
on the source intensity profile. The upper limit, SAU, is the 
maximum signal value when the particle passes through the peak 
of the source intensity profile. In the two dimensional space, 
shown in Figure 61, points are shown where the particle passes 
through the intensity peak, [log (SBU) , log (SAU) ] , and where the 
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particle passes through the edge of the intensity profile, 
[log (SBL) , log (SAL) ] , where the intensity is lowest. As 
particle size changes, this square region will move along a 
curve which describes the scatter for particles of a certain 
composition, as shown in Figure 62. The moving square will 
define a region, between the two blue lines which pass through 
the edges of the square. Real particles can only produce 
points within this region. Points outside this region can be 
rejected as noise points or artifact signals. This two 
dimensional count profile can be deconvolved using 2- 
dimensional image deconvolution techniques (as described 
previously) , because each square defines the outline of the 
two dimensional impulse response in Log(S) space. The two 
dimensional count profile is the concentration (counted points 
per unit area of log(SA) and log (SB) plot 2-dimensional space) 
of counted points at each coordinate in the log(SA) and 
log (SB) space. This count concentration could be plotted as 
the Z dimension of a 3-dimensional plot where the X and Y 
dimensions are log(SA) and log(SB), respectively. The two 
dimensional impulse function plotted in the Z dimension of 
this 3-dimensional space is determined from the product of 
functions as shown in Figure 60, one along each of the log{SA) 
and log (SB) axes of Figures 61 and 62. The 2-dimensional plot 
of Log (SI) vs. S1/S2 could also provide the particle size 
solution, as shown in Figure 39. All of the techniques for 
multi-dimensional analysis, described previously, apply to 
this case. In this case, the 2 dimensional space would 
deconvolved in the Log (SI) direction only, where the function 
broadening due to source intensity variation occurs. Very 
little broadening will occur in the S1/S2 function due to the 
ratio correction. These multi-dimensional views and the 
previously described methods apply to all combinations of 
signals (SI, S2 ,etc.), Log(signals) (Log(Sl), Log(S2), etc). 
Ratios of signals (S2/S1, S3/S1, etc), and/or any combinations 
of these (S1/S2 and Log(Sl), etc.). Where the signals SI, S2, 
etc- are functions of the measured scatter signal (such as 
peak value, integral, value at a certain time during the 
pulse, product of two scatter signals, etc.) 

The previous concept for ensemble particle systems uses 
particle counting to eliminate the particle size errors caused 
by background drift in the angular scattering signals, because 
the frequency content of the counted pulses is much higher 
than the background drift, and so the pulses can be detected 
by methods described previously by this inventor, without 
being effected by background drift. The local baseline is 
easily subtracted from each pulse because the background drift 
is negligible during the period of the pulse. However, this 
advantage can also be used with the integrators as shown in 
Figure 63 . The slowly varying baseline can be removed by high 
pass analog electronic filters with a cutoff frequency between 
the lowest frequency of the particle scatter pulse spectrum 
and the highest frequency of the background drift spectrum. 
The input to each of the integrators which follow each high 
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pass filter are the particle scatter pulses, without 
background which is attenuated by the filter. These pulses can 
be integrated and multiplexed into the same analog to digital 
converter as the higher scattering angle signals, which do not 
need the highpass filtering to remove the baseline drift. 
These integrators integrate over an extended period where many 
particles pass through the beam. In the case of smaller 
particles, there may be many particles in the beam at any 
instant in time. However, since the scatter signal from larger 
particles is much larger than that for smaller particles and 
with many smaller particles in the beam, these smaller 
particle signals will have very low fluctuations relative to 
the discrete pulses from the larger particles. So this high 
pass filtering will select the larger particles where the 
scatter signal fluctuations are large. This measurement could 
also be made with an RMS (root mean squared) module which only 
detects the higher frequency portion of the scatter signal for 
the lower angle detectors. All of these integrated signals, 
from low and high angle scatter detectors, are then inverted 
by techniques such as deconvolution . The detector signals 
could also be digitized directly; and the filtering and 
integration steps could be done digitally. However, the 
optical scatter model for the deconvolution must include the 
loss of the small particle contribution to these filtered 
signals, because as the number of particles in the beam 
increases, the higher frequency components will be attenuated 
due to overlap of pulses. Given this attenuation process and 
the fact that signals at the lowest scattering angles scale as 
the fourth power of particle' diameter, the smaller particle 
signals should not be significant in these filtered low angle 
signals. The correction to the model is very minor; 
essentially the small particle contribution to these filtered 
detector signals can be assumed to be small in the scattering 
model . 

These methods do not assume any particular number of lower 
angle zones. For example, DO, Dl, D2, and D3 could be handled 
with the techniques above. Essentially, any detectors with 
background drift problems should be handled with these 
methods . 

Normally all of the ADC scans of the multiplexer output are 
summed together and this sum is then inverted to produce the 
particle size distribution. But due to the large difference in 
scattering efficiency between large and small particles, 
smaller particles can be lost in the scatter signal of larger 
ones in this sum. This problem can be mitigated by shortening 
the integration time for each multiplexer scan and ADC cycle 
to be shorter than the period between pulses from the large 
particles. Then each multiplexer scan and subsequent 
digitization can be stored in memory and compared to each 
other for scattering angle distribution. ADC scans of similar 
angular distribution shape are summed together and inverted 
separately to produce multiple particle size distributions. 
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Then these resulting particle size distributions are summed 
together, each weighted by the amount of total integration 
time of its summed ADC scans. In this way, scans which contain 
larger particles will be summed together and inverted to 
produce the large particle size portion of the size 
distribution and scans which contain only smaller particles 
will be summed together and inverted to produce the small 
particle size portion of the size distribution, without errors 
caused by the presence of higher angle scatter from larger 
particles - 

Another method to measure larger particles is to place a 
sinusoidal target in an image plane of the sample cell on 
front of a scatter detector as described previously by this 
inventor. The dispersant flow could be turned off and particle 
settling velocity measured by the modulation frequency of the 
scatter signal from individual particles settling through the 
source beam. The hydrodynamic diameter of each particle can 
then be determined from the particle density, and dispersant 
density and viscosity. 

Finally the three size distributions from dynamic light 
scattering, ensemble scattering and counting are combined to 
produce one single distribution over entire size range of the 
instrument by scaling each size distribution to the adjacent 
distribution, using overlapping portions of the distribution. 
Then segments of each distribution can be concatenated 
together to produce the complete size distribution, with 
blending between adjacent distributions in a portion of each 
overlap region. This method works well but it does not make 
most effective use of the information contained in the data 
from the three sizing methods. Each inversion process for each 
of the three techniques would benefit from size information 
produced by other techniques which produce size information in 
its size range. This problem may be better solved by inverting 
all three data sets together so that each of the three methods 
can benefit from information generated by the others at each 
step during the iterative inversion process. For example, the 
logarithmic power spectrum (dynamic light scattering) , 
logarithmic angular scattering distribution and logarithmic 
count distribution could be concatenated into a single data 
vector and deconvolved using an impulse response of likewise 
concatenated theoretical data. However, in order to produce a 
single shift invariant function, the scale of the counting 
data must be changed to produce a scale which is linear with 
particle size. For example, the pulse heights on an angular 
detector array will scale nearly as a power function of 
particle size, but the power spectrum and ensemble angular 
scattering distributions shift along the log frequency and log 
angle axes linearly with particle size. So a function of the 
pulse heights must be used from the count data to provide a 
count function which shifts by the same amount (linear with 
particle size) as the dynamic light scattering and ensemble 
distributions . This function may vary depending upon the 
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particle size range, but for low scattering angles the pulse 
height would scale as the fourth power of the particle 
diameter, so that the log of the quarter power of the pulse 
heights should be concatenated into the data vector. This 
technique will work even though the concatenated vectors are 
measured verses different parameters (logarithm of frequency 
for dynamic light scattering, logarithm of scattering angle 
for ensemble scattering, and logarithm of pulse height or 
integral for counting), simply because each function will 
shift by the same amount, in its own space, with change in 
particle diameter. And so the concatenation of the three 
vectors will produce a single shift invariant function which 
can be inverted by powerful deconvolution techniques to 
determine the particle size distribution. This technique can 
also be used with any two of the measurement methods (for 
example: ensemble scattering and dynamic light scattering) to 
provide particle size over smaller size ranges than the three 
measurement process. This problem can also be formulated as a 
matrix equation, where the function variables can be Log(x) or 
X (where x is the variable frequency (dynamic light 
scattering) , scattering angle (ensemble angular scattering) or 
S (the counting parameter) ) . Again these functions can be 
concatenated into vectors and a matrix of theoretical 
concatenated vectors. And this single matrix equation, which 
contains the dynamic light scatter, the ensemble scatter and 
the count data, can be solved for the differential particle 
volume vs. size distribution, Vd. 

Fm = Ht*Vd 

Where Fm is the vector of measured values which consist of 
three concatenated data sets (dynamic, angular, and counting). 
Ht is the theoretical matrix, whose columns are the 
theoretical vectors which each represent the theoretical Fm of 
the size corresponding to the value Vd which multiplies that 
column. This matrix equation can be solved for Vd, given Fm 
and -Ht-. _ 

If the convolution form holds, then the equation becomes: 

Fm = Him e Vd 

Where Him is the Fm response at a single particle size and 9 
is the convolution operator. This equation can be solved for 
Vd, given Fm and Him. 

Another way to accomplish this is to constrain the inversion 
process for each technique (dynamic light scattering, ensemble 
scattering and counting) , to agree with size distribution 
results from the other two techniques in size regions where 
those other techniques are more accurate. This can be 
accomplished by concatenating the constrained portion of the 
distribution, Vc, onto the portion (Vk) which is being solved 
for by the inversion process during each iteration of the 
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inversion. The concatenated portion is scaled relative to the 
solved portion (AVc) , at each iteration, by a parameter A 
which is also solved for in the inversion process during the 
previous iteration. This can be done with different types of 
inversion methods (global search, Newton's method, Levenburg- 
Marquart, etc.) where the scaling parameter A is solved for as 
one additional unknown, along with the unknown values of the 
particle size distribution. This technique will work for any 
processes where data is inverted and multiple techniques are 
combined to produce a single result. 

Fn = Hnm^Vn (matrix equation describing the scattering 
model) 

Vn = Vk 1 AVc (concatenation of vectors Vk and AVc, n number 

of total values in Vn) 

Solve for k values of Vk and constant A 

Vn = Fn/Hnm (solution of the matrix equation by iterative 
techniques (not a literal division) ) 

k ^ n+1 

Another hybrid combination is particle settling, ensemble 
scattering, and dynamic light scattering as shown in figure 
64. As before, dynamic light scattering probes a portion of 
the particle dispersion flow stream, with the flow turned off. 
The ensemble scattering system uses a detector* array to 
measure the angular scattering distribution from groups of 
particles in the sample cell as the dispersion flows through' 
the cell. A particle settling measurement is used for the 
largest particles which have the highest settling velocities. 
The settling is measured by sensing the power spectrum of the 
scattered light as viewed through a set of sinusoidal or 
periodic masks, which are also referred to as a multi- 
frequency modulation transfer target. Some examples of these 
masks are shown in figures 16 and 66. The mask can be placed 
between lens 2 and lens 3 or on front of a group of detectors 
as shown in figure 64. The detector is placed in the back 
focal plane of lens 3, as shown previously by the inventor, to 
collect scattered light in separate ranges of scattering 
angle. A portion of the scattered light is split off by a beam 
splitter to an aperture in the focal plane of lens 3. This 
aperture can be an annular opening, which passes a certain 
range of scattering angles (signal rises as particle passes 
through a bright fringe) , or a pinhole centered to pass only 
the focused spot of the source (signal drops as particle 
passes through a bright fringe) in the back focal plane of 
lens 3. Where the fringe is defined in the sample cell plane 
as an image of each highly transmitting line in the multi- 
frequency target. The light passing through the aperture, also 
passes through a periodic mask, as described previously, which 
is in a plane conjugate to the sample cell. This mask contains 
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multiple regions, each with a different spatial frequency for 
the periodic absorption or reflection pattern. Behind each 
region in the mask is a separate detector which collects the 
light which only passes through that region. As particles pass 
through a region, the scattered light (for the annular 
aperture) or the attenuation (for the pinhole passing the 
source) of the beam are modulated by the motion of the 
particle's image across the absorption cycles of the mask. The 
particle dispersion flow pump is turned off and the particles 
are allowed to settle through the sample cell. The frequency 
of signal modulation for any particle is proportional to its 
settling velocity, which indicates the hydrodynamic size of 
the particle, given the particle and dispersant densities and 
the dispersant viscosity. The signal can be digitized and 
analyzed on an individual particle basis to count and size 
individual particles by measuring the settling velocity of 
each particle- In this case zero crossing measurement or , 
Fourier transform of the signal segments for each particle 
could be used. In the case where many particles are in the 
beam at each instant, the power spectrum of the signal could 
be measured over an extended time. This power spectrum would 
then be inverted to produce the particle size distribution 
using a matrix model as described previously. As identical 
particles pass though different focal planes (planes 
perpendicular to the optic axis) in the sample cell, the power 
spectrum will change because the sharpness of the image of the 
mask will be reduced as the particle moves farther from the 
image plane. Also if the source beam is focused into the 
sample cell, as shown previously, then the source intensity 
and the scatter signal will drop as the particle passes 
farther from the best focus plane of the source. These effects 
can be included in the counting system model which is inverted 
to produce the particle size. The H function (or H matrix) 
described previously will contain columns which describe the 
count vs. signal frequency from a group of identical 
particles, of the size corresponding to that matrix coliamn,r 
passing through every point in the sample cell. For the 
ensemble scattering system model, the H function (or H matrix) 
will contain columns which describe the integrated scatter 
signal vs. angle from a group of identical particles, of the 
size corresponding to that matrix column, passing through 
every point in the sample cell. 

The following list describes the various options for using 
scattered light to measure size. In each case, the following 
matrix equation must be solved to determine V from measurement 
of F: 

F = H^V 

This equation can be solved by many different methods. 
However, because this equation is usually ill-conditioned, the 

use of constraints on the values of V is recommended, using 
apriori knowledge- For example, constraining the particle 
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count or particle volume vs . size distributions to be positive 
is very effective. In some cases, as shown previously by this 
inventor, changing the abscissa scale (for example from linear 
to logarithmic) of F can produce a convolution relationship 
between F and V, which can be inverted by very powerful 
deconvolution techniques. 



F = H e V 



Particle Counting 

1) Angular scatter or attenuation due to scatter: 
V = particle count per size interval vs. size 

F = count per signal amplitude interval vs. signal amplitude 
where signal amplitude is either pulse peak value ^ integral 
of the pulse, or other function of these values 

H = matrix where each column is the F function for the 
particle size corresponding to that column 

Response broadening mechanisms in the H matrix: 

source intensity variation in x and y directions where 
particles can pass 

(mitigated by aperturing of intensity distribution at an image 

plane of the sample cell or using diffractive or absorptive 
optic beam shapers and apodizers) 

source intensity variation in z direction 

(mitigated by double" pulse sensing and detector aperture' at 
image plane of sample cell ) 

Advantages: high resolution and aerosol capability 
Disadvantages: counting statistic errors for low count 



2) Settling (hydrodynamic size) 

V = particle count per size interval vs. size 

F = count per signal frequency interval vs. signal frequency 
where signal frequency is the frequency of the scatter signal 
segment for the counted particle 
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H = matrix where each column is the F function for the 
particle size corresponding to that column 

Response broadening mechanisms in the H matrix: 

Finite length of modulated signal segment from each particle 

Brownian motion 

Variation of signal frequency along z direction 

Advantages: high size resolution, excellent detection of 
small particles mixed with large particles, excellent 
measurement of low tails in the size distribution 

Disadvantages: counting statistic errors for low count; and 
possible difficulty measuring large particles in aerosols due 
to very high settling velocities 

Ensemble scattering 

1) Angular scatter, or attenuation due to scatter 

V = particle volume per size interval vs. size 

F = scattered light flux per scattering angle interval vs . 
scattering angle 

H = matrix where each column is the F. function for the 
particle size corresponding to that column 

Response broadening mechanisms in the H matrix: 

The broad angular range of scatter from a single particle 

Advantages: excellent size reproducibility 

Disadvantages: low size resolution, poor detection of small 
particles mixed with large particles, poor measurement of low 
tails in the size distribution. 

2) Settling (hydrodynamic particle size) 

V = particle volume per size interval vs. size 

F = scattered light detector current power per frequency 
interval vs . frequency 

H = matrix where each column is the F function for the 
particle size corresponding to that column 

Response broadening mechanisms in the H matrix: 

Finite length of modulated signal segment from each particle 
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Brownian motion 

Variation of signal frequency along z direction 

Advantages: high size resolution, excellent detection of 
small particles mixed with large particles, excellent 
measurement of low tails in the size distribution 

Disadvantages: difficulty measuring large particles in 
aerosols due to very high settling velocities 

Another system for counting and sizing particles, using 
imaging, is shown in figure 67. There are two optical systems, 
using light source 1 and light source IB. The source 1 system 
measures larger particles by direct imaging of the particles, 
flowing through the sample cell, onto 2~dimensional detector 
array #1. The source IB system produces a cone shaped 
illiominating beam, using the source block on lens IB, which 
defines an illuminated focus volume in the particle 
dispersion. This focal volume is imaged by lens 3B though the 
mirror and beamsplitter onto the same 2-dimensional detector 
array #1. The focal volume is placed close to the sample cell 
window which is closest to the 2-dimensional array. Before and 
after the focal point of the source IB in the cell, the 
illumination beam is a hollow cone, which provides an un- 
illuminated volume through which lens 3B can view the focal 
point of source IB in the sample cell. The 2-dimensional array 
is multiplexed between the two systems by sequentially turning 
on either source 1 or source IB. Each source is pulsed so as 
to only illuminate the flowing particles during a travel 
distance which is less than the required imaging resolution. 
Alternately, the flow can be stopped during the exposure to 
eliminate any smearing of the image. 

The source 1 system can take many forms. In Figure 67, 
aperture 1 blocks the scattered light and passes the un- 
scattered light (approximately) so that particles will appear 
as dark on a bright background. In Figure 68, aperture 1 
blocks the un-scattered light and passes the scattered light 
(approximately) so that particles will appear as bright on a 
dark background. In Figure 68, pinhole 1 is used to remove any 
higher angle components of light source 1 which could create 
background light which can pass around the light block. In 
either case, contiguous detector array pixel values, which are 
above some threshold in Figure 68 or below some threshold in 
Figure 67, can be combined to determine the total light 
extinction of the particle (in Figure 67) or the light 
scattered in the acceptance angle of lens 3 in Figure 68. 
These values can be used to determine the size of particles 
throughout the size range, even for a particle whose image 
size is less than the size of single pixel, as described 
previously in this document. However, for very small 
particles, many particle coincidence counts may occur in the 
source 1 system and the signal to noise may drop below 
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acceptable levels . So the smaller particles are measured by 
the source IB system, which focuses the source to a higher 
irradiance in the sample cell and defines a much smaller 
interaction volume than the other system. The source block on 
lens IB creates a hollow cone of light, which is focused close 
to inner sample cell window wall which is closest to the 
detector array. This is shown in more detail in Figure 69, 
Lens 3B collects scattered light from the particles, with a 
field of view which falls inside of the hollow cone of light. 
Therefore, only particles in the source focal volume will 
contribute to the scattered light and the image formed on the 
detector array. Since the focus is close to the inner wall, 
few particles will block or rescatter light which is scattered 
by particles in the focal volume. However, this method will 
produce good results for any location of the lens IB focus, as 
long as lens 3B is focused to the same- location. Also 
particles in the extended illuminated portions of the hollow 
cone cannot contribute scatter to lens 3B due to the limited 
acceptance cone of this lens. The optical magnification is 
chosen such that the 2-dimensional array only sees the focal 
volume, without seeing any light from the hollow light cones 
on the input or exit of the focal volume. The scattered light, 
accepted by lens 3B, is reflected by a mirror and a 
beamsplitter, through lens 4, to the detector array. The 
lenses 3, 3B, and 4 are designed to work at infinite 
conjugates, however lens 4 could be removed and lens 3 and 
lens , 3B could be adjusted to create images of the particles 
directly onto the detector array at finite conjugates. In both 
Figures, the size of very large particles can be measured 
directly by the size of digitized image to avoid errors in the 
magnitude of the scattered light from these large objects 
which only scatter at very small scattering angles, where the 
background light is high. 

As shown in Figure 69, the. Source IB system also uses a 
concave inner surface whose center of curvature is 
approximately coincident with the focal volume. This design 
produces a focal volume which does not shift with change in 
the refractive index of the dispersing fluid in the sample 
cell. A concave surface could also be used on the opposing 
window to control the focal shift of lens 3B due to refractive 
index change of the dispersant, to maintain sharp focus of the 
scattered light in the particle image on the detector array. 
Again, the center of curvature would be coincident with the 
focal volume. However, Figure 6 9 shows another alternative 
which may be more flexible and provide better focus precision. 
Lens 3B is attached to a focus mechanism which can move the 
lens to various focal positions. This mechanism could be any 
appropriate mechanical means, including motor or piezoelectric 
drivers- The position of lens 3B, along the optical axis, is 
changed under computer control to maximize particle edge 
sharpness in the image on the detector array. This sharpness 
could be determined by many criteria which include the spatial 
derivative of the image. In the case where the depth of field 
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of lens 3B is shorter than the depth of the focal volume of 
lens IB, this focus adjustment can be used to only select 
particles which are in sharp focus, by measuring the edge 
sharpness of each particle in the field of view, at three 
different focal positions. Only particles, whose edge 
sharpness is maximum in the middle focal position, are sized 
and counted. In this way, only particles which are accurately 
sized are counted- The maximum edge sharpness for each 
particle may vary among different particles which may have 
soft edges. So by measuring the edge sharpness in three 
different planes, the particles which are located in the 
middle plane can be selected by choosing only the ones whose 
sharpness is maximum in that plane. The edge sharpness could 
be determined by the spatial derivative of the intensity 
profile at the edge of each particle. This could also be 
calculated from the maximum of the spatial derivative of the 
entire particle, because this usually occurs at the particle's 
edge. The derivative could also be calculated from a smoothed 
version of the image, if image noise is a problem. This 
comparison can be done while the particles are stationary or 
by using 3 successive source pulses with a detector array scan 
during each pulse. 

Also, the hollow cone source in Figure 69 could be replaced by 
a single focused light beam, which is focused through the 
focal volume and projected at an angle to the optical axis of 
lens 3B, such that it is not captured by lens 3B, as shown in 
Figure 70. The scattering interaction volume is the 
intersection of the viewing focus of lens SB and the source 
focus of lens IB. The 2-dimensional detector array sees the 
image of the particles at the focus of lens IB, using only 
light scattered from the particles. In this way, the 2- 
dimensional array only sees particles in a very small 
interaction volume. All of the other focusing mechanisms and 
options mentioned previously for Figure 69, also apply for 
figure 70. 

Another problem associated with counting techniques is the 
coincidence counting error. In some cases, pulses from 
individual particles will overlap as shown in figure 71, which 
shows three pulses - and the signal which represents the sum of 
those three pulses. In most cases, these pulses all have the 
same shape, but with different pulse amplitudes. For example, 
any particle passing through a Gaussian laser beam will 
produce a pulse with a Gaussian shape. The only difference 
between different pulses from different particles is the 
amplitude of the pulse and the position of the pulse in time. 
Therefore, the sum of the pulses is simply the convolution of 
a single pulse with three delta functions, each delta function 
centered at one of the different pulse positions. The general 
equation for this sum of pulses is: 

S(t) = H(t) e S (Ai*a(t-ti)) 
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Where : 

^ = sum over variable i 
t = time 

S(t) — the total signal from the overlapping pulses 
5(t) = the delta function 

ti = the time at the center of the ith pulse 
Ai = the amplitude of the ith pulse 
© is the convolution operator 

H(t) = the function describing a single particle pulse shape 

So the original pulses can be recovered from S (t) by- 
inverting the above equation, using H(t) as the impulse 
response in a Fourier trans fojnti deconvolution or in iterative 
deconvolution algorithms. As shown by figure 11, the 
individual pulse heights and areas cannot be determined from 
the sum of the pulses S(t). However, through deconvolution 
the pulses can be separated as shown in figure 72, which shows 
S(t) and the total deconvolved signal resulting after some 
degree of deconvolution of S(t). Due to signal noise in S(t), 
S(t) cannot be deconvolved down to separate delta functions, 
the deconvolution will usually stop at some point before 
artifacts are created, leaving separated pulses of finite 
width. However, these pulse heights will be proportional to 
the actual heights of the original separated pulses. So that 
by using a single scale factor on the deconvolved signal, all 
of the individual pulse heights will agree with those of the 
original separated pulses. Figure 73 shows this scaled 
deconvolved result along with the original separated pulses to 
show that the separated pulse heights are recovered by the 
decovolution process. This technique can be applied to any 
time signals which have overlapping pulses of the same shape, 
such as found in particle counting. For most laser beams, H(t) 
will be a Gaussian. However, in some cases, where the laser 
beam has been apodized or truncated to reduce the large 
intensity variation, H(t) will take on the functional form 
describing the signal vs-. time profile of a single particle 
passing through that beam, which may be flat-topped Gaussian, 
rectangular, etc. 

Figures 7 4 through 7 7 show particle shape measuring systems 
which combine the concepts of Figures 2 9 through 31, and 
Figures 42 through 45. Figure 74 shows one system in the 
first scattering plane of multiple scattering planes. A three 
scattering plane system is shown, but as mentioned before, any 
number of scattering planes may be needed to describe the 
shape of more complicated particles . The lens 4 and lens 5 
systems use multiple detector elements to measure scatter in 
each of the scattering planes on one detector array/lens 
assembly. The lens 7 and lens 8 systems are repeated in each 
of the scattering planes to measure the high angle scattered 
light. These detection systems are aligned as shown in figure 
77, so that each scattering plane element on multi-element 
detector A and multi-element detector B measures in the same 



93 



wo 2005/091970 



PCT/US2005/007308 



scattering plane as the corresponding leris 7/lens8 aperture 
openings. The configurations of the multi^element detectors 
are shown in Figure 76. This concept combines the two concepts 
described earlier in Figures 29-31 and Figures 42-45. There is 
one lens 4/lens5 system with both lens 4 and lens 5 (through 
beamsplitter 1) centered on the optical axis of the source to 
measure low angle scatter. Each detector element, shown in 
Figure 7 6, is aligned to view the same scattering plane as the 
corresponding element on the other multi-element detector, but 
all three scattering planes are measured by the same multi- 
element detector. The three Lens7/lens8 systems, which measure 
the higher angle scatter, each have only one aperture (instead 
of 3 elements for each of lens4 or lens5) for each lens with 
shapes like those in Figure 7 6 ( IC, 2C, 3C apertures have the 
same (or similar) shape as the detector elements in multi- 
element detector A and ID, 2D, 3D have shapes similar to 
those in multi-element detector B) . As before, the ratios of 
the corresponding detectors C to D and the ratio of each 
element on the multi-element detector A to the corresponding 
element on multi-element detector B, provide the particle 
dimension parameter for the corresponding scattering plane. 
These parameters are then combined using a lookup table, 
search algorithm, or regression algorithm to solve for the 
particle shape. As before, to solve for the dimensions of a 
rectangular particle with arbitrary orientation, parameters in 
3 scattering planes must be measured. The' optimal separation 
of these planes in the plane of Figure 77 is approximately 120 
degrees, with one of the planes being parallel to the flow 
direction (because many particles will align with the flow 
direction) . The search algorithm will take an initial guess at 
the width, length, and orientation angle, and then calculate 
the three parameters for that guess and then compare those 
parameters to the measured ones to generate a change in the 
width, length, and orientation for the next guess and then go 
through the same loop again. As this loop is repeated, the 
change in the width, length, and orientation diminishes as the 
algorithm approaches -the - true width, length, and orientation 
of the particle. Algorithms such as Newton's method, global 
optimization, or Levenburg Marquardt could be used. The three 
dimension parameters could also be replaced by the full set of 
12 detector signal values (6 for detectors A and B and 6 for 
detector sets C and D) for input to these search or 
optimization algorithms. This would be the case for 3 
scattering planes. For particles with more complicated shapes, 
measurements in more scattering planes would be needed to 
solve for the shape parameters and the arbitrary orientation, 
but the same methods would be used to search for the solution. 
The corresponding elements in detector C and D and in detector 
A and B could also be detector segments which view different 
ranges of scattering angle instead of different angular 
weightings (as shown in Figure 7 6) of the same range of 
scattering angles. However, the use of the use of different 
angular weightings may provide larger size range because when 
the particles become very large, very little light may fall on 
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the higher angle detectors and the ratio of high to low angl-e 
signals will become multi-valued. The ratio of 2 detectors 
with different angular weightings (figure 7 6) will have a 
smooth monotonic size dependence over a large size range. 

Figure 75 shows the use of spherical window shape on the 
sample cell to avoid focal shift of the focused source spot 
and of the focal viewing spot of the collection optics as the 
refractive index of the dispersing fluid is changed. The 
center of curvature for each surface on each spherical window 
is at the beam focal position in the sample cell. The bottom 
portion of Figure 7 5 shows a spherical cell with an inlet 
tube, which ends just above the focal spot of the source beam. 
This cell is placed into a flow loop as shown in Figure 13, 
where the pump pulls the dispersion from the outlet and 
returns dispersion to the inlet. In this way, homogeneous 
dispersion passes through the source beam directly from the 
inlet. Regions far from the inlet, in the spherical cell, may 
have inhomogeneous particle dispersion, which may not be 
representative of the entire particle sample. Also when the 
flow loop is drained, this orientation of the spherical cell 
will drain completely without leaving residual particles to 
contaminate the next particle sample. 

Another concept for measuring the shape and size of small 
particles is shown in Figure 78. This system consists of two 
scatter collection subsystems: the first subsystem using lens 
3 and detector array 1 and the second subsystem using a 
segment of a nonspherical mirror and detector array 2. 
Detector array 1 measures scattered light at low scattering 
angles and detector 2 measures light scattered at high 
scattering angles. The light source is spatially filtered by 
aperture 1 and lens 1. The spatially filtered beam is then 
focused, by lens 2, into a spherical sample cell (see figure 
75) which contains the flowing particle dispersion. If the 
source beam already has appropriately attenuated components at 
higher divergence angles, then Lens 1 and aperture 1 can be 
eliminated and lens 2 can focus the source directly into the 
sample cell. The sample cell should have spherical shaped 
windows (also see Figure 75) to minimize the focal shift of 
the light beam focal spot, due to changes in the refractive 
index of the dispersant and to reduce Fresnel reflections. If 
this focal shift or Fresnel reflections are not a problem, 
planar windows can be used on the sample cell as shown 
previously. The scattered light from the particles is focused 
through aperture 2 by lens 3. Aperture 2 is in the image plane 
of the focal spot of the source inside of the sample cell. As 
shown previously in Figures 43 and 45, aperture 2 will 
restrict the size of the scattering interaction volume which 
can be seen by detector array 1 so that the probability of 
detecting more than one particle in the scattering interaction 
volume is small- The light passing through aperture 2 is 
detected by a detector array, as shown in Figure 79, for 
example. The detector array contains 18 separate detector 
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elements (numbered 1 through 18 in Figure 7 9) . Element 1 
measures the approximate optical flux of tKe unscattered 
light. Element 2 measures the low angle scatter for all 
scattering planes. Elements 3 through 18 measure the scattered 
light in 8 different scattering planes (2 detection sides per 
scattering plane for a total of 16 detectors) , For example, 
detector element 4 and detector element 12 measure the 
positive and negative scattering angular ranges for a single 
scattering plane. Actually, each scattering plane is the sum 
of scattering over a small range of scattering planes, around 
the center scattering plane for that wedge segment. Sometimes 
the positive and negative angular ranges will be the same, if 
the intensity is uniform across the particle and if the 
crossection of the particle in the scattering plane has 
rotational symmetry about an axis perpendicular to the 
scattering plane. Then only detectors 3 through 10 (half of 
the detector array) would be needed to cover all of the 
scattering planes. This source uniformity could be insured by 
using an appropriate attenuation profile across the beam at 
aperture 1 or diffractive beam shaper as discussed previously. 
If the intensity uniformity of the source focused spot or 
particle crossectional symmetry cannot be insured, the 
positive and negative sides of each scattering plane should be 
measured separately as shown in Figure 79. The scattering 
angle range for detector 1 is limited by the size of lens 3. 
Typically, a single lens can measure up to scattering angles 
of approximately 60 degrees. The nonspherical mirror segment 
collects light scattered at higher angles and focuses this 
light through aperture 3, which is also in an image plane of 
the source focal spot in the sample cell. The shape of the 
nonspherical mirror is designed to minimize the aberrations 
between the source focal plane in the sample cell and aperture 
3. For example, the nonspherical mirror could be a segment of 
an ellipsoid of revolution, where the source focus in the 
sample cell and aperture 3 are each located at different foci 
of the ellipsoid. This aperture 3 defines a restricted 
scattering volum'e as" shown previously in Figures 43 and 45. 
The light passing through the aperture is projected onto a 
second detector array, which is similar to that shown in 
Figure 79. However, elements 1 and 2 will not be needed for 
detector array 2 because they are only effective at low 
scattering angles, where the scattered intensity does not 
change significantly for various scattering planes. Detector 
arrays 1 and 2 are oriented so that the bisector of element 3 
is parallel to the particle flow direction. Then each element 
on detector array 2 will provide the scattered light signal at 
higher scattering angles for the same scattering plane of the 
corresponding detector element in detector array 1 . For each 
scattering plane, the signals from elements 1 and 2, the 2 
elements in that scattering plane from detector array 1, and 
the 2 elements in that scattering plane from detector array 2 
will determine the ''effective dimension'' in that plane for 
that particle. Since these calculated ^^effective dimensions" 
are not totally independent of each other, they must be 
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calculated from a set of simultaneous equations, one equation 
from each scattering plane. However, the advantage using many 
scattering planes is that the directions of the minimum and 
maximum dimension can be found quickly by comparing the ratio 
of the high and low angle scattering for each scattering 
plane. Fewer scattering planes would require much more 
computation time to solve for the dimensions of a randomly 
oriented particle, using iterative inversion of the equations. 
However, when many scattering planes are measured, the major 
and minor axes, and orientation, of the randomly oriented 
particle can be found quickly by inspection of flux ratios 
(see later in Figure 8 9) . One of the scattering planes should 
be parallel to the flow of the particles, because the major or 
minor axis of each particle is more likely to be parallel to 
the flow direction, particularly in accelerating flow which 
may occur from a crossectional area change in the sample cell. 
This crossectional area change may be designed into the flow 
path to provide the flow acceleration and particle orientation 
parallel to one of the measured scattering planes. 

In order to increase the range of dimensions which can be 
measured, more scattering angular ranges must be measured. For 
example, Figure 8 0 shows a detector array which measures two 
scattering angle ranges for each of the positive and negative 
scattering sides of each scattering plane. For example, 
elements 4 and 2 0 are low and high angle ranges for the 
positive scattering side and elements 12 and 28 are low and 
high angle ranges for the negative scattering side of the same 
scattering plane. These types of detector arrays are easily 
fabricated in custom silicon photodetector arrays. When very 
staall particles are measured, silicon photodetector arrays may 
not have sufficient signal to noise to detect the very small 
scattering intensities. In this case, photomultipliers or 
avalanche photodiodes may be used, but at much greater expense 
for manufacture of custom arrays. One solution to this problem 
is to replace the custom detector array with an array of 
dif f ractive, Fresnel, or binary lenses and rout the light from 
each element of the diffractive optic array to a separate 
element of a commercially available (inexpensive) linear or 2- 
dimensional detector array, which could be made of PMT 
elements. However, it is claimed that any diffractive optic 
array could be replaced by a detector array with elements of 
the same shape, if that option is affordable. Arrays of 
conventional curved surface lenses, diffractive lenses, 
Fresnel lenses, or binary lenses are all included in the terms 
optic array or diffractive optic array used in this 
disclosure. The diffractive optic array would consist of a 
separate diffractive lens structure covering the aperture of 
each detector element shape in Figures 79 or 80. Each 
diffractive lens would have a separate optical axis. 
Therefore, each diffractive lens element (or segment) would 
focus the scattered light captured by the aperture shape of 
that lens element to a separate point behind the lens array. 
Figure 81 shows a lens array, where each detection element 
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section contains a lens with a different optical axis . This 
idea will work for any types of lens arrays: spherical, 
nonspherical, diffractive, binary, and Fresnel lenses. Figure 
81 shows a front view of the diffractive or binary lens array, 
where the curved lines inside each detection element segment 
represent the diffractive optic structure, whose center is the 
optical axis of that lens element. The optical axis, of the 
element corresponding to detector element 1 in Figure 8 0, is 
approximately in the center of that element. The optical axis, 
of the element corresponding to detector element 2 in Figure 
80, is located off of center of the array to shift the optical 
axis away from that of element 1. 

To demonstrate this concept, consider the elements with 
optical axes marked with ''x"% in the front view (Figure 81) 
and side view (Figure 82) of Figure 81, and the corresponding 
element numbers in Figure 80. Figure 82 shows the paths for 
light rays which pass through aperture 2 in Figure 78. The 
aperture in- Figure 82 could be aperture 2 or aperture 3 in 
Figure 78- The lens element for detection element 1 collects 
the light from the unscattered light beam and focuses that 
light into fiber optic 1, as shown in figure 82. The annular 
lens element for detection element 2 collects, the scattered 
light at the lower scattering angles and focuses that light 
into fiber optic 2, as shown in figure 82. All of the light 
from the annular segment of element 2 is focused to one fiber 
optic, because that annular section is an annular section of 
one lens element with an optical axis which is shifted from 
the center of the lens array. Figure 82 also shows the 
scattered light focused by, other lens elements in the same 
scattering plane into fiber opticis 3, 4, 5, and 6. Each lens 
element has a separate optical axis so that light passing 
through aperture 2 (Figure 78) will be ^ focused into a separate 
fiber optic for each lens array element, which are shaped to 
collect the scattered light over the appropriate range of 
scattering angles for each scattering plane. Each lens element 
focuses "the light "from that element into a separate fiber 
optic, which carries that light to a separate element of a 
detector array, which may be any type of detector array, 
including 2 dimensional array or linear arrays . Each fiber 
optic might also carry light to separate detectors, which are 
not in an array. Therefore, the detector array does not need 
to have the shapes of the lens elements so that commercially 
available (non-custom) detector arrays can be used. And also 
the detector elements can be much smaller than the lens 
elements, providing much lower noise and lower cost. The cost 
of fabricating the custom lens array is much less than the 
cost of fabricating a custom photomultiplier detector array or 
silicon detector array. Spherical lens arrays can be molded 
into plastic or glass. And diffractive lens arrays can be 
molded or patterned (lithography) into nearly planar plastic 
or glass plates. The choice of positions for the optical axis 
for each element should be optimized to reduce optical 
aberrations. In some cases, the optical axis may lie outside 
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of the lens element. Also, the optical axes could be arranged 

so that the focused array of spots conforms to the 
configuration of a 2 dimensional detector array so that the 
spots can be focused directly onto the detector array, 
eliminating the fiber optics. This could also be done with 
linear arrays by creating a linear array of optical axes in 
the lens array. 

The same technique can be used for aperture 3 by replacing 
detector array 2 (see figure 78) with a lens array whose 
elements are coupled to a separate detector array through 
fiber optics (or directly coupled to the detector array, 
without fiber optics, as described above) . In this case, 
elements 1 and 2 may not be needed, because the scattered 
light distribution from the non-spherical mirror segment has a 
hole in the center where low scattered light is not captured. 
The high angle scattering should be separated into different 
scattering planes due to the high degree of asymmetry in the 
scattering pattern at higher scattering angles, so the 
scattered light in annular segments would not be useful. 

This lens array idea is most effective for large numbers of 
detection elements. For smaller numbers of elements, each 
element could have a separaite wedge prism behind it to divert 
the light to a lens which would focus it onto a particular 
fiber optic or detector element. But still the point is to 
eliminate the need for a custom detector array, to reduce the 
detector element size to reduce noise, and to allow use of 
highly sensitive detectors such as photomultipliers , which 
have limited customization. 

Quadrant detectors are commercially available for most 
detector types, including silicon photodiodes and 
photomultipliers. Figure 83 shows a method to use two quadrant 
detectors to measure scattered light from 8 adjacent 
scattering planes by using a mask which is positioned on top 
of- the- quadrant detector. The first quadrant detector replaces 
each of the detector arrays in Figure 78. The second quadrant 
detector captures scattered light at the same distance from 
the aperture 1 (or aperture 2) as the first quadrant detector 
by diverting a portion of the beam with a beam splitter placed 
on front of the first quadrant detector. Each mask is designed 
and each quadrant detector is oriented as shown in Figure 83, 
so that all 8 scattering planes are measured by the two 
quadrant detectors. In this way, two inexpensive quadrant 
detectors can measure the equivalent scattering planes 
measured by one expensive custom 8 element custom detector 
array. Typically, the central portion of the mask is designed 
to block the unscattered light from the source beam and to 
define a minimum scattering angle for each detector element. 
The detection concepts described in Figures 79, 80, 81, 82, 
and 8 3 can also be used in the other particle shape measuring 
optical systems, which were described previously in this 
disclosure . 
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Figure 84 shows a diffractive optic, where different segments 
consist of linear diffractive gratings which are designed to 
diffract nearly all of the light into one diffraction order. 
This diffractive optic is used in a hybrid 

dif f ractive/conventional lens system as shown in figure 85. 
The diffractive optic is normal to the optical axis of the 
optical system. The aperture in Figure 8 5 is the same as 
aperture 2 in figure 78. And this design can also be used with 
aperture 3 in Figure 78, without the need for the central hole 
and annular ring segment^ because aperture 3 passes only high 
angle scatter. Consider the case with aperture 2 in Figure 78. 
The scattered light is collected by the lens in figure 85. The 
diffractive optic (as shown in figure 84) in the back focal 
plane of the lens, diffracts the light incident on each 
segment of the diffractive optic into various directions. When 
the aperture hole is small and only single particles are 
measured, the diffractive optic does not need to be placed in 
the back focal plane of the lens, but the diffractive lens 
array must be scaled to match the marginal rays of the lens in 
the plane in which it resides . The light from each segment is 
nominally focused to the image plane of the aperture. However, 
the light from each segment is focused to a different position 
on that plane, because the linear grating structure in each 
segment is oriented in a different direction from the grating 
structures in other segments. Each wedge segment of the 
diffractive optic is split into 2 sections, each with a 
different grating line spacing and diffraction angle. These 2 
wedge shaped segments collect scattered light in 2 different 
angular ranges in each of various scattering planes. For 
example, scattering planes 1, 3, and 5 are shown in Figure 84. 
In addition, there is central hole and a surrounding annular 
region. The central hole passes the unscattered light beam and 
the surrounding annular region collects scattered light at 
very low scattering angles in all scattering planes. At low 
scattering angles, the scattered light is not strongly 
dependent upon scattering plane and this annular segment 
-provides a measure of low angle scattering to be used as a 
basis for each of the scattering planes. Figure 85 shows the 
light rays in one of the scattering planes. Figures 84 and 85 
show the case where the spatial frequency of the lower 
scattering angle segments is higher than the spatial frequency 
of the higher scattering angle segments. Hence the separated 
beams cross each other after the diffractive optic, because 
the diffraction angle of the lower scattering angle segments 
is higher than that for the higher scattering angle segments. 
This could also be reversed, where the spatial frequencies are 
higher for the higher scattering angle segments. In this case, 
the focused beams would still be separated at the plane of the 
fiber optics, but the beams may not cross each other. 

The shape of detector array or optic array elements is not 
limited to wedge shape. Other shapes such as linear shapes 
shown in Figure 33 could be used. 
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This method can also be used to measure ^^equivalent particle 
diameter'' without any shape determination. In this case, a 
diffractive optic as shown in Figure 8 6 could be used. This 
diffractive optic consists of annual rings which define 
multiple scattering angle ranges. The signals from these 
annual rings are independent of particle orientation and will 
only provide ^^equivalent particle diameter". 

When using a photomultiplier (PMT) , one must prevent the 
detector from producing large current levels which will damage 
the detector. This damage could be avoided by using a feedback 
loop which reduces the anode voltage of the PMT when the anode 
current reaches damaging levels. In order to avoid non-linear 
behavior over the useful range of the detector, the change of 
anode voltage should be relatively sharp at the current damage 
threshold level, with, very little change below that level. The 
response time of the feedback loop should be sufficiently 
short to prevent damage. The feedback signal could also be 
provided by a premeasuring system as shown in Figures 49, 49B, 
and 51. 

Another important point is that any of the static scattering 
techniques described in this disclosure can be- applied to 
particles which are prepared on a microscope slide, which is 
scanned through the interaction volume instead of flowing a 
particle dispersion through the interaction volume. This 
provides some advantages: the particles are confined to a thin 
layer reducing the number of coincidence counts and the 
detection system could integrate scattered light signal for a 
longer time from smaller particles, to improve signal to 
noise, by stopping the scan stage or reducing the scan speed 
when smaller particles are. detected (this can also be 
accomplished by slowing the flow rate in the dispersion case) . 
However, preparing a slide of the particles for analysis, 
greatly increases the sample preparation time and the 
potential for sample inhomogeneity . When a cover slip is 
placed onto the dispersion, the smaller particles are forced 
farther from their original positions, distorting the 
homogeneity of the sample. 

The scattering signal currents from elements on these detector 
arrays are digitized to produce scattered signal vs. time for 
each detector element- The data from each detector element is 
analyzed to produce a single value from that element per 
particle. This analysis may involve determining the time of 
maximum peak of the detection element with largest scatter 
signal and then using the same time sample for all of the 
other detection elements. Also the peak can be integrated for 
all detection elements to produce a single value for each 
element. Also the methods described previously for pulse 
analysis can be employed, including the methods (Figure 28) 
for eliminating events which are not particles . The final 
single value for each detector element represents the 
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scattered light flux collected by that element. In the 
following analysis, the following definitions will be used: 

The integral of F(x) between x == xl and x = x2 is given by: 

INT(F(x), xl, x2) 

The sum of terms of Fi(x) over index i from i = n to i = m is 
given by: 

SUM(Fi(x), i=n, i=m) 

Where Fi(x) is the ith term 

Each scattering angle corresponds to a radius, measured from 
the center of the source beam, in the detection plane. For the 
case shown in figure 78, z is the distance between aperture 2 
and detector 1 (or between aperture 3 and detector array 2) . 
Then the relationship between scattering angle 9 and the 
radius r on the detector is given by: 

r = z* tan(M0) 

where M is the angular magnification of the optical system 
(lens 3 for detector array 1 and the nonspherical mirror 
segment for detector array 2) . In the case of aperture 3 and 
detector array 2, the 9 in the above equation is related to 
the actual scattering angle through a simple equation which 
describes the angular transformation of the nonspherical 
mirror segment . 

Consider the case of rectangular particles, with dimensions 
■da and db . Then the scattered flux (flux integrated over the 
pulse, peak flux value of pulse, etc.) collected by the ith 
detector element can be described by: 

Fij = i]siT( I(da, db, a , 0i, r) *Wij (r) , rli j , r2ij ). 

I is the scattered intensity. 0i is the bisecting angle of the 
intersection of the ith scattering plane with the detector 
plane in the detector plane, as illustrated in Figure 8 4 for 
the S"""^ scattering plane. The scattering plane (which contains 
the scattered ray and the incident light beam) is 
perpendicular to the detector plane (which is the plane of 
Figures 7 9, 8 0, and 8 3 for example) . Wij (r) is the weighting 
function of scattered light as a function of r for the jth 
detector aperture, with limits between rlij and r2ij, of the 
detector for the ith scattering plane. 

Let Fl be the flux measured detector element 1 in Figure 80 
and let F2 be the flux measured detector element 2 in Figure 
80. Then the following sets of simultaneous equations can be 
formed to solve for the dimensions da , db, and random 
orientation, a, of the particle, where Ri is a ratio of two 
fluxes . 
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Fi = 


Fij 


Ri = 


Fi j /Fik 


Ri = 


Fij /F2 


Ri = 


Fij/Fl 


Ri = 


Fij/Fkj 


Ri = 


FijA/FijB 



Recognize that i indicates the ith scattering plane and ij 
indicates the jth detector in the ith scattering plane, Fij A 
and FijB are the corresponding detector elements from two 
different detector arrays (A and B) , each with a different 
Wij (r) , as shown in detector pairs A and B in Figure 74, and 
diffractive optics 1 and 2 in Figure 87. This equation also 
holds for more than two detectors (Ri = FijC/FijB, etc.). 
Also, for all other Ri equations, Wij can be the same or 
different for the numerator F and the denominator F in the 
ratio. The above simultaneous equations can be formed from the 
Fij values measured from each particle. These equations are 
solved for da, db, and oc by using various search, 
optimization, and regression methods. In most cases, the 
equations will be non-linear functions of these unknowns, 
requiring iterative methods for solution. The computation of 
theoretical values for scattered intensity I (da, db, a , 0i, 
r) , of nonspherical particles, requires long computer time. 
This is particularly problematic as this computation must be 
accomplished for each counted particle. Also numerical 
integration of these functions to produce Fij values during 
the iterative optimization process requires far too much 
computer time. This computer time can be reduced by fitting a 
series of explicitly integratible functions to each of the 
theoretical I (da, db, oc , 0i, r) functions. For example 
consider the following power series form for the I (da, db, ot , 
0i, r) and Wij (r) functions: 

Wij = SUM(Qp(j)^( r^p) , p=0, p=pmax) 

li = SUM(Cm(da, db, a , 0i)* (r^m) , m=0, m=mmax) 

Where x^p == x to the pth power and x*y = product of x and y. 

Fij = INT( I (da, db, a , 0i, r)^Wij(r), rli j , r2ij) 

Fij = INT (SUM (Cm(da, db, oc , 0i) * (r-^m) , m=0, m=mmax) * 
SUM(Qp ( j ) * (r^p)-, p=0, p-pmax) , rlij , r2ij) 



103 



wo 2005/091970 



PCT/US2005/007308 



Fij = SUM(Bq(da, db, a , 0i)*( r2ij'^ (q+1) ) / (q+l) ) , q=0, q - 
pmax+mmax ) - 

SUM(Bq{da, db, oc , 0i) * ( r li j ^ (q+1 ) ) / ( q+1 ) ) , q=0, q = 
pmax+mmax ) 

Then the previously listed sets of simultaneous equations can 
be formed from these equations for Fij . Where Bq are 
coefficients which are products of values of Cm and Qp, which 
are all known functions of da, db, a , and 0i . This concept 
can easily be extended to other particle types (just assume da 
and db to be the major and minor axes for an elliptical 
particle) and particles with more dimensions, such as 
pentagons, etc. In each case, the model must expand to account 
for the added dimensions dc, dd, de, etc. 

Fij = INT( I(da, db, dc, dd, de, , 0i, r ) -^Wi j (r ) , rli j , 

r2ij) 

In order to solve for particle types with larger number of 
dimensions (i.e. octagons., etc.), sufficient scattering planes 
and detectors must be used to provide a fully determined set 
.of simultaneous equations. In other words, the number of 
equations should be greater than or equal to the number 
unknowns, which include the dimensions, da,db,dc, ...etc. and 
the particle orientation, a. However, as indicated before, the 
solution of these equations can be computationally time 
consuming. The measurement of scattered light in many 
scattering planes can reduce the computational time, because 
the extrema of the dimension function can be found quickly. 
For example, take the case of the rectangular particle, where 
da, db, and of could be solved for with only measurements in 3 
scattering planes. The solution of these 3 equations may^ 
require iterative search and much computer time per particle. 
However, if scattering is measured in many more scattering 
planes, the major and minor axes of the rectangl.e can be 
determined immediately, eliminating the requirement for 
determining-"the particle orientation a. or reducing the range 
of a for the search solution. For example, if we interpolate a 
plot of Ri vs. 0i (or Rijk vs. i) , we will obtain a, function 
(or functions) with a maximum at 0max and a minimum at 0min, 
as shown in Figure 89. If Ri is the ratio of the high angle 
flux divided by the low angle flux in the ith scattering 
plane, then the direction 0min will provide major axis and the 
direction 0max will provide minor axis of the particle; and 
the R values for those directions will provide the particle 
dimensions in those directions. This technique can be extended 
to particles with more dimensions, as shown in Figure 90, 
where multiple extrema indicate the 0 direction and effective 
particle dimension in that direction. Since these dimensions 
may not be completely independent, this computation may 
require iterative minimization of the RMS error between 
Fij (measured) and Fij (theoretical) using various search, 
optimization, and regression methods. These methods use a 
first quests to the unknowns (all the dimensions and a) from 
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which the Fij (theoretical) is calculated," Then the difference 

between Fij (measured) and Fij (theoretical) is used to refine 

the next guess for the unknowns from which the next 

Fij (theoretical) is calculated. And this loop is repeated 

until the fit between all Fij (measured) ' s and 

Fij (theoretical) ' s is sufficient, where i=scattering plane 

index and j=detection array element or lens array element 

index in a certain scattering plane. This iterative loop is 

run many times until one of the following RMS errors are 

minimized . 

Error = SUM ( (Fij (measured) -Fij (theoretical) ) ^2) summed over 
all ij 

Error = SUM ( (Rij k (measured) -Rij k (theoretical) ) ^^2) summed over 
all ijk 

All of the detector, mask, and diffractive optic 
configurations shown in this disclosure are only examples. 
This disclosure claims the measurement of any number of 
scattering angular ranges in each of any number of scattering 
planes, as required to determine the shape and size of each 
particle . 

All segmented detector arrays or lens arrays could be replaced 
by 2-dimensional detector arrays. In this case the inverse 2- 
dimensional Fourier Transform of the spatial distribution of 
detector element flux values would produce a direct 2- 
dimensional function of the particle shape, in the Fraunhofer 
approximation. The dimensions of the contour plot (perhaps by 
choosing the 50% contour of the peak contour) of this 2- 
dimensional inverse Fourier Transform function will provide 
the outline of the particle directly, for particles which are 
modeled by the Fraunhofer approximation. However, the size 
range of this type of array (per number of detector elements) 
is not as efficient as the wedge shaped scattering plane 
arrays described - previously , because the detector elements in 
the 2-dimensional array are all the same size. When the 
particles become large enough to produce scattered light in 
only the lowest scattering angle detector element, the size 
determination becomes difficult. This problem could be solved 
by using a custom array, where low scattering angle elements 
are smaller than larger scattering angle elements. This 
progression of element size with increasing scattering angle 
can also be accomplished with a linear array which follows a 
lens or diffractive optic with non-linear distortion. The 
lens or diffractive optic distorts the scattering pattern so 
that the pattern is spread out near the center and compressed 
more at higher radii (larger scattering angles) in the 
pattern. In this way, detector elements closer to the center 
of the pattern will subtend a smaller scattering angular width 
than elements farther from the center. This would increase the 
size range of the detector array and still allow use of 
standard CCD type linear arrays. However, due to the limited 
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dynamic range of most CCD arrays, a singxe t^ixiT or ul-iicj. a-do-^^i^ 
dynamic range pre-sensor, placed upstream of the CCD array, 
could provide some indication of scattering signal level 
before the particle arrives in the view of the CCD array, 
similar to the systems shown in Figures 67 and 68. Then the 
source power level or CCD electronic gain could be adjusted 
during the CCD data collection to optimize the signal to noise 
and fully utilize the range of the analog to digital converter 
on the CCD. 

This wider size range can also be obtained by using different 
weighting functions (Wij) for two different detector elements 
which view the same range of scattering angle in the same 
scattering plane, as described previously. As long as the Wij 
functions are different for the two measurements, the ratio of 
those scattered flux values will be size dependent over a 
large range of particle size and will be relatively 
insensitive to position of the particle in the beam. The Wij 
function can be implemented in the detector element shape, as 
shown in Figure 7 6, and/or by placing an attenuation mask, 
which has varying attenuation along the direction of changing 
scattering angle, over the detector elements of one the arrays 
or with different Wij functions on each of two or more arrays 
which see the same scattering angle ranges and scattering 
planes. This idea can be implemented in the system shown in 
Figure 78, using the detection module shown in Figure 87. The 
aperture in Figure 87 could be either aperture 2 or aperture 3 
in Figure 78. In this case, the scattered light (and incident 
beam when used with aperture 2) is focused by a lens, through 
a beamsplitter, onto two optical arrays, which are diffractive 
arrays using conventional lens/ diffractive array hybrid as in 
Figure 85 (or conventional molded or diffractive optics if no 
single conventional lens is used) . The lens arrays (examples 
are shown in Figures 79, 80,, 84, 86) are used to divert the 
scattered light from different array elements to different 
positions in the image plane of the beam focus in the cell so 
that the flux from each element can be collected separately by 
fiber optics or detectors. The lens array in Figure 87 is 
similar to that in Figure 79, just for illustration, but any 
of the previous lens arrays could be used in this design. 
There is a mask on front of each lens array which provides 
selection of center hole (element 1 in Figure 78) or the 
annular ring (element 2 in Figure 78) . This mask also provides 
light attenuation which varies along the radius (or scattering 
angle) of the mask, as shown by the transmission function, for 
example only, in Figure 88. This transmission function can 
take many forms, but as long as the radial transmission 
function of mask 1 is different from the radial transmission 
function of mask 2, the ratio of flux from corresponding lens 
elements in diffractive optic 1 and diffractive optic 2 (as 
captured by the fiber optics or detector elements) will 
provide dimensional information for the particle in the 
scattering plane of that element. For example transmission 
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functions for mask 1^ Tlij, and for mask 2, T2ij, could have 
many cases including: 

Case 1: Tlij = 1 (constant transmission), T2ij = r 

Case 2: Tlij = T2ij = r^2 

where r is the radius along the lens array 

Then the effective weighting functions are the product of the 
transmission functions and the weighting function, Wijs, of 
the segment or element shape in the detector array or optic 
array. For wedge shaped segments the shape weighting function 
is : 

Wijs = r 

Hence the effective weighting functions are: 
Wlij = Wlijs ,*Tlij 
W2ij = W2ijs .*T2ij 



R12ij = Flij/F2ij 

Where Flij = flux from the jth detector aperture in the ith 
scattering plane of detector array 1 or optic array 1 

Where F2ij - flux from the jth detector aperture in the ith 
scattering plane of detector array 2 or optic array 2 

Many combinations of transmission functions will work. This 
invention disclosure claims the use of any two different Wlij 
and W2ij functions (using transmission and/or shape weighting 
functions) and using flux from corresponding array elements, 
one with Wlij and the other with W2ij . Diffractive optic array 
1 and diffractive optic array 2 are identical and they are 
rotated so that each array segment collects scattered light 
from the same scattering plane as the corresponding element in 
the other diffractive optic. Then ratio (R12ij) of light flux 
values from these two elements are is used to determine the 
"^^effective dimension^' of the particle in that scattering 
plane. However, as described previous, these dimensions are 
not independent and the full set of simultaneous Rl2ij 
equations must be used to solve for the actual dimensions. 

In general, we have two cases for measuring 2 dimensional 
scattering distributions. The detector array can be a set of 
radial extensions in various scattering planes (as in Figure 
84) or a 2-dimensional array. In first case, we define the 
array as a function of radius, r, and angle 0 in the detector 
array or optic array plane, as shown in Figure 84. In the 
second case the array is a standard 2-dimensional array (i.e. 
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CCD array) with elements arranged in columns and rows which 
are defined as a function of x and y. The choice between these 
cases will be determined from the available computer speed, 
shape constraints, and particle size range. In order to 
measure particle shape statistics of a particle sample, an 
enormous number of particles must be characterized. For 
example, consider the following case. We are to count 
particles and characterize them into 25 different dimensional 
classes. Each class could vary by size or shape. If we want to 
measure the content of each class to within 5% we need to 
measure 400 particles per class, if the counting process is 
Poisson. Therefore we need to measure and characterize 10000 
particles. If we want the entire analysis process to take 5 
minutes, then we have 30 milliseconds per particle to measure 
and digitize the scattered light, and solve equations for the 
shape and dimensions. We could do a full global search for the 
particle orientation and dimensions using the equations given 
above. All of the parameters of those equations can be solved 
from the 2-dimensional scattered intensity distribution in the 
plane of the detector array or optic array. The theoretical 2- 
dimensional scattered intensity distribution is calculated 
using known methods, such as T-matrix and Discrete Dipole 
Approximation, (see ^'Light Scattering , by Nonspherical 
Particles'% M. Mishchenko, et al.). Then this theoretical 
intensity distribution is integrated over the areas of each 
element of detector array or optic array. These values for Fij 
are calculated for various* particle orientations, a, and 
various particle dimensions. A global search routine would 
search among these theoretical Fij sets to find the one which 
best fits to the measured set of Fij . Then the dimensions of 
that set would be accepted for that particle. While this 
process will take far longer than 30 milliseconds per 
particle, it would produce the most accurate determination of 
particle dimensions and shape. This full search method is 
claimed by this invention because in some cases, users will be 
willing to collect data on the optical system computer and 
then transfer the' raw data off-line to a set of parallel 
processors which continuously crunch data sets 24 hours per 
day. And also future computers will be capable of doing these 
computations in the required time. However, for the users, who 
want their results in real-time (i.e. within 5 minutes) with 
present personal computers some of the shortcuts, as 
described previously, must be employed. In the case of process 
control, users need results quickly in order to adjust the 
process parameters in nearly real-time. And also, many times 
the determined shape and size does not need to be precise 
because these results only need to correlate to the quality of 
their product. In general a particle with multiple flat edges 
will produce a scattering pattern with a radial projection for 
each edge into a scattering plane which is perpendicular to 
that edge. So when the ratio data, Rijk, is plotted vs. 0i (or 
vs. scattering plane) as shown in Figure 90, one will obtain a 
maximum in Rijk for each side of a multi-sided particle. So 
evaluation of Rijk vs. 0i will very quickly determine the 



108 



wo 2005/091970 



PCT/US2005/007308 



orientation of the particles sides and the values of Rijk are 
then used in a more limited global search routine, which does 
not have to search over all possible orientations of the 
particle and over all possible number of particle sides. This 
will dramatically reduce the search time. Also, the 
'Mimension" of the particle in the direction perpendicular to 
each side can be determined approximately from Rijk values in 
the scattering plane which is perpendicular to that side. 
These values could be used directly to determine ''approximate 
dimensions"' of the multisided particle, without a search 
algorithm, because the approximate dimension can be calculated 
directly from a theoretical function of Rijk. As the particle 
model becomes simpler (i.e. rectangle) the orientation and 
dimensions are calculated immediately by finding the 0max and 
0min in the Rijk vs. 0i, and then using the Rijk values in 
those two scattering planes to calculate the dimensions in 
those scattering planes. If only a few scattering planes are 
measured, the Rijk function could be interpolated or fit to a 
theoretical function for the rectangle case, to calculate the 
exact orientation, which may be between two adjacent 
scattering planes. 

These concepts can be combined with imaging systems to record 
the image of selected particles after they have passed through 
the interaction volume. An imaging system could be placed 
downstream of the interaction volume, with a pulsed light 
source which is triggered to fire at the correct delay, 
relative to the scatter pulse time, so that the particle has 
flowed into the center of the imaging beam during image 
capture. The pulsed light source has a very short pulse period 
so that the moving particle has very little motion during the 
illumination and image capture on a CCD array. The particle is 
imaged onto the CCD array at high magnification with a lens 
(microscope objective would be a good choice) . In this way, 
particles which meet certain criteria, can be imaged to 
determine their morphology. 

The alignment of aperture 2 and aperture 3 in Figure 78, and 
other related figures, could be accomplished by running a 
medium concentration sample of sub-micron particles through 
the sample cell. The concentration is chosen so that a large 
number of particles are in the interaction volume at the same 
time so that a constant large scatter signal is seen on the 
detectors. Then the x, y, and z position of each aperture is 
adjusted to maximize the scatter signal on the detectors. 

Most of the concepts in this disclosure can accommodate 
aerosol particle samples, by removing the sample cell and by 
flowing the aerosol through the interaction volume of the 
incident beam. The effective scattering angles may change due 
to the change in refractive index of the dispersant. 

In some cases, the dynamic range of detectors will not be 
sufficient to cover the entire range of scatter signals from 
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the particles. In particular, particles in the Rayleigh 
scattering range will produce scatter signals proportional to 
the 6^^ power of the particle diameter. Photomultipliers can 
also be damaged by large levels of light. Figure 93 shows an 
optical system which uses upstream scatter measurement to 
control the laser power or detection gain (or anode voltage) 
for a system down stream in particle flow to protect 
photomultipliers to maximize the signal to noise, or to avoid 
detector saturation. Two light sources, light source 1 and 
light source 2, are combined by beamsplitter 1 and focused 
into the center of the sample cell to two different locations 
along the particle flow path. Light source 2 could be 
magnifiied to produce a larger spot size in the sample cell for 
detecting larger particles than the spot from light source 1. 
However, the other purpose of light source 2 is to detect an 
oncoming particle before it reaches the focused spot from 
light source 1. The spot from light source 2 is upstream from 
the light source 1 spot. The light from both sources and the 
light scattered from both sources pass through lens 3, which 
images both light spots to the planes of three apertures, 1, 2 
and 3. Apertures 2 and 3,^ which receive light reflected from 
beamsplitter 2, block light from the interaction volume of 
light source 1 but pass light from the interaction volume of 
light source 2. Likewise, aperture 1^, which receives light 
transmitted by beamsplitter 1, blocks light from the 
interaction volume of light source 2 but passes light from the 
interaction volume of light source 1. Therefore multi-element 
detector 1 sees only light scattered from light source 1; and 
multi-element detector 2 and detector 3 see only light 
scattered from light source 2 . These multi-element detectors 
can also be replaced by the optic array systems described 
previously. Multi-element detector 2 is operated at much lower 
sensitivity than Multi-element detector 1, which is operated 
at , maximum sensitivity to detect the smallest particles. 
Whenever a particle, which would saturate and/or damage Multi- 
element detector 1, passes through light source 2 spot, 
detector (s) from Multi-element detector 2 or detector 3 will 
measure the larger amount of scattered light and trigger a 
circuit to lower the power level of light source 1 or lower 
the, gain (or anode voltage) of Multi-element detector 1 so 
that Multi-element detector 1 will not be saturated and/or 
damaged when that same particle passes through the interaction, 
volume for light source 1. The level of adjustment can be 
variable depending upon the light level measured by Multi- 
element detector 2 or detector 3. Light source 1 is normally 
run at maximum power to detect the smallest particles. Light 
source 1 power is only reduced after a calculated delay time 
after Multi-element detector 2 detects a potentially damaging 
or saturating particle scatter signal- After that particle 
passes through the light source 1 spot, the source 1 power (or 
gains, anode voltages etc) is reset to maximum. The time delay 
is calculated from the spacing between the two source spots in 
the sample cell and the particle flow velocity. The multi- 
element detector arrays are in the focal plane of lens 4 or 
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lens 5. Notice that the blue light rays indicate that multi- 
element detector arrays 1 and 2 are effectively, at infinity, 
in the back focal planes of lenses 4 and 5. These lenses place 
the detectors at infinity so that the effects of finite 
pinhole size in apertures 1 and 2 do not cause smearing in the 
scatter pattern. Typically the pinhole sizes are small enough 
so that lens 4 and lens 5 are not needed. The advantage of 
using lens 2 for both sources is not only the cost of 
manufacture. This design allows the two source spots in the 
sample cell to be very close to each other^. insuring that all 
particles which flow through the light source 1 spot will have 
also previously flowed through the light source 2 spot and 
been detected by multi-element detector 2, even in the event 
of any flow anomalies in the cell. To provide a large size 
dynamic range, the sample cell spot size for light source 2 
may be much larger than the spot size for light source 1, in 
order to measure much larger particles - Then many particles 
which pass through source 1 spot will flow around source 2 
spot. In this case, a third system is added using beamsplitter 
3 and aperture 3 . The pinhole of aperture 3 is smaller than 
aperture 2 to only pass scatter from the portion of Source 2 
spot which is directly above the Source 1 spot. In this way, 
detector 3 will only see the particles which will eventually 
pass through the source 1 spot. So detector 3 is used to set 
the source 1 power level or detector 1 gain (or anode voltage) 
using the time delay described above. Multi-element detector 2 
measures the size of larger particles from a much larger 
interaction volume as defined by aperture 2. For example, the 
spot sizes in the sample cell could be 20 microns for light 
source 1 and 500 microns for light source 2, but aperture 3 
would only allow detector 3 to see scatter from a 30 micron 
portion, of the Source 2 spot, which is directly above the 20 
micron spot of Source 1. The 30 micron portion could be 
slightly larger than the 20 micron spot to accommodate slight 
flow direction misalignment, because the larger size will only 
trigger the source 1 power drop more times than needed, but it 
will guarantee that no particle scatter will saturate and/or 
damage multi-element detector 1. 

Figure 96 shows a variation of the optical system shown in 
figure 78, which provides heterodyne detection and measures 
scattered light over multiple scattering angle ranges and in 
multiple scattering planes. Some light is split off from the 
source light by beamsplitter 1 and focused into a fiber optic. 
This fiber optic passes through an optional optical phase or 
frequency shifter to provide an optical frequency shifted 
local oscillator for the detection of scattered light. This 
shifter could be an acousto-optic device or moving diffraction 
grating. The frequency shift can also be provided by a scanned 
optical phase shifter (moving mirror or piezoelectric fiber 
stretcher) whose optical phase is ramped by a sawtooth 
function to produce an effective optical frequency shift 
during each period of phase ramp. Fiber optic coupler 1 splits" 
off a portion of the light after the phase shifter and passes 
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this light to detector 3 which monitors the fluctuations in 
the light intensity (12) which may be due to laser noise or 
amplitude modulation from the optical frequency or phase 
shifter. Then the source light is finally split into two 
fibers by fiber optic coupler 2 to provide local oscillators 
for both detection systems . The light exiting the one fiber 
from fiber optic coupler 1 is expanded by negative lens 7 and 
then focused by lens 6 through aperture 2 (via beamsplitter 2) 
to be mixed with the scattered light on detector array 1 (or 
optic array 1) . Likewise light, from the other output fiber of 
fiber optic coupler 2, is expanded by negative lens 5 and then 
focused by lens 4 through aperture 3 (via beamsplitter 3) to 
be mixed with the scattered light on detector array 2 (or 
optic array 2). In this way, heterodyne detection is 
accomplished with laser noise reduction using the equation and 
method described previously; 



Idiff = Ilnb - (R/K)*I2nb = Sqrt (R^ ( 1-R) *S ) *Ioc*COS ( F*t + 
A) " 

Idiff only contains the heterodyne signal. The common mode 
noise in the local oscillator and the heterodyne signal is 
removed by this differential measurement (see the previous 
description of the method) . Heterodyne detection provides very 
high signal to noise, if the laser noise is removed by this 
equation and method. However, if the heterodyne frequency is 
only due to Doppler shift of the scattered light from particle 
motion, then the frequency of the heterodyne beat frequency 
will depend upon scattering angle and scattering plane (for 
example, the scattering plane, which is perpendicular to the 
particle flow, will show zero Doppler frequency shift of the 
scattered light) . The addition of the optical phase or 
frequency shifter provides a much higher heterodyne frequency 
which is nearly equal for all scattering angles and scattering 
planes, allowing heterodyne detection of particle size and 
shape. The only problem presented by the frequency shifter is 
that all light that hits the detector, by scatter or 
reflection, will be frequency shifted. Without the frequency 
shifter, only scatter from moving particles will contribute to 
the heterodyne signal at the beat frequency, so background 
light can be distinguished from particle scatter based upon 
signal frequency. So when the frequency shifter is used, a 
background scatter heterodyne signal should be recorded 
without particles and this background signal should be 
subtracted from the scatter heterodyne signal with particles 
present. Addition of an optical frequency shifter also 
provides a higher beat frequency and phase sensitive detection 
capability. The signals due to the particle motion and the 
Doppler effect have random phase for each particle. So the 
other advantage of the frequency shifter is that the 
heterodyne signal will have a known phase (same as the 
frequency shifter) , which could allow for phase sensitive 
detection (lock in amplifier) . As shown before, the signals 
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from this system will consist of a sinusoidal signal with an 
envelope function from a particle's passage through the 
intensity profile of the source beam. So all of the techniques 
described previously for processing these signals can be 
applied to this case. 

Also the heterodyne beat frequency will be different for each 
scattering plane. One method to eliminate this dependence is 
to design the system to have minimal particle motion induced 
phase shift by choosing appropriate angles, between the 
particle motion direction and the beam, and using the optical 
frequency or phase shifter (Figure 96 ) to provide the phase 
modulation at a high frequency, instead of the phase shift due 
to particle motion. Then all detector heterodyne signals will 
have almost the same frequency and phase, with the high signal 
to noise provided by heterodyne detection. Also the phase 
modulation signal could be used with phase sensitive detection 
(lock m amplifier) to detect the heterodyne signal. However, 
two disadvantages of this system are light reflections and 
amplitude modulation due to the phase modulator. Without the 
phase modulator, only scattered light from moving particles 
will create a heterodyne beat signal. However, with the phase 
modulator and appropriate angles between the particle motion 
direction and the beam, all light reaching the detector from 
scattering or reflections will be at the beat frequency and 
will be passed by the band pass filter. Also the phase or 
frequency modulator will produce some small amount of 
amplitude modulation in the beam, which may completely 
overwhelm the particle scattering signal, even after it is 
removed using detector 3 (figure 96) and the differential 
detection described previously. The severity of these problems 
will be determined by the characteristics of the phase or 
frequency modulator and the level of light reflections in the 
optical system. 

The system in Figure 96 could also be designed without fiber 
optics . Each fiber optic coupler would be replaced by. a 
beamsplitter and the light beams could be routed to lens 5 and 
lens 7 using mirrors and lenses to create a beam focus at 
aperture 2 and aperture 3 through beamsplitter 2 and 
beamsplitter 3, respectively. The negative lenses 5 and 7 may 
be needed to expand the beam to fill the angular range of 
light on each detector array or optic array. 

Figure 97 shows another variation of the optical system shown 
in figure 78. The source beam is split into two beams which 
cross each other in the interaction volume in the sample cell. 
The two beams will create a fringe pattern at their 
intersection, which will modulate the scattered intensity as a 
particle passes through the intersection, as shown previously 
in Figure 18. Lens 1 focuses the source light into the center 
of the sample cell. However, the light beam is split into two 
beams by beamsplitter 1. These two beams are reflected by 
mirror 1 and mirror 2 to cross in the center of the sample 
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cell. One of the beams may pass through an optical frequency 
shifter to provide a beat signal of known phase and/or higher 
frequency. The advantage of this method is that only particles 
passing through the fringe pattern at the intersection of 
these dual beams will produce signals at the beat frequency. 
This intersection, and aperture 2 or aperture 3, define a 
small interaction volume which reduces coincidence particle 
counts. Addition of an optical frequency shifter provides a 
higher beat frequency and phase sensitive detection 
capability. The signals due to the periodicity of the fringe 
pattern have random phase for each particle. So the other 
advantage of the frequency shifter is that the heterodyne 
signal will have a known phase (same as the frequency 
shifter) , which could allow for phase sensitive detection 
(lock in amplifier) . As shown before, the signals from this 
system will consist of a sinusoidal signal with an envelope 
function from a particle's passage through the intensity 
profile of the fringe pattern. So all of the techniques 
described previously for processing these signals can be 
applied to this case. In this case, the scattered signal, at 
each position on the detector array, is the square root of the 
product of the scattered intensities from each of the crossed 
light beams,. Hence the scattered light intensities, from a 
particular position on the detector array, in the simultaneous 
equations shown previously, must be replaced with the square 
root of the product of the scattered intensities from each of 
the crossed light beams at that same detector array position. 

The method shown in Figure 97 could also be implemented using 
the optical system in Figure 78, by placing a periodic mask in 
the plane of aperture 1 in Figure 78. The image of this mask 
in the sample cell would provide a periodic intensity profile 
which would modulate the scatter signal from a particle 
passing through the intensity profile- The mask could have a 
sinusoidal or square wave transmission profile as shown in a 
single section of the mask in Figure 66 (for example) . The 
mask could also be fabricated with a Barker code or other code 
which has a very sharp autocorrelation function to use 
correlation between detector signals (also see below) . 
However, in this case the scattering signals are not the 
square root of the product of scattering signals from two 
different scattering angles. The signal is only from the 
actual scattering angles defined before for detector array 1 
or detector array 2 . 

In both Figures 96 and 97, the system can be designed to 
provide beat signals on all detector elements with nearly the 
same phase and frequency. Hence high level signals can be 
multiplied times lower level signals to retrieve the lower 
level signals from noise. This method can also be used with 
the system shown in Figure 78, where each particle produces a 
single pulse. The integral of the product of the largest 
scatter signal with a smaller signal, which needs to be 
recovered from the noise, will improve the signal to noise of 
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the integral over the pulse for that smaller signal. This 
could also be accomplished by only integrating the lower 
signal while the larger signal is above some threshold, as 
described previously. This method will also work for signals 
which are not modulated. The integral of the product of single 
pulses will also improve the signal to noise of lower level 
signals when multiplied by higher level signals which have the 
same pulse shape. This could be accomplished with the 
following equation which calculates a more accurate estimate 
to the integral of 32 by using correlation with a higher level 
signal of the same functional shape (single pulse, amplitude 
modulated heterodyne signal, etc.): 

S2I = INT(S1.^S2, tl,t2) /INT(Sl,tl,t2) 

Where SI is the high level signal and 32 is the low level 
signal which has high correlation to SI, And tl and t2 are the 
start and stop times of the particle scatter pulse. The value 
of INT{S1.*S2, tl,t2) could be used directly, without 
normalization, in simultaneous equations, lookup tables, or 
search models as long as the theoretical model is calculated 
for 1NT(S1.*S2, tl,t2). 



In both Figures 96 and 97, the detector arrays can be replaced 
by optic arrays as shown in Figures 79, 80, 81 83, '84, and 86, 
for example. These optic arrays can be configured in a 
detection system as shown in Figures 82, 85, and 87, for 
example . 

All drawings of optical systems in this disclosure are for 
illustrative purposes and do not necessarily describe the 
actual size of lens apertures, lens surface shapes, lens 
designs, lens numerical apertures, and beam divergences. All 
lens and mirror designs should be optimized for their optical 
conjugates and design requirements of that optical system 
using known lens and mirror design methods . The drawings are 
designed to describe the concept; and so some beam divergences 
are exaggerated in order to clearly show beam focal planes and' 
image planes within the optical system. If these drawings were 
made to scale, certain aspects of the invention could not be 
illustrated. And in particular, the source beam divergence 
half angle must be smaller than the lowest scattering angle 
which will be measured from particles in that beam. For low 
scattering angles (larger particles) the source beam would 
have a very small divergence angle, which could not be seen on 
the drawing . 

Also in this disclosure, where ever an aperture is used to 
pass scattered light and that aperture is in an image plane of 
the scattering particle, a lens can be placed between the 
detector (s) (or optic array) and that aperture to reduce 
smearing of the scattering pattern due to the finite size of 
the aperture. The detectors could be placed in the back focal 
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plane of said lens, where each point in the focal plane 

corresponds to the same scattering angle from any point in the 
interaction volume . The detector (s) (or optic array) would be 
placed in the back focal plane of said lens to effectively 
place the detector at infinity, where the angular smearing is 
negligible, as shown and discussed in Figure 24. 

Also in this document, any use of the term ^^scattering angle" 
will refer to a range of scattering angles about some mean 
scattering angle. The angular range is chosen to optimize the 
performance of the measurement in each case. For example the 
use of the terms ^'low scattering angle" or ^^high scattering 
angle" refer to two different ranges of scattering angles, 
because each detector measures scattered light over a certain 
range of scattering angles 

Note that ail optic arrays described in this disclosure can be 
constructed from segments of conventional spherical and 
aspherical lenses, diffractive optics, binary optics, and 
Fresnel lenses. 



Ensemble particle size measuring systems gather data from a 
large group of particles and then invert the scattering 
information from the large particle group to determine the 
particle size distribution. This scatter data usually consists 
of a scatter signal vs. time (dynamic scattering) or scatter 
signal vs. scattering angle (static scattering). The data is 
collected in data sets, which are then combined into a single 
larger data record for processing and inversion to produce the 
particle size distribution. Inversion techniques such as 
deconvolution and search routines have been used. The data set 
for dynamic light scattering consists of a digital record of 
the detector signal over a certain time, perhaps 1 second. The 
power spectra or autocorrelation functions of the data sets 
are usually combined to produce the combined input to the 
inversion algorithm for dynamic light scattering to invert the 
power spectrum or autocorrelation function into a particle 
size distribution. Also the data sets can be combined by 
concatenation, or by windowing and concatenation, to produce 
longer data sets prior to power spectrum estimation or 
autocorrelation. Then these power spectra or autocorrelation 
functions are averaged (the values at each frequency or delay 
are averaged over the data sets) to produce a single function 
for inversion to particle size. Like wise for angular 
scattering, the angular scatter signals from multiple 
detectors are integrated over a short interval. These angular 
scattering data sets are combined by simply averaging data 
values at each scattering angle over multiple data sets. 

Since the inverse problem for these systems is usually ill- 
conditioned, detecting small amounts of large particles mixed 
in a sample of smaller particles may be difficult because all 
of the particle signals from the particle sample are inverted 
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as one signal set. If the signal, from only a few larger 
particles, is mixed with the signals from all of the other 
smaller particles, the total large particle scatter signal may 
be less than 0,01 percent of the total and be lost in the 
inversion process . However, in the single short data set which 
contained the larger particle's scattered light, the larger 
particle scatter may make up 50% to 90% of the total signal. 
The larger particle will easily be detected during inversion 
of these individual data sets. 

Users of these systems usually want to detect small numbers of 
large particles in a much larger number of smaller particles, 
because these larger particles cause problems in the use of 
the particle sample. For example, in lens polishing slurries, 
only a few larger particles can damage the optical surface 
during the polishing process. In most cases these larger 
particles represent a very small fraction of the sample on a 
number basis. Therefore, if many signal sets (a digitized 
signal vs. time for dynamic scattering or digitized signal vs. 
scattering angle for static scattering) are collected, only a 
few sets will include any scattered signals from larger 
particles . An algorithm could sort out all of the data sets 
which contain signals from larger particles and invert them 
separately, in groups, to produce multiple size distributions, 
which are then weighted by their total signal time and then 
combined to form the total particle size distribution. The 
data sets may also be sorted into groups of similar 
characteristics, and then each group is inverted separately to 
produce multiple size distributions, which are then weighted 
by their total signal time and then summed over each size 
channel to form the total particle size distribution. In this 
way, the larger particles are found easily and the smaller 
particle data sets are not distorted by scatter signals from 
the larger particles . ^ 

This technique works better when many short pieces of data are 
analyzed separately, because then the best" discrimination and 
detection of particles is obtained. However, this also . 
requires much pre-inversion analysis of a large number of data 
sets. The key is that these data sets can be categorized with 
very little analysis. In the case of angular light scattering, 
comparison of signal values from a few scattering angles from 
each signal set is sufficient to determine which signal sets 
include signals from larger particles or have specific 
characteristics- In the case of dynamic light scattering, the 
spectral power in certain frequency bands, as measured by fast 
Fourier transform of the data set or by analog electronic 
bandpass filters could be used to categorize data sets. 
Consider a dynamic scattering system where the scattering 
signal from the detector (in heterodyne or homodyne mode) is 
digitized by an analog to digital converter for presentation 
to a computer inversion algorithm. In addition, the signal is 
connected to multiple analog filters and RMS circuits, which 
are sequentially sampled by the analog to digital converter to 
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append each digitized data set with values of total power in 
certain appropriate frequency bands which provide optimal 
discrimination for larger particles. The use of analog filters 
may shorten the characterization process when compared to the 
computation of the Fourier transform. These frequency band 
power values are then used to sort the data sets into groups 
of similar characteristics. Since larger particles will 
usually produce a large signal pulse, both signal amplitude 
and/or frequency characteristics can be used to sort the data 
sets. The total data from each formed group is then processed 
and inverted separately to produce an individual particle size 
distribution. These particle size distributions are suinmed 
together after each distribution is weighted by the total time 
of the data collected for the corresponding group. 

The use of analog filters is only critical when the computer 
speed is not sufficient to calculate the power spectriom of 
each data set. Otherwise the power spectra could be calculated 
from each data set first, and then the power values in 
appropriate frequency bands, as determined from the computed 
power spectrum^ could be used to sort the spectra into groups , 
before the total data from each group is then processed and 
inverted separately to produce an individual particle size 
distribution. For example the rati,o of the power in two 
different frequency bands can indicate the presence of large 
particles. These particle size distributions are summed 
together after each distribution is weighted by the total time 
of the data collected for the corresponding group. This 
process could also be accomplished using the autocorrelation 
function instead of the, power spectrum of the scatter signal. 
Then the frequency would be replaced by time delay of the 
autocorrelation function and different bands of time delay 
would be analyzed to sort the data sets before creating data 
groups . 

In angular scattering, a group of detectors measure scattered 
light from the particles over a different angular range for 
each detector. These detector signals are integrated over a 
certain measurement interval and then the integrals are 
sampled by multiplexer and an analog to digital converter- In 
this case, the angular scattering values at appropriate 
angles, which show optimal discrimination for larger 
particles, could be used to sort the angular scattering data 
sets into groups before the total data from each group is then 
processed and inverted separately to produce an individual 
particle size distribution for that group. These resulting 
particle size distributions are summed together after each 
distribution is weighted by the total time of the data 
collected for the corresponding group. 

These sorting techniques can also be used to eliminate certain 
data sets from any data set group which is inverted to produce 
the particle size distribution. For example, in dynamic 
scattering, larger particles may occasionally pass through the 
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interaction volume of the optical system and produce a large 
signal with non-Brownian characteristics which would distort 
the results for the data set group to which this defective 
data set would be added. Large particles^ which are outside of 
the instrument size range, may also cause errors in the 
inverted size distribution for smaller particles when their 
data sets are combined. Therefore, such defective data sets 
should be rejected and not added to any group. This method 
would also be useful in conventional dynamic light scattering 
systems, where multiple groups are not used, to remove bad 
data sets from the final grouped data which is inverted. By 
breaking the entire data record into small segments and 
sorting each segment, the bad data segments can be found and 
discarded prior to combination of the data power spectra or 
autocorrelation functions and final data inversion. This 
method would also be useful in static ' angular scattering to 
eliminate data sets from particles which are outside of the 
instrument size range. 

In some cases, a large number of categories for sorted groups 
are appropriate to obtain optimal separation and 
characterization of the particle sample. The number of 
categories is only limited by the cumulated inversion time for 
all of the sorted groups. The total inversion time may become 
too long for a large number of groups, because a separate 
inversion must be done for each group. However, after the 
information is sorted, abbreviated inversion techniques may be 
used because the high accuracy of size distribution tails 
would not be required to obtain high accuracy in the final 
combined particle size distribution. In many cases, only two 
groups are necessary to separate out the largest particles or 
to eliminate defective data sets. 

This disclosure claims sorting of data sets for any 
characteristics of interest (not only lar.ge particles) and for 
any applications where large data sets can be broken up into 
smaller segments and sorted prior to individual analysis or 
inversion of each individual set. Then the resulting 
distributions are combined to create the final result. This 
includes applications outside of particle size measurement. 

Another application is Zeta potential measurement. Low 
scattering angles are desirable in measurement of mobility of 
particles to reduce the Doppler broadening due to Brownian 
motion. However, large particles scatter much more at small 
angles than small particles do; and so the scatter from any 
debris in the sample will swamp the Doppler signal from the 
motion of the smaller charged particles in the electric field. 
This inventor has disclosed methods of measuring Dynamic light 
scattering from small interaction volumes created by 
restricting the size of the illuminating beam and the 
effective viewing volume. When only scattered light from a 
very small sample volume is measured, the scatter from large 
dust particles will be very intermittent, due to their small 
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counr per unit volume. So the techniques outlined above can 
used to eliminate the portions of th2 signal vS timHeco^d 
which contain large signal bursts due to passage oft llral 
particle. In this way, Zeta potential measurements can be made 

?^ni°H ^^^tt^fi^? ^^gl^^ without the scattering interference 
from dust contaminants. v=j.ciiL-e 

In optical systems which need to count very small particles 
ii! h'^k'^''"'®^ "^^^^ shorter wavelengths may be preferred due 'to 
the higher scattering efficiencies (or scattering 
crossections) at shorter wavelengths. 

^5^^.?^^?^^''^ °^ angular scattering systems can drift due to 
drift of the source beam position or changes in the wedge 
between the sample cell windows. Figure 100 shows an optical 
system which measures angular scattering distribution from a 
■ ^tm'r}^J^^^^''^^°'^J^ ^ sample cell. The angular scattering 
distribution IS used to determine the size distribution of ?he 
particles. If the optical wedge between the sample cell 
windows changes or if the dispersant refractive index is 
changed, the system will go out of alignment due to refraction 
m the optical wedge. The laser focus spot will move to 
another position on the detector array, saturating detector 
elements which should be measuring scattered light only. This 
invention uses a retroref lector and a source-detector module 
to provide stable alignment against the drift sources 
described above. Lens 1, pinhole, and lens 2 form a spatial 
filter. The light beam from this spatial filter passes through 
the sample cell and dispersion. The beam and scattered light 
are then retroref lected back through the same sample cell, 
where more scatter occurs on the second pass. All of the 
scattered light is collected by lens 3, which focuses it onto 
a detector array in the back focal plane of lens 3. As long as 
the components are rigidly mounted to a common base in the 
source-detector module, the system will maintain alignment 
after initial alignment at the factory. After initial 
alignment, this system will maintain alignment over a large 
range of beam alignment drift, dispersant refractive index 
change, or sample cell wedge drift. If the sample cell drift 
or dispersant refractive index drift are not problems, the 
optics can be arranged so that the beam between lens 2 and the 
retroreflector does not pass through the sample cell, while 
the beam does pass through the sample cell between the 
retroreflector and lens 3. 

It will be understood that the embodiments of the present 
invention system and method are merely exemplary and that a 
person skilled in the art can make many modifications to the 
embodiments using functionally equivalent components. All such 
modifications, variations and alternate embodiments are 
intended to be included within the scope of the present 
invention as setforth by the claims. 
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